Optimization of oxidative folding methods for cysteine‐rich peptides: a study of conotoxins containing three disulfide bridges

The oxidative folding of small, cysteine‐rich peptides to selectively achieve the native disulfide bond connectivities is critical for discovery and structure‐function studies of many bioactive peptides. As the propensity to acquire the native conformation greatly depends on the peptide sequence, numerous empirical oxidation methods are employed. The context‐dependent optimization of these methods has thus far precluded a generalized oxidative folding protocol, in particular for peptides containing more than two disulfides. Herein, we compare the efficacy of optimized solution‐phase and polymer‐supported oxidation methods using three disulfide‐bridged conotoxins, namely µ‐SIIIA, µ‐KIIIA and ω‐GVIA. The use of diselenide bridges as proxies for disulfide bridges is also evaluated. We propose the ClearOx‐assisted oxidation of selenopeptides as a fairly generalized oxidative folding protocol. Copyright © 2010 European Peptide Society and John Wiley & Sons, Ltd.     

Principles of carbopeptoid folding: a molecular dynamics simulation study.

The conformational spaces of five oligomers of tetrahydrofuran-based carbopeptoids in chloroform and dimethyl sulfoxide were investigated through nine molecular dynamics simulations. Prompted by nuclear magnetic resonance experiments that indicated various stable folds for some but not all of these carbopeptoids, their folding behaviour was investigated as a function of stereochemistry, chain length and solvent. The conformational distributions of these molecules were analysed in terms of occurrence of hydrogen bonds, backbone torsional-angle distributions, conformational clustering and solute configurational entropy. While a cis-linkage across the tetrahydrofuran ring favours right-handed helical structures, a trans-linkage results in a larger conformational variability. Intra-solute hydrogen bonding is reduced with increasing chain length and with increasing solvent polarity. Solute configurational entropies confirm the picture obtained: they are smaller for cis- than for trans-linked peptides, for chloroform than for dimethyl sulfoxide as solvent and for shorter peptide chains. The simulations provide an atomic picture of molecular conformational variability that is consistent with the available experimental data.     

β‐Hairpin folding and stability: molecular dynamics simulations of designed peptides in aqueous solution

The structural properties of a 10‐residue and a 15‐residue peptide in aqueous solution were investigated by molecular dynamics simulation. The two designed peptides, SYINSDGTWT and SESYINSDGTWTVTE, had been studied previously by NMR at 278 K and the resulting model structures were classified as 3:5 β‐hairpins with a type I + G1 β‐bulge turn. In simulations at 278 K, starting from the NMR model structure, the 3:5 β‐hairpin conformers proved to be stable over the time period evaluated (30 ns). Starting from an extended conformation, simulations of the decapeptide at 278 K, 323 K and 353 K were also performed to study folding. Over the relatively short time scales explored (30 ns at 278 K and 323 K, 56 ns at 353 K), folding to the 3:5 β‐hairpin could only be observed at 353 K. At this temperature, the collapse to β‐hairpin‐like conformations is very fast. The conformational space accessible to the peptide is entirely dominated by loop structures with different degrees of β‐hairpin character. The transitions between different types of ordered loops and β‐hairpins occur through two unstructured loop conformations stabilized by a single side‐chain interaction between Tyr2 and Trp9, which facilitates the changes of the hydrogen‐bond register. In agreement with previous experimental results, β‐hairpin formation is initially driven by the bending propensity of the turn segment. Nevertheless, the fine organization of the turn region appears to be a late event in the folding process. Copyright © 2004 European Peptide Society and John Wiley & Sons, Ltd.     

One membrane protein,two structures and six environments: a comparative molecular dynamics simulation study of the bacterial outer membrane protein PagP

PagP is a bacterial outer membrane protein consisting of an 8 stranded transmembrane β-barrel and an N-terminal α-helix. It is an enzyme which catalyses transfer of a palmitoyl chain from a phospholipid to lipid A. Molecular dynamics simulations have been used to compare the dynamic behaviour in simulations starting from two different structures (X-ray vs. NMR) and in six different environments (detergent micelles formed by dodecyl phosphocholine and by octyl glucoside, vs. four species of phospholipid bilayer). Analysis of interactions between the protein and its environment reveals the role played by the N-terminal α-helix, which interacts with the lipid headgroups to lock the PagP molecule into the bilayer. The PagP β-barrel adopts a tilted orientation in lipid bilayers, facilitating access of lipid tails into the mouth of the central binding pocket. In simulations starting from the X-ray structure in lipid bilayer, the L1 and L2 loops move towards one another, leading to the formation of a putative active site by residues H33, D76 and S77 coming closer together.     

