Use of restrained molecular dynamics in water to determine three-dimensional protein structure: prediction of the three-dimensional structure of Ecballium elaterium trypsin inhibitor II

Synthetic Multivalent Disulfide-Constrained Peptide Agonists Potentiate Wnt1/β-Catenin Signaling via LRP6 Coreceptor Clustering

Wnt ligands are critical for tissue homeostasis and form a complex with LRP6 and frizzled coreceptors to initiate Wnt/β-catenin signaling. Yet, how different Wnts achieve various levels of signaling activation through distinct domains on LRP6 remains elusive. Developing tool ligands that target individual LRP6 domains could help elucidate the mechanism of Wnt signaling regulation and uncover pharmacological approaches for pathway modulation. We employed directed evolution of a disulfide constrained peptide (DCP) to identify molecules that bind to the third β-propeller domain of LRP6. The DCPs antagonize Wnt3a while sparing Wnt1 signaling. Using PEG linkers with different geometries, we converted the Wnt3a antagonist DCPs to multivalent molecules that potentiated Wnt1 signaling by clustering the LRP6 coreceptor. The mechanism of potentiation is unique as it occurred only in the presence of extracellular secreted Wnt1 ligand. While all DCPs recognized a similar binding interface on LRP6, they displayed different spatial orientations that influenced their cellular activities. Moreover, structural analyses revealed that the DCPs exhibited new folds that were distinct from the parent DCP framework they were evolved from. The multivalent ligand design principles highlighted in this study provide a path for developing peptide agonists that modulate different branches of cellular Wnt signaling.

Heterochiral Knottin Protein: Folding and Solution Structure

Homochirality is a general feature of biological macromolecules, and Nature includes few examples of heterochiral proteins. Herein, we report on the design, chemical synthesis, and structural characterization of heterochiral proteins possessing loops of amino acids of chirality opposite to that of the rest of a protein scaffold. Using the protein Ecballium elaterium trypsin inhibitor II, we discover that selective β-alanine substitution favors the efficient folding of our heterochiral constructs. Solution nuclear magnetic resonance spectroscopy of one such heterochiral protein reveals a homogeneous global fold. Additionally, steered molecular dynamics simulation indicate β-alanine reduces the free energy required to fold the protein. We also find these heterochiral proteins to be more resistant to proteolysis than homochiral l-proteins. This work informs the design of heterochiral protein architectures containing stretches of both d- and l-amino acids.

Solid phase protein chemical synthesis

The chemical synthesis of peptides or small proteins is often an important step in many research projects and has stimulated the development of numerous chemical methodologies. The aim of this review is to give a substantial overview of the solid phase methods developed for the production or purification of polypeptides. The solid phase peptide synthesis (SPPS) technique has facilitated considerably the access to short peptides (<50 amino acids). However, its limitations for producing large homogeneous peptides have stimulated the development of solid phase covalent or non-covalent capture purification methods. The power of the native chemical ligation (NCL) reaction for protein synthesis in aqueous solution has also been adapted to the solid phase by the combination of novel linker technologies, cysteine protection strategies and thioester or N,S-acyl shift thioester surrogate chemistries. This review details pioneering studies and the most recent publications related to the solid phase chemical synthesis of large peptides and proteins.

Selenopeptide analogs of EETI-II retain potent trypsin inhibitory activities

Three-disulfide-bridged Ecballium elaterium trypsin inhibitor II (EETI-II) is a 28-residue peptide that belongs to the squash family of canonical trypsin inhibitors. Herein, we report synthesis and biological activity of three EETI-II analogs. In each of analog, a pair of cysteine residues forming a native disulfide bridge was individually replaced by a pair of selenocysteine residues. All selenopeptide analogs were chemically synthesized using the Fmoc protocol and subsequently folded in the presence of oxidized and reduced glutathione. The analogs containing a diselenide bridge displayed association constants with trypsin that ranged from 2.6 x 10(9) to 5.1 x 10(9)  [M(-1) ]. Our results suggest that the selenopeptide analogs retained low nanomolar inhibitory potencies, and only the diselenide bridge adjacent to the inhibitory binding loop weakened the interactions with trypsin by approximately fivefold. We discuss these findings in the context of a broader use of selenopeptide analogs as proxies to study cysteine-rich peptides.

Conopeptide characterization and classifications: an analysis using ConoServer