The N‐terminal of annexin A1 as a secondary membrane binding site: A molecular dynamics study

Annexin A1 has been shown to cause membrane aggregation and fusion, yet the mechanism of these activities is not clearly understood. In this work, molecular dynamics simulations were performed on monomeric annexin A1 positioned between two negatively charged monolayers using AMBER's all atom force field to gain insight into the mechanism of fusion. Each phospolipid monolayer was made up of 180 DOPC molecules and 45 DOPG molecules to achieve a 4:1 ratio. The space between the two monolayers was explicitly solvated using TIP3P waters in a rectilinear box. The constructed setup contained up to 0.14 million atoms. Application of periodic boundary conditions to the simulation setup gave the desired effect of two continuous membrane bilayers. Nonbonded interactions were calculated between the N‐terminal residues and the bottom layer of phospholipids, which displayed a strong attraction of K26 and K29 to the lipid head‐groups. The side‐chains of these two residues were observed to orient themselves in close proximity (～3.5 Å) with the polar head‐groups of the phospholipids. Proteins 2014; 82:2936–2942. © 2014 Wiley Periodicals, Inc.     

Transient stability of the helical pattern of region F19-L22 of the N-terminal domain of p53: a molecular dynamics simulation study

Two molecular dynamics simulations of the region E17-N29 of p53 (p53(17-29)) at different temperatures were performed for a total time of 0.2 micros, to study the conformational landscape of this region. Previous studies have suggested that this region displays different structural motifs, such as helix of a double beta-turn, and that its secondary structure might be transiently stable. Interestingly, in this study it was found that the region F19-L25, and particularly its fragment F19-L22, display a stable, transient helical pattern at sub-microsecond periods. The region F19-L22, which contains one of the most important residues needed for the interaction of p53 with MDM2, seems to be formed and stabilized by the existence of one hydrophobic and one aromatic cluster. The main function of these clusters is to help their surrounding area to desolvate, to allow the hydrogen bond network, therefore favoring the formation of a stable helix. This preliminary study would be useful for a better understanding of the structure and function of the N-terminal domain of p53 and its implications for the control of different types of cancer.     

Conformational dynamics of a short antigenic peptide in its free and antibody bound forms gives insight into the role of β‐turns in peptide immunogenicity

Earlier immunological experiments with a synthetic 36‐residue peptide (75‐110) from Influenza hemagglutinin have been shown to elicit anti‐peptide antibodies (Ab) which could cross‐react with the parent protein. In this article, we have studied the conformational features of a short antigenic (Ag) peptide (⁹⁸YPYDVPDYASLRS¹¹⁰) from Influenza hemagglutinin in its free and antibody (Ab) bound forms with molecular dynamics simulations using GROMACS package and OPLS‐AA/L all‐atom force field at two different temperatures (293 K and 310 K). Multiple simulations for the free Ag peptide show sampling of ordered conformations and suggest different conformational preferences of the peptide at the two temperatures. The free Ag samples a conformation crucial for Ab binding (β‐turn formed by “DYAS” sequence) with greater preference at 310 K while, it samples a native‐like conformation with relatively greater propensity at 293 K. The sequence “DYAS” samples β‐turn conformation with greater propensity at 310 K as part of the hemagglutinin protein also. The bound Ag too samples the β‐turn involving “DYAS” sequence and in addition it also samples a β‐turn formed by the sequence “YPYD” at its N‐terminus, which seems to be induced upon binding to the Ab. Further, the bound Ag displays conformational flexibility at both 293 K and 310 K, particularly at terminal residues. The implications of these results for peptide immunogenicity and Ag–Ab recognition are discussed. Proteins 2015; 83:1352–1367. © 2015 Wiley Periodicals, Inc.     

On the stability of different experimental dimeric structures of the SL1 sequence from the genomic RNA of HIV-1 in solution: a molecular dynamics simulation and electrophoresis study