Cone snails are carnivorous marine gastropods that have evolved potent venoms to capture their prey. These venoms comprise a rich and diverse cocktail of peptide toxins, or conopeptides, whose high diversity has arisen from an efficient hypermutation mechanism, combined with a high frequency of post-translational modifications. Conopeptides bind with high specificity to distinct membrane receptors, ion channels, and transporters of the central and muscular nervous system. As well as serving their natural function in prey capture, conopeptides have been utilized as versatile tools in neuroscience and have proven valuable as drug leads that target the nervous system in humans. This paper examines current knowledge on conopeptide sequences based on an analysis of gene and peptide sequences in ConoServer (http://www.conoserver.org), a specialized database of conopeptide sequences and three-dimensional structures. We describe updates to the content and organization of ConoServer and discuss correlations between gene superfamilies, cysteine frameworks, pharmacological families targeted by conopeptides, and the phylogeny, habitat, and diet of cone snails. The study identifies gaps in current knowledge of conopeptides and points to potential directions for future research.

Structure and folding of disulfide-rich miniproteins: insights from molecular dynamics simulations and MM-PBSA free energy calculations

The fold of small disulfide-rich proteins largely relies on two or more disulfide bridges that are main components of the hydrophobic core. Because of the small size of these proteins and their high cystine content, the cysteine connectivity has been difficult to ascertain in some cases, leading to uncertainties and debates in the literature. Here, we use molecular dynamics simulations and MM-PBSA free energy calculations to compare similar folds with different disulfide pairings in two disulfide-rich miniprotein families, namely the knottins and the short-chain scorpion toxins, for which the connectivity has been discussed. We first show that the MM-PBSA approach is able to discriminate the correct knotted topology of knottins from the laddered one. Interestingly, a comparison of the free energy components for kalata B1 and MCoTI-II suggests that cyclotides and squash inhibitors, although sharing the same scaffold, are stabilized through different interactions. Application to short-chain scorpion toxins suggests that the conventional cysteine pairing found in many homologous toxins is significantly more stable than the unconventional pairing reported for maurotoxin and for spinoxin. This would mean that native maurotoxin and spinoxin are not at the lowest free energy minimum and might result from kinetically rather than thermodynamically driven oxidative folding processes. For both knottins and toxins, the correct or conventional disulfide connectivities provide lower flexibilities and smaller deviations from the initial conformations. Overall, our work suggests that molecular dynamics simulations and the MM-PBSA approach to estimate free energies are useful tools to analyze and compare disulfide bridge connectivities in miniproteins.

Knottin cyclization: impact on structure and dynamics

Background

Present in various species, the knottins (also referred to as inhibitor cystine knots) constitute a group of extremely stable miniproteins with a plethora of biological activities. Owing to their small size and their high stability, knottins are considered as excellent leads or scaffolds in drug design. Two knottin families contain macrocyclic compounds, namely the cyclotides and the squash inhibitors. The cyclotide family nearly exclusively contains head-to-tail cyclized members. On the other hand, the squash family predominantly contains linear members. Head-to-tail cyclization is intuitively expected to improve bioactivities by increasing stability and lowering flexibility as well as sensitivity to proteolytic attack.

Results

In this paper, we report data on solution structure, thermal stability, and flexibility as inferred from NMR experiments and molecular dynamics simulations of a linear squash inhibitor EETI-II, a circular squash inhibitor MCoTI-II, and a linear analog lin-MCoTI. Strikingly, the head-to-tail linker in cyclic MCoTI-II is by far the most flexible region of all three compounds. Moreover, we show that cyclic and linear squash inhibitors do not display large differences in structure or flexibility in standard conditions, raising the question as to why few squash inhibitors have evolved into cyclic compounds. The simulations revealed however that the cyclization increases resistance to high temperatures by limiting structure unfolding.

Conclusion

In this work, we show that, in contrast to what could have been intuitively expected, cyclization of squash inhibitors does not provide clear stability or flexibility modification. Overall, our results suggest that, for squash inhibitors in standard conditions, the circularization impact might come from incorporation of an additional loop sequence, that can contribute to the miniprotein specificity and affinity, rather than from an increase in conformational rigidity or protein stability. Unfolding simulations showed however that cyclization is a stabilizing factor in strongly denaturing conditions. This information should be useful if one wants to use the squash inhibitor scaffold in drug design.

Barnase Fusion as a Tool to Determine the Crystal Structure of the Small Disulfide-rich Protein McoEeTI

We present a fusion system suited to determine the crystal structure of small disulfide-rich proteins. McoEeTI, a hybrid inhibitor cystine knot microprotein, was produced as a soluble fusion to a catalytically inactive variant of the RNAse barnase in Escherichia coli. Functioning as a versatile tag, barnase facilitated purification, crystallization and high-resolution structure determination. Flexibility of the linker region allows for different relative orientations of barnase and the fusion partner in two crystallographically independent molecules and may thereby facilitate crystal packing. Nevertheless, the linker region is well ordered in both molecules. This system may prove more generally useful to determine the crystal structure of peptides and small proteins.

A fusion protein system for the recombinant production of short disulfide bond rich cystine knot peptides using barnase as a purification handle

The inhibitor cystine knot (ICK) structural motif has been found in several small proteins and peptides from plants, insects, marine molluscs, and also in human. It is defined by a triple beta-sheet that is held together by three intramolecular disulfide bonds built by six conserved cysteine residues that generate a highly rigid and stable fold. We describe a procedure for the production of ICK peptides with correct disulfide bond connectivities via expression in Escherichia coli as fusion proteins with an enzymatically inactive variant of the Bacillus amyloliquefaciens RNAse barnase. Barnase directs the fused peptide to the culture medium and the fusion protein can be isolated by combined cation exchange/reverse-phase chromatography. The ICK peptides are released from the barnase expression and purification handle either by cyanogen bromide or by protease cleavage to give pure and correctly folded cystine knot peptides.

The KNOTTIN website and database: a new information system dedicated to the knottin scaffold

The KNOTTIN website and database organize information about knottins or inhibitor cystine knots, small disulfide-rich proteins with a knotted topology. Thanks to their small size and high stability, knottins provide appealing scaffolds for protein engineering and drug design. Static pages present the main historical and recent results about knottin discoveries, sequences, structures, folding, functions, applications and bibliography. Database searches provide dynamically generated tabular reports or sequence alignments for knottin three-dimensional structures or sequences. BLAST/HMM searches are also available. A simple nomenclature, based on loop lengths between cysteines, is proposed and is complemented by a uniform numbering scheme. This standardization is applied to all knottin structures in the database, facilitating comparisons. Renumbered and structurally fitted knottin PDB files are available for download. The standardized numbering is used for automatic drawing of two-dimensional Colliers de Perles. The KNOTTIN website and database are available at http://knottin.cbs.cnrs.fr and http://knottin.com.

Structural analysis of free and enzyme-bound amaranth alpha-amylase inhibitor: classification within the knottin fold superfamily and analysis of its functional flexibility

The three-dimensional structure of the amaranth alpha-amylase inhibitor (AAI) adopts a knottin fold of abcabc topology. Upon binding to alpha-amylase, it adopts a more compact conformation characterized by an increased number of intramolecular hydrogen bonds, a decreased volume and in addition a trans to cis isomerization of Pro20. A systematic analysis of the 3-D structural databanks revealed that similar proteins and domains share with AAI the characteristic presence of proline residues, many of which are in a cis backbone conformation. As these proteins fulfil a variety of functional roles and are expressed in very different organisms, we conclude that the structure of the knottin fold, including the propensity of the cis bond, are the result of convergent evolution.

The cystine knot motif in toxins and implications for drug design

The cystine knot structural motif is present in peptides and proteins from a variety of species, including fungi, plants, marine molluscs, insects and spiders. It comprises an embedded ring formed by two disulfide bonds and their connecting backbone segments which is threaded by a third disulfide bond. It is invariably associated with nearby beta-sheet structure and appears to be a highly efficient motif for structure stabilization. Because of this stability it makes an ideal framework for molecular engineering applications. In this review we summarize the main structural features of the cystine knot motif, focussing on toxin molecules containing either the inhibitor cystine knot or the cyclic cystine knot. Peptides containing these motifs are 26-48 residues long and include ion channel blockers, haemolytic agents, as well as molecules having antiviral and antibacterial activities. The stability of peptide toxins containing the cystine knot motif, their range of bioactivities and their unique structural scaffold can be harnessed for molecular engineering applications and in drug design. Applications of cystine knot molecules for the treatment of pain, and their potential use in antiviral and antibacterial applications are described.

Solution structure of the major alpha-amylase inhibitor of the crop plant amaranth

alpha-Amylase inhibitor (AAI), a 32-residue miniprotein from the Mexican crop plant amaranth (Amaranthus hypochondriacus), is the smallest known alpha-amylase inhibitor and is specific for insect alpha-amylases (Chagolla-Lopez, A., Blanco-Labra, A., Patthy, A., Sanchez, R., and Pongor, S. (1994) J. Biol. Chem. 269, 23675-23680). Its disulfide topology was confirmed by Edman degradation, and its three-dimensional solution structure was determined by two-dimensional 1H NMR spectroscopy at 500 MHz. Structural constraints (consisting of 348 nuclear Overhauser effect interproton distances, 8 backbone dihedral constraints, and 9 disulfide distance constraints) were used as an input to the X-PLOR program for simulated annealing and energy minimization calculations. The final set of 10 structures had a mean pairwise root mean square deviation of 0.32 A for the backbone atoms and 1.04 A for all heavy atoms. The structure of AAI consists of a short triple-stranded beta-sheet stabilized by three disulfide bonds, forming a typical knottin or inhibitor cystine knot fold found in miniproteins, which binds various macromolecular ligands. When the first intercystine segment of AAI (sequence IPKWNR) was inserted into a homologous position of the spider toxin Huwentoxin I, the resulting chimera showed a significant inhibitory activity, suggesting that this segment takes part in enzyme binding.

Molecular dynamics study of substance P peptides partitioned in a sodium dodecylsulfate micelle

Two neuropeptides, substance P (SP) and SP-tyrosine-8 (SP-Y8), have been studied by molecular dynamics (MD) simulation in an explicit sodium dodecylsulfate (SDS) micelle. Initially, distance restraints derived from NMR nuclear Overhauser enhancements (NOE) were incorporated in the restrained MD (RMD) during the equilibration stage of the simulation. It was shown that when SP-Y8 was initially placed in an insertion (perpendicular) configuration, the peptide equilibrated to a surface-bound (parallel) configuration in approximately 450 ps. After equilibration, the conformation and orientation of the peptides, the solvation of both the backbone and the side chain of the residues, hydrogen bonding, and the dynamics of the peptides were analyzed from trajectories obtained from the RMD or the subsequent free MD (where the NOE restraints were removed). These analyses showed that the peptide backbones of all residues are either solvated by water or are hydrogen-bonded. This is seen to be an important factor against the insertion mode of interaction. Most of the interactions come from the hydrophobic interaction between the side chains of Lys-3, Pro-4, Phe-7, Phe-8, Leu-10, and Met-11 for SP, from Lys-3, Phe-7, Leu-10, and Met-11 in SP-Y8, and the micellar interior. Significant interactions, electrostatic and hydrogen bonding, between the N-terminal residues, Arg-Pro-Lys, and the micellar headgroups were observed. These latter interactions served to affect both the structure and, especially, the flexibility, of the N-terminus. The results from simulation of the same peptides in a water/CCl4 biphasic cell were compared with the results of the present study, and the validity of using the biphasic system as an approximation for peptide-micelle or peptide-bilayer systems is discussed.

Solution structure of R-elafin, a specific inhibitor of elastase

The solution structure of r-elafin, a specific elastase inhibitor, has been determined using NMR spectroscopy. Characterized by a flat core and a flexible N-terminal extremity, the three-dimensional structure is formed by a central twisted beta-hairpin accompanied by two external segments linked by the proteinase binding loop. A cluster of three disulfide bridges connects the external segments to the central beta-sheet and a single fourth disulfide bridge links the binding loop to the central beta-turn. The same spatial distribution of disulfide bridges can be observed in both domains of the secretory leukocyte protease inhibitor (SLPI), another elastase inhibitor. The structural homology between r-elafin and the C-terminal domain of SLPI confirms the former as a second member of the chelonianin family of proteinase inhibitors. Based on the homology between the two proteins and recent results obtained for elastase binding mutants of the bovine pancreatic trypsin inhibitor (BPTI), we define the segment 22 to 27 as the binding loop of elafin, with the scissile peptide bond between Ala24 and Met25. In our solution structures, this loop is extended and solvent-exposed, and exhibits a large degree of flexibility. This mobility, already observed for the binding loop in other protease inhibitors in solution, might be an important feature for the interaction with the corresponding protease.

Role of a nonnative interaction in the folding of the protein G B1 domain as inferred from the conformational analysis of the alpha-helix fragment

BACKGROUND

The role of local interactions in protein folding and stability can be investigated by the conformational analysis of protein fragments. The hydrophobic staple and Schellman motifs have been described at the N and C terminus, respectively, of protein alpha-helices. These motifs are characterized by an interaction between two hydrophobic residues, one outside the helix and one within the helix, and their importance for helix stability has been analyzed in model peptides. In the alpha-helix of the protein G B1 domain, only the Schellman motif is formed--the hydrophobic staple motif is absent despite the favourable sequence pattern. We have experimentally analyzed the solution conformation of the 19-41 fragment of protein G. This peptide comprises the helical residues and contains both the hydrophobic staple and Schellman motif sequences.

RESULTS

In the isolated peptide in water, the hydrophobic staple motif is formed and stabilizes the helical structure as compared with a shorter peptide lacking it, but the Schellman motif is not formed. In 30% aqueous TFE, the helix is more stable than in pure water and both motifs are formed.

CONCLUSIONS

The results suggest that the importance of each motif for the folding and stability of protein G is different. The nonnative hydrophobic staple interaction can help to nucleate the helix at the beginning of folding but has later to be disrupted. The Schellman motif, while not providing enough energy for substantial helix stabilization in the unfolded state, could be important for determining the local fold of the sequence in the context of the rest of the protein.

Synthesis and conformational analysis by 1H NMR and restrained molecular dynamics simulations of the cyclic decapeptide [Ser-Tyr-Ser-Met-Glu-His-Phe-Arg-Trp-Gly]

The design of enzyme mimics with therapeutic and industrial applications has interested both experimental and computational chemists for several decades. Recent advances in the computational methodology of restrained molecular dynamics, used in conjunction with data obtained from two-dimensional 1H NMR spectroscopy, make it a promising method to study peptide and protein structure and function. Several issues, however, need to be addressed in order to assess the validity of this method for its explanatory and predictive value. Among the issues addressed in this study are: the accuracy and generizability of the GROMOS peptide molecular mechanics force field; the effect of inclusion of solvent on the simulations; and the effect of different types of restraining algorithms on the computational results. The decapeptide Ser-Tyr-Ser-Met-Glu-His-Phe-Arg-Trp-Gly, which corresponds to the sequence of ACTH1-10, has been synthesized, cyclized, and studied by two-dimensional 1H NMR spectroscopy. Restrained molecular dynamics (RMD) and time-averaged restrained molecular dynamics (TARMD) simulations were carried out on four different distance-geometry starting structures in order to determine and contrast the behavior of cyclic ACTH1-10 in vacuum and in solution. For the RMD simulations, the structures did not fit the NOE data well, even at high values of the restraining potential. The TARMD simulation method, however, was able to give structures that fit the NOE data at high values of the restraining potential. In both cases, inclusion of explicit solvent molecules in the simulation had little effect on the quality of the fit, although it was found to dampen the motion of the cyclic peptide. For both simulation techniques, the number and size of the NOE violations increased as the restraining potential approached zero. This is due, presumably, to inadequacies in the force field. Additional TARMD vacuum-phase simulations, run with a larger memory length or with a larger sampling size (16 additional distance-geometry structures), yielded no significantly different results. The computed data were then analyzed to help explain the sparse NOE data and poor chymotryptic activity of the cyclic peptide. Cyclic ACTH1-10, which contains the functional moieties of the catalytic triad of chymotrypsin, was evaluated as a potential mimic of chymotrypsin by measurement of the rate of hydrolysis of esters of L- and D-phenylalanine. The poor rate of hydrolysis is attributed to the flexibility of the decapeptide, the motion of the side chains, which result in the absence of long-range NOEs, the small size of the macrocycle relative to that of the substrate, and the inappropriate orientation of the Gly, His, and Ser residues. The results demonstrate the utility of this method in computer-aided molecular design of cyclic peptides and suggest structural modifications for future work based on a larger and more rigid peptide framework.

Structure and multiple conformations of the kunitz-type domain from human type VI collagen alpha3(VI) chain in solution

BACKGROUND

The Kunitz-type inhibitor motif is found at the C terminus of the human collagen alpha3(VI) chain. This 76-residue module (domain C5) was prepared in recombinant form and showed high stability against proteases; however, it lacked any inhibitory activity against trypsin, thrombin, kallikrein and several other proteases. We have undertaken the determination of the three-dimensional (3D) structure of domain C5 in solution, by nuclear magnetic resonance (NMR), in order to establish the structural basis for the properties of this protein.

RESULTS

The 7 N-terminal and 12 C-terminal residues of domain C5 are disordered in the solution structure. The 55-residue core, which shows high homology to bovine pancreatic trypsin inhibitor, retains the characteristic fold of all members of the Kunitz-type inhibitor family. 24 residues of this main structural body show more than one resonance, symptomatic of multiple conformations slowly exchanging on the NMR time scale. In addition, significant proton chemical exchange line broadening is observed for residues in the vicinity of the disulfide bridge between residues 20 and 44: this indicates interconversion, on the micro- to millisecond time scale, between multiple conformations.

CONCLUSION

The NMR study demonstrates that domain C5 is a highly dynamic molecule at temperatures studied (between 10 and 30 degrees C). Indeed, some 44% of the main body structure of C5 showed multiple conformations. The existence of multiple conformations was not necessarily expected in view of the conformational constraints imposed by the 3D structure of proteins as rigid as C5; it should therefore be considered in the interpretation of its structural and dynamical properties. The accessibility of the inhibitory binding loop (Gly18 [P4] to Leu25 [P4']) should be relatively unaffected by this conformational exchange and thus would not explain the unusual specificity of C5. Most serine proteinase inhibitors that, like C5, have an arginine at the P1 position inhibit trypsin; the lack of trypsin inhibition of C5 must therefore arise from the amino-acid side-chain composition of the adjoining positions in the binding loop.

Folding of the squash trypsin inhibitor EETI II. Evidence of native and non-native local structural preferences in a linear analogue

A peptide, corresponding to the entire sequence of the squash trypsin inhibitor EETI II (Ecballium elaterium trypsin inhibitor) in which the six cysteines, engaged in three disulphide bridges in native EETI II, have been replaced by six serines, has been synthesised. CD, Fourier-transform infrared spectroscopy (FTIR) and 1H-NMR studies of this peptide revealed that some secondary structures present in native EETI II are still populated in the absence of disulphide bonds. Native-like secondary structures were observed for segments 10-15 (helix), 16-19 and 22-25 (reverse turns) but no native tertiary interaction was detected. However, a non-native local interaction between the aromatic ring of Phe26 and the amide group of Gly28 was observed. It is hypothesised that the 10-15, 16-19 and 22-25 native-like local conformations could play a major role in the early folding of EETI II.

Comparative modeling of the three-dimensional structure of type II antifreeze protein

Type II antifreeze proteins (AFP), which inhibit the growth of seed ice crystals in the blood of certain fishes (sea raven, herring, and smelt), are the largest known fish AFPs and the only class for which detailed structural information is not yet available. However, a sequence homology has been recognized between these proteins and the carbohydrate recognition domain of C-type lectins. The structure of this domain from rat mannose-binding protein (MBP-A) has been solved by X-ray crystallography (Weis WI, Drickamer K, Hendrickson WA, 1992, Nature 360:127-134) and provided the coordinates for constructing the three-dimensional model of the 129-amino acid Type II AFP from sea raven, to which it shows 19% sequence identity. Multiple sequence alignments between Type II AFPs, pancreatic stone protein, MBP-A, and as many as 50 carbohydrate-recognition domain sequences from various lectins were performed to determine reliably aligned sequence regions. Successive molecular dynamics and energy minimization calculations were used to relax bond lengths and angles and to identify flexible regions. The derived structure contains two alpha-helices, two beta-sheets, and a high proportion of amino acids in loops and turns. The model is in good agreement with preliminary NMR spectroscopic analyses. It explains the observed differences in calcium binding between sea raven Type II AFP and MBP-A. Furthermore, the model proposes the formation of five disulfide bridges between Cys 7 and Cys 18, Cys 35 and Cys 125, Cys 69 and Cys 100, Cys 89 and Cys 111, and Cys 101 and Cys 117.(ABSTRACT TRUNCATED AT 250 WORDS)

An 1H NMR determination of the three-dimensional structures of mirror-image forms of a Leu-5 variant of the trypsin inhibitor from Ecballium elaterium (EETI-II)

The 3-dimensional structures of mirror-image forms of a Leu-5 variant of the trypsin inhibitor Ecballium elaterium (EETI-II) have been determined by 1H NMR spectroscopy and simulated annealing calculations incorporating NOE-derived distance constraints. Spectra were assigned using 2-dimensional NMR methods at 400 MHz, and internuclear distances were determined from NOESY experiments. Three-bond spin-spin couplings between C alpha H and amide protons, amide exchange rates, and the temperature dependence of amide chemical shifts were also measured. The structure consists largely of loops and turns, with a short region of beta-sheet. The Leu-5 substitution produces a substantial reduction in affinity for trypsin relative to native EETI-II, which contains an Ile at this position. The global structure of the Leu-5 analogue studied here is similar to that reported for native EETI-II (Heitz A, Chiche L, Le-Nguyen D, Castro B, 1989, Biochemistry 28:2392-2398) and to X-ray and NMR structures of the related proteinase inhibitor CMTI-I (Bode W et al., 1989, FEBS Lett 242:285-292; Holak TA et al., 1989a, J Mol Biol 210:649-654; Holak TA, Gondol D, Otlewski J, Wilusz T, 1989b, J Mol Biol 210:635-648; Holak TA, Habazettl J, Oschkinat H, Otlewski J, 1991, J Am Chem Soc 113:3196-3198). The region near the scissile bond is the most disordered part of the structure, based on geometric superimposition of 40 calculated structures. This disorder most likely reflects additional motion being present in this region relative to the rest of the protein. This motional disorder is increased in the Leu-5 analogue relative to the native form and may be responsible for its reduced trypsin binding. A second form of the protein synthesized with all (D) amino acids was also studied by NMR and found to have a spectrum identical with that of the (L) form. This is consistent with the (D) form being a mirror image of the (L) form and not distinguishable by NMR in an achiral solvent (i.e., H2O). The (D) form has no activity against trypsin, as would be expected for a mirror-image form.

Solution structure of a pyrimidine.purine.pyrimidine DNA triplex containing T.AT, C+.GC and G.TA triples

BACKGROUND

Under certain conditions, homopyrimidine oligonucleotides can bind to complementary homopurine sequences in homopurine-homopyrimidine segments of duplex DNA to form triple helical structures. Besides having biological implications in vivo, this property has been exploited in molecular biology applications. This approach is limited by a lack of knowledge about the recognition by the third strand of pyrimidine residues in Watson-Crick base pairs.

RESULTS

We have therefore determined the solution structure of a pyrimidine.purine.pyrimidine (Y.RY) DNA triple helix containing a guanine residue in the third strand which was postulated to specifically recognize a thymine residue in a Watson-Crick TA base pair. The structure was solved by combining NMR-derived restraints with molecular dynamics simulations conducted in the presence of explicit solvent and counter ions. The guanine of the G-TA triple is tilted out of the plane of its target TA base pair towards the 3'-direction, to avoid a steric clash with the thymine methyl group. This allows the guanine amino protons to participate in hydrogen bonds with separate carbonyls, forming one strong bond within the G-TA triple and a weak bond to an adjacent T.AT triple. Dramatic variations in helical twist around the guanine residue lead to a novel stacking interaction. At the global level, the Y.RY DNA triplex shares several structural features with the recently solved solution structure of the R.RY DNA triplex.

CONCLUSIONS

The formation of a G.TA triple within an otherwise pyrimidine.purine.pyrimidine DNA triplex causes conformational realignments in and around the G.TA triple. These highlight new aspects of molecular recognition that could be useful in triplex-based approaches to inhibition of gene expression and site-specific cleavage of genomic DNA.

Conformational analysis of protein structures derived from NMR data

A study is presented of the conformational characteristics of NMR-derived protein structures in the Protein Data Bank compared to X-ray structures. Both ensemble and energy-minimized average structures are analyzed. We have addressed the problem using the methods developed for crystal structures by examining the distribution of phi, psi, and chi angles as indicators of global conformational irregularity. All these features in NMR structures occur to varying degrees in multiple conformational states. Some measures of local geometry are very tightly constrained by the methods used to generate the structure, e.g., proline phi angles, alpha-helix phi,psi angles, omega angles, and C alpha chirality. The more lightly restrained torsion angles do show increased clustering as the number of overall experimental observations increases. phi, psi, and chi 1 angle conformational heterogeneity is strongly correlated with accessibility but shows additional differences which reflect the differing number of observations possible in NMR for the various side chains (e.g., many for Trp, few for Ser). In general, we find that the core is defined to a notional resolution of 2.0 to 2.3 A. Of real interest is the behavior of surface residues and in particular the side chains where multiple rotameric states in different structures can vary from 10% to 88%. Later generation structures show a much tighter definition which correlates with increasing use of J-coupling information, stereospecific assignments, and heteronuclear techniques. A suite of programs is being developed to address the special needs of NMR-derived structures which will take into account the existence of increased mobility in solution.

Solution structure of the calcium channel antagonist omega-conotoxin GVIA

The three-dimensional solution structure is reported for omega-conotoxin GVIA, which is a potent inhibitor of presynaptic calcium channels in vertebrate neuromuscular junctions. Structures were generated by a hybrid distance geometry and restrained molecular dynamics approach using interproton distance, torsion angle, and hydrogen-bonding constraints derived from 1H NMR data. Conformations of GVIA with low constraint violations converged to a common peptide fold. The secondary structure in the peptide is an antiparallel triple-stranded beta-sheet containing a beta-hairpin and three tight turns. The NMR data are consistent with the region of the peptide from residues S9 to C16 being more dynamic than the rest of the peptide. The peptide has an amphiphilic structure with a positively charged hydrophilic side and an opposite side that contains a small hydrophobic region. Residues that are thought to be important in binding and function are located on the hydrophilic face of the peptide.

Computational challenges for macromolecular structure determination by X-ray crystallography and solution NMR-spectroscopy

Macromolecular structure determination by X-ray crystallography and solution NMR spectroscopy has experienced unprecedented growth during the past decade.

Characterization and 2D NMR study of the stable [9-21, 15-27] 2 disulfide intermediate in the folding of the 3 disulfide trypsin inhibitor EETI II

The three disulfide Ecballium elaterium trypsin inhibitor II (EETI II) reduction with dithiothreitol (DTT) and reoxidation of the fully reduced derivative have been examined. A common stable intermediate has been observed for both processes. Isolation and sequencing of carboxymethylated material showed that the intermediate lacks the [2-19] bridge. The NMR study showed a very strong structural conservation as compared to the native EETI II, suggesting that the bridges are the [9-21] and [15-27] native ones. The differences occurred in sections 2-7 (containing the free cysteine 2 and the Arg 4-Ile 5 active site) and 19-21 (containing the second free cysteine). Distance geometry calculations and restrained molecular dynamics refinements were also in favor of a [9-21, 15-27] arrangement and resulted in a well-conserved (7-28) segment.

Molecular dynamics studies of a DNA-binding protein: 2. An evaluation of implicit and explicit solvent models for the molecular dynamics simulation of the Escherichia coli trp repressor

Although aqueous simulations with periodic boundary conditions more accurately describe protein dynamics than in vacuo simulations, these are computationally intensive for most proteins. Trp repressor dynamic simulations with a small water shell surrounding the starting model yield protein trajectories that are markedly improved over gas phase, yet computationally efficient. Explicit water in molecular dynamics simulations maintains surface exposure of protein hydrophilic atoms and burial of hydrophobic atoms by opposing the otherwise asymmetric protein-protein forces. This properly orients protein surface side chains, reduces protein fluctuations, and lowers the overall root mean square deviation from the crystal structure. For simulations with crystallographic waters only, a linear or sigmoidal distance-dependent dielectric yields a much better trajectory than does a constant dielectric model. As more water is added to the starting model, the differences between using distance-dependent and constant dielectric models becomes smaller, although the linear distance-dependent dielectric yields an average structure closer to the crystal structure than does a constant dielectric model. Multiplicative constants greater than one, for the linear distance-dependent dielectric simulations, produced trajectories that are progressively worse in describing trp repressor dynamics. Simulations of bovine pancreatic trypsin were used to ensure that the trp repressor results were not protein dependent and to explore the effect of the nonbonded cutoff on the distance-dependent and constant dielectric simulation models. The nonbonded cutoff markedly affected the constant but not distance-dependent dielectric bovine pancreatic trypsin inhibitor simulations. As with trp repressor, the distance-dependent dielectric model with a shell of water surrounding the protein produced a trajectory in better agreement with the crystal structure than a constant dielectric model, and the physical properties of the trajectory average structure, both with and without a nonbonded cutoff, were comparable.

Theoretical determination of conformational paths in citrate synthase

Two methods are developed for the theoretical determination of a conformational path between two well-documented forms, a closed form and the open form [Remington et al. (1982) J. Mol. Biol. 158, 111-152] of pig heart citrate synthase, a dimeric enzyme of 2 x 437 residues. The first method uses the minimization of the sum of the potential energies at a set of equidistant points, according to Elber and Karplus [(1987) Chem. Phys. Lett. 139, 375-380]. The initialization of the algorithm is modified to account for large-angle rotations of many groups by performing the interpolations in the space of internal polar coordinates of a set of generalized Jacobi vectors earlier introduced by Durup [(1991) J. Phys. Chem. 95, 1817-1829] and by carefully testing all choices of directions of rotation for determining the initialized midpoint between the known forms. The path includes intermediate points, created by successive splittings of each interval into two equal parts, with a partial energy minimization performed after each splitting. The minimization encounters the well-known local-minima problem, which here is handled by low-temperature molecular dynamics annealing. It is shown that the best ratio of potential energy decrease to rms deviation is achieved by running the dynamics at 50 K, as compared to 100 K and above. The main character of the path obtained is the occurrence of strong to-and-fro variations of some dihedral angles at specific stages along the path. The second method, which we name directed dynamics, uses only low-temperature molecular dynamics simulations by starting trajectories from each of the two known forms with initial velocities directed toward the other one. The procedure is iterated by restarting trajectory pairs after the points of closest approach of the preceding pair. The two half-paths thus built eventually meet after 70 iterations. This method provides a second path with strong similarities, as well as some differences, with respect to the path obtained by the first method.

Inclusion of solvation free energy with molecular mechanics energy: alanyl dipeptide as a test case

A combined force field of molecular mechanics and solvation free energy is tested by carrying out energy minimization and molecular dynamics on several conformations of the alanyl dipeptide. Our results are qualitatively consistent with previous experimental and computational studies, in that the addition of solvation energy stabilizes the C5 conformation of the alanyl dipeptide relative to the C7.

Natural protein proteinase inhibitors and their interaction with proteinases

The substrate-like 'canonical' inhibition by the 'small' serine proteinase inhibitors and the product-like inhibition by the carboxypeptidase inhibitor have provided the only atomic models of protein inhibitor--proteinase interactions for about 15 years. The recently published structures of cystatin/stefin--papain complexes and of hirudin--thrombin complexes reveal novel non-substrate-like interactions. In addition, the structure of pro-carboxypeptidase shows a model of inactivation which bears resemblance to proteinase/protein inhibitor systems. Considerable progress in understanding the transition between native and cleaved states of the serpins has also been made by several recent structural studies.

Conformation of endothelin in aqueous ethylene glycol determined by 1H-NMR and molecular dynamics simulations

The solution conformation of a 21-residue vasoconstrictor peptide endothelin-1 (ET-1) in water-ethylene glycol has been determined by two-dimensional 1H-NMR spectroscopy and constrained molecular dynamics simulations. The N-terminus (residues 1-4) appears to undergo conformational averaging and no single structure consistent with the NMR constraints could be found for this region. Residues 5-8 form a turn, and residues 9-16 exist in a helical conformation. A flexible 'hinge' between residues 8-9 allows various orientations of the turn relative to the helix. Another 'hinge' at residue 17 connects the extended C-terminus to the bicyclic core region (residues 1-15). Residues important for binding and biological activity form a contiguous surface on one side of the helix, with the two disulfides extending from the other side of the helix.

Protein model structure evaluation using the solvation free energy of folding

A systematic study of solvation free energy of folding for proteins with known crystallographic structures is presented. There is a linear relationship between the solvation free energy of folding and the protein size. This relationship, which can be rationalized by a simple model of chain folding, allows prediction of the solvation free energy of folding for proteins for which no high resolution structures are available. All misfolded structures analyzed show solvation free energies higher than predicted; however, some of the misfolded structures have values close enough to the predicted values so that one must be very careful when using such a criterion to check the correctness of a protein model.