SL1 is a stem-loop RNA sequence from the genome of HIV-1 thought to be the initiation site for the dimerization of the retroviral genomic RNA. The aim of this study is to check the stability in solution of different experimental dimeric structures available in the literature. Two kinds of dimer have been evidenced: an extended duplex looking like a double helix with two internal bulges and a kissing complex in which the monomers with a stem/loop conformation are linked by intermolecular loop-loop interactions. Two divergent experimental structures of the kissing complex from the Lai isolate are reported in the literature, one obtained from NMR (Mujeeb et al., Nature Structural Biology, 1998, Vol. 5, pp. 432-436) and the other one from x-ray crystallography (Ennifar et al., Nature Structural Biology, 2001, Vol. 8, pp. 1064-1068). A crystallographic structure of the Mal isolate was also reported (Ennifar et al., Nature Structure Biology, 2001, Vol. 8, pp. 1064-1068). Concerning the extended duplex, a NMR structure is available for Lai (Girard et al., Journal of Biomolecular Structure and Dynamics, 1999, Vol. 16, pp. 1145-1157) and a crystallographic structure for Mal (Ennifar et al., Structure, 1999, Vol. 7, pp. 1439-1449). Using a molecular dynamics technique, all these experimental structures have been simulated in solution with explicit water and counterions. We show that both extended duplex structures are stable. On the contrary, the crystallographic structures of the Lai and Mal kissing complexes are rapidly destabilized in aqueous environment. Finally, the NMR structure of the Lai loop-loop kissing complex remains globally stable over a 20 ns MD simulation, although large rearrangements occur at the level of the stem/loop junctions that are flexible, as shown from free energy calculations. These results are compared to electrophoresis experiments on dimer formation.     

Exploring the dynamic information content of a protein NMR structure: Comparison of a molecular dynamics simulation with the NMR and X‐ray structures of Escherichia coli ribonuclease HI

The multiconformer nature of solution nuclear magnetic resonance (NMR) structures of proteins results from the effects of intramolecular dynamics, spin diffusion and an uneven distribution of structural restraints throughout the molecule. A delineation of the former from the latter two contributions is attempted in this work for an ensemble of 15 NMR structures of the protein Escherichia coli ribonuclease HI (RNase HI). Exploration of the dynamic information content of the NMR ensemble is carried out through correlation with data from two crystal structures and a 1.7‐ns molecular dynamics (MD) trajectory of RNase HI in explicit solvent. Assessment of the consistency of the crystal and mean MD structures with nuclear Overhauser effect (NOE) data showed that the NMR ensemble is overall more compatible with the high‐resolution (1.48 Å) crystal structure than with either the lower‐resolution (2.05 Å) crystal structure or the MD simulation. Furthermore, the NMR ensemble is found to span more conformational space than the MD simulation for both the backbone and the sidechains of RNase HI. Nonetheless, the backbone conformational variability of both the NMR ensemble and the simulation is especially consistent with NMR relaxation measurements of two loop regions that are putative sites of substrate recognition. Plausible side‐chain dynamic information is extracted from the NMR ensemble on the basis of (i) rotamericity and syn‐pentane character of variable torsion angles, (ii) comparison of the magnitude of atomic mean‐square fluctuations (msf) with those deduced from crystallographic thermal factors, and (iii) comparison of torsion angle conformational behavior in the NMR ensemble and the simulation. Several heterogeneous torsion angles, while adopting non‐rotameric/syn‐pentane conformations in the NMR ensemble, exist in a unique conformation in the simulation and display low X‐ray thermal factors. These torsions are identified as sites whose variability is likely to be an artifact of the NMR structure determination procedure. A number of other torsions show a close correspondence between the conformations sampled in the NMR and MD ensembles, as well as significant correlations among crystallographic thermal factors and atomic msf calculated from the NMR ensemble and the simulation. These results indicate that a significant amount of dynamic information is contained in the NMR ensemble. The relevance of the present findings for the biological function of RNase HI, protein recognition studies, and previous investigations of the motional content of protein NMR structures are discussed. Proteins 1999;36:87–110. © 1999 Wiley‐Liss, Inc.     

Study of the interactions between a proline‐rich protein and a flavan‐3‐ol by NMR: Residual structures in the natively unfolded protein provides anchorage points for the ligands

Astringency is one of the major organoleptic properties of food and beverages that are made from plants, such as tea, chocolate, beer, or red wine. This sensation is thought to be due to interactions between tannins and salivary proline‐rich proteins, which are natively unfolded proteins. A human salivary proline‐rich protein, namely IB‐5, was produced by the recombinant method. Its interactions with a model tannin, epigallocatechin gallate (EGCG), the major flavan‐3‐ol in green tea, were studied here. Circular dichroism experiments showed that IB‐5 presents residual structures (PPII helices) when the ionic strength is close to that in saliva. In the presence of these residual structures, IB‐5 undergoes an increase in structural content upon binding to EGCG. NMR data corroborated the presence of preformed structural elements within the protein prior to binding and a partial assignment was proposed, showing partial structuration. TOCSY experiments showed that amino acids that are involved in PPII helices are more likely to interact with EGCG than those in random coil regions, as if they were anchorage points for the ligand. The signal from IB‐5 in the DOSY NMR spectrum revealed an increase in polydispersity upon addition of EGCG while the mean hydrodynamic radius remained unchanged. This strongly suggests the formation of IB‐5/EGCG aggregates. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 745–756, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Cyclic pentapeptides of chiral sequence DLDDL as scaffold for antagonism of G‐protein coupled receptors: Synthesis,activity and conformational analysis by NMR and molecular dynamics of ITF 1565 a substance P inhibitor

Under the hypotheses of a structurally related binding site for antagonists of G‐protein coupled receptors and the ability of cyclic pentapeptides of chiral sequence D ¹L ²D ³D ⁴L ⁵ to form rigid structures with which probe the pharmacophoric specificity of these receptors, inhibitors of substance P were designed based on available structure–activity relationships. ITF 1565, cyclo[D ‐Trp¹‐Pro²‐D ‐Lys³‐D ‐Trp⁴‐Phe⁵], antagonized substance P activity mediated by type 1 neurokinin receptor (NK1) whereas it acted weakly against NK2 and did not inhibit endothelin at all. The preferential conformation of the peptide was obtained from nmr spectroscopy and computer calculations, and shown to contain the same βII‐turn and γ′‐turn found in other cyclic pentapeptides with the same chiral sequence. The structure of the peptide was compared with that of the β‐D ‐glucose molecule that has been proposed as a semirigid scaffold for antagonists of G‐protein coupled receptors. The γ′‐turn of the cyclic peptide superimposed well with β‐D ‐glucose in the chair conformation. Furthermore, when the side chains were considered, the aromatic groups of the two molecules were found to generally overlap. These results support the view of G‐protein coupled receptors as possessing structurally similar binding sites for antagonists and suggest that cyclic pentapeptides of chiral sequence D ¹L ²D ³D ⁴L ⁵ may be useful as semirigid scaffolds for the design of antagonists of this family of receptors. © 1999 John Wiley & Sons, Inc. Biopoly 50: 211–219, 1999     

Insights into the binding specificity of wild type and mutated wheat germ agglutinin towards Neu5Acα(2‐3)Gal: a study by in silico mutations and molecular dynamics simulations

Wheat germ agglutinin (WGA) is a plant lectin, which specifically recognizes the sugars NeuNAc and GlcNAc. Mutated WGA with enhanced binding specificity can be used as biomarkers for cancer. In silico mutations are performed at the active site of WGA to enhance the binding specificity towards sialylglycans, and molecular dynamics simulations of 20 ns are carried out for wild type and mutated WGAs (WGA1, WGA2, and WGA3) in complex with sialylgalactose to examine the change in binding specificity. MD simulations reveal the change in binding specificity of wild type and mutated WGAs towards sialylgalactose and bound conformational flexibility of sialylgalactose. The mutated polar amino acid residues Asn114 (S114N), Lys118 (G118K), and Arg118 (G118R) make direct and water mediated hydrogen bonds and hydrophobic interactions with sialylgalactose. An analysis of possible hydrogen bonds, hydrophobic interactions, total pair wise interaction energy between active site residues and sialylgalactose and MM‐PBSA free energy calculation reveals the plausible binding modes and the role of water in stabilizing different binding modes. An interesting observation is that the binding specificity of mutated WGAs (cyborg lectin) towards sialylgalactose is found to be higher in double point mutation (WGA3). One of the substituted residues Arg118 plays a crucial role in sugar binding. Based on the interactions and energy calculations, it is concluded that the order of binding specificity of WGAs towards sialylgalactose is WGA3 > WGA1 > WGA2 > WGA. On comparing with the wild type, double point mutated WGA (WGA3) exhibits increased specificity towards sialylgalactose, and thus, it can be effectively used in targeted drug delivery and as biological cell marker in cancer therapeutics. Copyright © 2014 John Wiley & Sons, Ltd.