Protein folding determinants: structural features determining alternative disulfide pairing in alpha- and chi/lambda-conotoxins

Alpha-conotoxins isolated from Conus venoms contain 11-19 residues and preferentially fold into the globular conformation that possesses a specific disulfide pairing pattern (C1-3, C2-4). We and others isolated a new family of chi-conotoxins (also called lambda conotoxins) with the conserved cysteine framework of alpha-conotoxins but with alternative disulfide pairing (C1-4, C2-3) resulting in the ribbon conformation. In both families, disulfide pairing and hence folding are important for their biological potency. By comparing the structural differences, we identified potential structural determinants responsible for the folding tendencies of these conotoxins. We examined the role of conserved proline in the first intercysteine loop and the conserved C-terminal amide on folding patterns of synthetic analogues of ImI conotoxin by comparing the isoforms with the regiospecifically synthesized conformers. Deamidation at the C-terminus and substitution of proline in the first intercysteine loop switch the folding pattern from the globular form of alpha-conotoxins to the ribbon form of chi/lambda-conotoxins. The findings are corroborated by reciprocal folding of CMrVIA chi/lambda-conotoxins. Substitution of Lys-6 from the first intercysteine loop of CMrVIA conotoxin with proline, as well as the inclusion of an amidated C-terminal shifted the folding preference of CMrVIA conotoxin from its native ribbon conformation toward the globular conformation. Binding assays of ImI conotoxin analogues with Aplysia and Bulinus acetylcholine binding protein indicate that both these substitutions and their consequent conformational change substantially impact the binding affinity of ImI conotoxin. These results strongly indicate that the first intercysteine loop proline and C-terminal amidation act as conformational switches in alpha- and chi/lambda-conotoxins.     

lambda-conotoxins, a new family of conotoxins with unique disulfide pattern and protein folding. Isolation and characterization from the venom of Conus marmoreus

Conotoxins are multiple disulfide-bonded peptides isolated from marine cone snail venom. These toxins have been classified into several families based on their disulfide pattern and biological properties. Here, we report a new family of Conus peptides, which have a novel cysteine motif. Three peptides of this family (CMrVIA, CMrVIB, and CMrX) have been purified from Conus marmoreus venom, and their structures have been determined. Their amino acid sequences are VCCGYK-LCHOC (CMrVIA), NGVCCGYKLCHOC (CMrVIB), and GICCGVSFCYOC (CMrX), where O represents 4-trans-hydroxyproline. Two of these peptides (CMrVIA and CMrX) have been chemically synthesized. Using a selective protection and deprotection strategy during disulfide bond formation, peptides with both feasible cysteine-pairing combinations were generated. The disulfide pattern (C(1)-C(4), C(2)-C(3)) in native toxins was identified by their co-elution with the synthetic disulfide-isomeric peptides on reverse-phase high pressure liquid chromatography. Although cysteine residues were found in comparable positions with those of alpha-conotoxins, these toxins exhibited a distinctly different disulfide bonding pattern; we have named this new family "lambda -conotoxins." CMrVIA and CMrX induced different biological effects when injected intra-cerebroventricularly in mice; CMrVIA induces seizures, whereas CMrX induces flaccid paralysis. The synthetic peptide with lambda-conotoxin folding is about 1150-fold more potent in inducing seizures than the mispaired isomer with alpha-conotoxin folding. Thus it appears that the unique disulfide pattern, and hence the "ribbon" conformation, in lambda-conotoxins is important for their biological activity.     

A new level of conotoxin diversity,a non-native disulfide bond connectivity in alpha-conotoxin AuIB reduces structural definition but increases biological activity

alpha-Conotoxin AuIB and a disulfide bond variant of AuIB have been synthesized to determine the role of disulfide bond connectivity on structure and activity. Both of these peptides contain the 15 amino acid sequence GCCSYPPCFATNPDC, with the globular (native) isomer having the disulfide connectivity Cys(2-8 and 3-15) and the ribbon isomer having the disulfide connectivity Cys(2-15 and 3-8). The solution structures of the peptides were determined by NMR spectroscopy, and their ability to block the nicotinic acetylcholine receptors on dissociated neurons of the rat parasympathetic ganglia was examined. The ribbon disulfide isomer, although having a less well defined structure, is surprisingly found to have approximately 10 times greater potency than the native peptide. To our knowledge this is the first demonstration of a non-native disulfide bond isomer of a conotoxin exhibiting greater biological activity than the native isomer.     

The ribbon of hydrogen bonds and the pseudomolecule in the three-dimensional structure of globular proteins. III. Bovine pancreatic ribonuclease A and bovine seminal ribonuclease

The model of the three-dimensional structure of globular proteins, which is based on a ribbon of hydrogen bonds along the whole of the backbone, is now applied to the comparison between monomeric bovine pancreatic ribonuclease A and dimeric bovine seminal ribonuclease. Some waters are involved in the hydrogen bonding of the ribbon, and the protein molecule plus these waters forms a pseudomolecule. The conformations of the three backbones are essentially identical and the three ribbons of hydrogen bonds are conserved with greater than 90% accuracy. We suggest that the conservation of the backbone conformations of the two molecules is a consequence of the conservation of the ribbons of hydrogen bonds. There are 16 simple mutations between the two molecules, of which 15 involve only side-chain groups with no more than one hydrogen bond to the backbone. Such mutations are not sufficient to change the ribbon of hydrogen bonds and hence there is no change in the backbone conformation. Generalizing this result, we suggest that the conservation of the ribbon is the reason why single point mutations rarely change the conformation of the backbone of the globular proteins.     

The pseudomolecule method and the structure of globular proteins. II. The example of ribonuclease F1 and T1

The pseudomolecule approach to the structure of globular proteins in which a small number of water molecules are incorporated into the "molecule" is tested again by comparing the ribbon of hydrogen bonds in two proteins, ribonuclease F1 and T1. These two molecules are 59% homologous and have the same backbone conformation both globally and locally. The two ribbons of hydrogen bonds that cover the whole of the backbone are conserved with an accuracy of some 95% providing that allowance is made for the intrusion into one of the pair of such extra factors as the presence of adducts or metal ions, the insertions and the absence of a few water molecules from one of the x-ray data sets. Without these corrections, the conservation of the ribbon is some 85%. There are 35 conserved hydrogen-bonding residues, nearly all of which show many unions to the backbone or interactions with the active site. There are 36 point mutations that involve one or two hydrogen-bonding side chains and nearly all of these have either none or one hydrogen bond to the backbone. These are minor contributors to the ribbon of hydrogen bonds. Of the 71 residues involved in these two categories, all but six fit into the pseudomolecular picture of the structure of globular proteins. The remaining 30 residues almost all contain conserved hydrocarbon side chains that may have a second order effect on the structure through their space filling effects.     

The Ribbon of Hydrogen Bonds in the Three-Dimensional Structure of Globular Proteins

We report quantum mechanical computations and experimental evidence which suggest that the backbone conformation of globular proteins depends generally on the conservation of that part of the hydrogen bond network or ribbon which is joined, in general, directly to the backbone and is largely independent of the remainder of this whole network of hydrogen bonds. The familiar hydrogen bonds of the helix and the sheet form about one-half of this ribbon of hydrogen bonds. Both water molecules and hydrogen bonding side chain groups are involved in the formation of the ribbon.This view of the three-dimensional structure of globular proteins in terms of the `molecule' allows us to deal with the non-secondary structure as well as with the familiar secondary structure. It also suggests that the ribbon contains approximately the same number of hydrogen bonds within all three structures – the helix, the sheet and the coil – and that this is the reason for the ease of interconversion of these three structures.The quantum mechanical computations on hydrogen bonding suggest that delocalised water molecules which have substantial mobility are an essential part of the ribbon. This situation arises because the hydrogen bonding groups of the protein molecule are not free to move to optimise the hydrogen bonding geometries as are the oxygen atoms in the waters and ices. Such delocalised water molecules either have high B values or are invisible in the X-ray data and yet are able to form a structure which is as strong as a normal hydrogen bond.The experimental data on the point mutations of the THRI57 residue of the T4 phage lysome provides an initial test of this model. Both the local backbone conformation and the ribbon of hydrogen bonds are conserved throughout all the mutations of residue 157,providing that the delocalised water molecules are accepted as a genuine part of the structure. These mutations include the introduction of hydrocarbon side chains at position 157 when water molecules or other side chain groups take over the formation of the hydrogen bonds.We suggest that, provided steric effects are not important, many point mutations succeed because they leave the ribbon of hydrogen bonds (and so the backbone conformation) largely unchanged.     

Stable monomeric intermediate with exposed Cys-119 is formed during heat denaturation of beta-lactoglobulin

The role of the free sulfhydryl group of beta-lactoglobulin in the formation of a stable non-native monomer during heat-treatment of beta-lactoglobulin solutions was investigated. Two concomitant events occurred at the earlier stage of heating: unfolding of native globular monomer and intramolecular sulfhydryl/disulfide exchange reaction. Thus, two denatured monomeric species were formed: a non-native monomer with exposed Cys-121 (Mcys121) which became reversible after cooling, and a stable non-native monomer with exposed Cys-119 (Mcys119) which exhibited both a larger hydrodynamic conformation than native monomer and low solubility at pH 4.7. The results also show that the formation of these monomeric species throughout heat-induced denaturation of native beta-lg monomers is faster than their subsequent aggregation. A mechanism describing the behavior of beta-lg denaturation/aggregation during heat-treatment under selected conditions (5.8 mg/ml, low ionic strength, pH 6.6, 85 degrees C) is presented.     

Thermodynamics of interactions of urea and guanidinium salts with protein surface: relationship between solute effects on protein processes and changes in water-accessible surface area.

To interpret effects of urea and guanidinium (GuH(+)) salts on processes that involve large changes in protein water-accessible surface area (ASA), and to predict these effects from structural information, a thermodynamic characterization of the interactions of these solutes with different types of protein surface is required. In the present work we quantify the interactions of urea, GuHCl, GuHSCN, and, for comparison, KCl with native bovine serum albumin (BSA) surface, using vapor pressure osmometry (VPO) to obtain preferential interaction coefficients (Gamma(mu3)) as functions of nondenaturing concentrations of these solutes (0-1 molal). From analysis of Gamma(mu3) using the local-bulk domain model, we obtain concentration-independent partition coefficients K(nat)(P) that characterize the accumulation of these solutes near native protein (BSA) surface: K(nat)(P,urea)= 1.10 +/- 0.04, K(nat)(P,SCN(-)) = 2.4 +/- 0.2, K(nat)(P,GuH(+)) = 1.60 +/- 0.08, relative to K(nat)(P,K(+)) identical with 1 and K(nat)(P,Cl(-)) = 1.0 +/- 0.08. The relative magnitudes of K(nat)(P) are consistent with the relative effectiveness of these solutes as perturbants of protein processes. From a comparison of partition coefficients for these solutes and native surface (K(nat)(P)) with those determined by us previously for unfolded protein and alanine-based peptide surface K(unf)(P), we dissect K(P) into contributions from polar peptide backbone and other types of protein surface. For globular protein-urea interactions, we find K(nat)(P,urea) = K(unf)(P,urea). We propose that this equality arises because polar peptide backbone is the same fraction (0.13) of total ASA for both classes of surface. The analysis presented here quantifies and provides a physical basis for understanding Hofmeister effects of salt ions and the effects of uncharged solutes on protein processes in terms of K(P) and the change in protein ASA.     

Excluded volume in protein side-chain packing

The excluded volume occupied by protein side-chains and the requirement of high packing density in the protein interior should severely limit the number of side-chain conformations compatible with a given native backbone. To examine the relationship between side-chain geometry and side-chain packing, we use an all-atom Monte Carlo simulation to sample the large space of side-chain conformations. We study three models of excluded volume and use umbrella sampling to effectively explore the entire space. We find that while excluded volume constraints reduce the size of conformational space by many orders of magnitude, the number of allowed conformations is still large. An average repacked conformation has 20 % of its chi angles in a non-native state, a marked reduction from the expected 67 % in the absence of excluded volume. Interestingly, well-packed conformations with up to 50 % non-native chi angles exist. The repacked conformations have native packing density as measured by a standard Voronoi procedure. Entropy is distributed non-uniformly over positions, and we partially explain the observed distribution using rotamer probabilities derived from the Protein Data Bank database. In several cases, native rotamers that occur infrequently in the database are seen with high probability in our simulation, indicating that sequence-specific excluded volume interactions can stabilize rotamers that are rare for a given backbone. In spite of our finding that 65 % of the native rotamers and 85 % of chi(1) angles can be predicted correctly on the basis of excluded volume only, 95 % of positions can accommodate more than one rotamer in simulation. We estimate that, in order to quench the side-chain entropy observed in the presence of excluded volume interactions, other interactions (hydrophobic, polar, electrostatic) must provide an additional stabilization of at least 0.6 kT per residue in order to single out the native state.     

Lysine 5 and phenylalanine 9 of the factor IX omega-loop interact with phosphatidylserine in a membrane-mimetic environment

The binding of factor IX to cell membranes requires a structured N-terminal omega-loop conformation that exposes hydrophobic residues for a highly regulated interaction with a phospholipid. We hypothesized that a peptide comprised of amino acids Gly4-Gln11 of factor IX (fIX(G4)(-)(Q11)) and constrained by an engineered disulfide bond would assume the native factor IX omega-loop conformation in the absence of Ca(2+). The small size and freedom from aggregation-inducing calcium interactions would make fIX(G4)(-)(Q11) suitable for structural studies for eliciting details about phospholipid interactions. fIX(G4)(-)(Q11) competes with factor IXa for binding sites on phosphatidylserine-containing membranes with a K(i) of 11 microM and inhibits the activation of factor X by the factor VIIIa-IXa complex with a K(i) of 285 microM. The NMR structure of fIX(G4)(-)(Q11) reveals an omega-loop backbone fold and side chain orientation similar to those found in the calcium-bound factor IX Gla domain, FIX(1-47)-Ca(2+). Dicaproylphosphatidylserine (C(6)PS) induces HN, Halpha backbone, and Hbeta chemical shift perturbations at residues Lys5, Leu6, Phe9, and Val10 of fIX(G4)(-)(Q11), while selectively protecting the NHzeta side chain resonance of Lys5 from solvent exchange. NOEs between the aromatic ring protons of Phe9 and specific acyl chain protons of C(6)PS indicate that these phosphatidylserine protons reside 3-6 A from Phe9. Stabilization of the phosphoserine headgroup and glycerol backbone of C(6)PS identifies that phosphatidylserine is in a protected environment that is spatially juxtaposed with fIX(G4)(-)(Q11). Together, these data demonstrate that Lys5, Leu6, Phe9, and Val10 preferentially interact with C(6)PS and allow us to correlate known hemophilia B mutations of factor IX at Lys5 or Phe9 with impaired phosphatidylserine interaction.     

Molecular dynamics simulations of Trp side-chain conformational flexibility in the gramicidin A channel

Gramicidin A (gA) is prototypical peptide antibiotic and a model ion channel former. Configured in the solid-state NMR beta(6.5)-helix channel conformation, gA was subjected to 1-ns molecular dynamics (MD) gas phase simulations using the all-atom charmm22 force field to ascertain the conformational stability of the Trp side chains as governed by backbone and neighboring side-chain contacts. Three microcanonical trajectories were computed using different initial atomic velocities for each of twenty different initial structures. For each set, one of the four Trp side chains in each monomer was initially positioned in one of the five non-native conformations (A. E. Dorigo et al., Biophysical Journal, 1999, Vol. 76, 1897-1908), the other Trps being positioned in the native state, o1. In three additional control simulations, all Trps were initiated in the native conformation. After equilibration, constraints were removed and subsequent conformational changes of the initially constrained Trp were measured. The chi(1) was more flexible than chi(2.1). The energetically optimal orientation, o1 (Dorigo et al., 1999), was the most stable in all four Trp positions (9, 11, 13, 15) and remained unchanged for the entire 1 ns simulation in 19 of 24 trials. Changes in chi(1) from each of the 5 suboptimal states occur readily. Two of the non-native conformations reverted readily to o1, whereas the other three converted to an intermediate state, i2. There were frequent interconversions between i2 and o1. We speculate that experimentally observed Trp stability is caused by interactions with the lipid-water interface, and that stabilization of one of the suboptimal conformations in gA, such as i2, by lipid headgroups could produce a secondary, metastable conformational state. This could explain recent experimental studies of differences in the channel conductance dispersity between gA and a Trp-to-Phe gA analog, gramicidin M (gM, J. C. Markham et al., Biochimica et Biophysica Acta, 2001, Vol. 1513, 185-192).     

Influences of disulfide connectivity on structure and antimicrobial activity of tachyplesin I

Tachyplesin I is a potent antimicrobial peptide with broad spectrum of antimicrobial activity. It has 2 disulfide bonds and can form 3 disulfide bond isomers. In this study, the structure and antimicrobial activity of 3 tachyplesin I isomers (tachyplesin I, 3C12C, 3C7C) were investigated using molecular dynamic simulations, circular dichroism structural study, as well as antimicrobial activity and hemolysis assay. Our results suggest that in comparison to the native peptide, the 2 isomers (3C12C, 3C7C) have substantial structural and activity variations. The native peptide is in the ribbon conformation, while 3C12C and 3C7C possess remarkably different secondary structures, which are referred as “globular” and “beads” isomers, respectively. The substantially decreased hemolysis effects for these 2 isomers is accompanied by significantly decreased anti‐gram‐positive bacterial activity.     

alpha4/7-conotoxin Lp1.1 is a novel antagonist of neuronal nicotinic acetylcholine receptors

Cone snails comprise approximately 700 species of venomous molluscs which have evolved the ability to generate multiple toxins with varied and exquisite selectivity. alpha-Conotoxin is a powerful tool for defining the composition and function of nicotinic acetylcholine receptors which play a crucial role in excitatory neurotransmission and are important targets for drugs and insecticides. An alpha4/7 conotoxin, Lp1.1, originally identified by cDNA and genomic DNA cloning from Conus leopardus, was found devoid of the highly conserved Pro residue in the first intercysteine loop. To further study this toxin, alpha-Lp1.1 was chemically synthesized and refolded into its globular disulfide isomer. The synthetic Lp1.1 induced seizure and paralysis on freshwater goldfish and selectively reversibly inhibited ACh-evoked currents in Xenopus oocytes expressing rat alpha3beta2 and alpha6alpha3beta2 nAChRs. Comparing the distinct primary structure with other functionally related alpha-conotoxins could indicate structural features in Lp1.1 that may be associated with its unique receptor recognition profile.     

Biostructural chemistry of magnesium ion: characterization of the weak binding sites on tRNA(Phe)(yeast). Implications for conformational change and activity

The thermodynamics and kinetics of magnesium binding to tRNA(Phe)(yeast) have been studied directly by 25Mg NMR. In 0.17 M Na+(aq), tRNA(Phe) exists in its native conformation and the number of strong binding sites (Ka greater than or equal to 10(4)) was estimated to be 3-4 by titration experiments, in agreement with X-ray structural data for crystalline tRNA(Phe) (Jack et al., 1977). The set of weakly bound ions were in slow exchange and 25Mg NMR resonances were in the near-extreme-narrowing limit. The line shapes of the exchange-broadened magnesium resonance were indistinguishable from Lorentzian form. The number of weak magnesium binding sites was determined to be 50 +/- 8 in the native conformation and a total line-shape analysis of the exchange-broadened 25 Mg2+ NMR resonance gave an association constant Ka of (2.2 +/- 0.2) X 10(2) M-1, a quadrupolar coupling constant (chi B) of 0.84 MHz, an activation free energy (delta G*) of 12.8 +/- 0.2 kcal mol-1, and an off-rate (koff) of (2.5 +/- 0.4) X 10(3) s-1. In the absence of background Na+(aq), up to 12 +/- 2 magnesium ions bind cooperatively, and 73 +/- 10 additional weak binding sites were determined. The binding parameters in the nonnative conformation were Ka = (2.5 +/- 0.2) X 10(2) M-1, chi B = 0.64 MHz, delta G* = 13.1 +/- 0.2 kcal mol-1, and koff = (1.6 +/- 0.4) X 10(3) s-1. In comparison to Mg2+ binding to proteins (chi B typically ca. 1.1-1.6 MHz) the lower chi B values suggest a higher degree of symmetry for the ligand environment of Mg2+ bound to tRNA. A small number of specific weakly bound Mg2+ appear to be important for the change from a nonnative to a native conformation. Implications for interactions with the ribosome are discussed.     

Pc16a, the first characterized peptide from Conus pictus venom, shows a novel disulfide connectivity

A novel conotoxin, pc16a, was isolated from the venom of Conus pictus. This is the first peptide characterized from this South-African cone snail and it has only 11 amino acid residues, SCSCKRNFLCC*, with the rare cysteine framework XVI and a monoisotopic mass of 1257.6Da. Two peptides were synthesized with two possible conformations: globular (pc16a_1) and ribbon (pc16a_2). pc16a_1 co-eluted with the native peptide, which indicates a disulfide connectivity I-III, II-IV. The structure of pc16a_1 was determined by NMR. Both synthetic peptides were used to elucidate the biological activity. Bioassays were performed on crickets, ghost shrimps, larvae of the mealworm beetle and mice, but no effect was seen. Using two-electrode voltage clamp, a range of voltage-gated ion channels (Na(v) and K(v)) and nicotinic acetylcholine receptors were screened, but again no activity was found. Hence, the specific target of pc16a still remains to be discovered.     

Structural and dynamic characteristics of a partially folded state of ubiquitin revealed by hydrogen exchange mass spectrometry

Structural and dynamic properties of a partially folded conformation (A-state) of ubiquitin are studied using amide hydrogen exchange in solution (HDX) and mass spectrometric detection. A clear distinction between the native state of the protein and the A-state can be made when HDX is carried out in a semicorrelated regime. Convoluted exchange patterns are interpreted with the aid of HDX simulations in a three-state system (highly structured, partially unstructured, and fully unstructured states). The data clearly indicate a highly dynamic character of the non-native state. Furthermore, combination of HDX and protein ion fragmentation in the gas phase [by means of collision-induced dissociation (CAD)] is used to evaluate the conformational stability of various protein segments specifically in the molten globular state. Chain flexibility appears to be distributed very unevenly in this non-native conformation. The highest degree of structural disorder is displayed by the C-terminal segment (Gly(53)-Gly(76)), which was previously suggested to form a transient alpha-helix. The least dynamic segment of ubiquitin in the A-state is Thr(9)-Glu(18) (which was previously suggested to form a stable nativelike beta-strand), with the adjacent segments exhibiting somewhat diminished conformational stability. The study also demonstrates the power of mass spectrometry as a tool in providing conformer-specific information about the structure and dynamics of both native and non-native protein states coexisting in solution under equilibrium.     

mr1e, a conotoxin from Conus marmoreus with a novel disulfide pattern

Conotoxins are well known for their highly variable structures and functions. Here we report the identification of a novel conotoxin named mr1e from Conus marmoreus . mr1e is composed of 11 amino acid residues cross-linked by two disulfide bonds (CCHSSWCKHLC). The spacing of intercysteine loops in mr1e is exactly the same as that in α4/3 conotoxins. However, the native mr1e peptide co-eluted on reverse-phase HPLC with the regioselectively synthesized ribbon disulfide linkage isomer (C¹-C⁴, C²-C³) but not the globular linkage isomer (C¹-C³, C²-C⁴). Although this peptide has the same disulfide connectivity as the χ-conotoxins, their sequences do not share significant homology. Thus, mr1e could be defined as a novel conotoxin family. By intracranial injection into mice, mr1e showed an excitatory effect. The characterization of mr1e certainly enriches our understanding of conotoxins, and also opens an avenue for further structural and functional investigation.     

Synthesis, crystal structure, and molecular conformation of peptide N-Boc-L-Pro-dehydro-Phe-L-Gly-OH

The peptide N-Boc-L-Pro-dehydro-Phe-L-Gly-OH was synthesized by the usual workup procedure and finally coupling the N-Boc-L-Pro-dehydro-Phe to glycine. The peptide crystallizes in monoclinic space group P2(1) with a = 8.951(4) A, b = 5.677(6) A, c = 21.192(11) A, beta = 96.97(4) degrees, V = 1069(1) A3, Z = 2, dm = 1.295(5) Mgm-3, and dc = 1.297(4) Mgm-3. The structure was determined by direct methods using SHELXS86. The structure was refined by the block-diagonal least-squares procedure to an R value of 0.074 for 1002 observed reflections. The C alpha 2-C beta 2 distance of 1.33(2) A is an appropriate double bond length. The angle C alpha 2-C beta 2-C gamma 2 is 133(1) degrees. The peptide backbone torsion angles are theta 1 = -167(1) degrees, omega 0 = 179(1) degrees, phi 1 = -48(1) degrees, psi 1 = 137(1) degrees, omega 1 = 175(1) degrees, phi 2 = 65(2) degrees, psi 2 = 15(2) degrees, omega 2 = -179(1) degrees, and phi 3 = -166(1) degrees. These values show that the Boc group has a trans-trans conformation while the peptide backbone adopts a beta-turn II conformation, which is stabilized by an intramolecular hydrogen bond of length of 3.05(1) A. The structures of dehydro-Phe containing peptides suggest that the dehydro-Phe promotes the beta-turn II conformation. The five-membered pyrrolidine ring of the Pro residue adopts an ideal C gamma-exo conformation with torsion angles chi 1(1) = -24(1) degrees, chi 2(1) = 34(1) degrees, chi 3(1) = -30(1) degrees, chi 4(1) = 15(1) degrees, and theta 0(1) = 6(1) degrees. The side-chain torsion angles in dehydro-Phe are chi 1(2) = -1(2) degrees, chi 2,1(2) = -176(1) degrees, and chi 2,2(2) = 8(2) degrees. The plane of C alpha 2-C beta 2-C gamma 2 is rotated with respect to the plane of the phenyl ring at 7(1) degrees, which indicates that the atoms of the side chain of dehydro-Phe are essentially coplanar. The molecules form a 2(1) screw axis related hydrogen-bonded rows along the b axis.     

A new computational model for protein folding based on atomic solvation.

A new model for calculating the solvation energy of proteins is developed and tested for its ability to identify the native conformation as the global energy minimum among a group of thousands of computationally generated compact non-native conformations for a series of globular proteins. In the model (called the WZS model), solvation preferences for a set of 17 chemically derived molecular fragments of the 20 amino acids are learned by a training algorithm based on maximizing the solvation energy difference between native and non-native conformations for a training set of proteins. The performance of the WZS model confirms the success of this learning approach; the WZS model misrecognizes (as more stable than native) only 7 of 8,200 non-native structures. Possible applications of this model to the prediction of protein structure from sequence are discussed.     

Statistical conformation of human plasma fibronectin

Fibronectin is a multifunctional glycoprotein (molecular mass, M = 530 kg/mol) of the extra cellular matrix (ECM) having a major role in cell adhesion. In physiological conditions, the conformation of this protein still remains debated and controversial. Here, we present a set of results obtained by scattering experiments. In "native" conditions, the radius of gyration (R(g) = 15.3 +/- 0.3 nm) was determined by static light scattering as well as small-angle neutron scattering. The hydrodynamic radius (R(H) = 11.5 +/- 0.1 nm) was deduced from quasi-elastic light scattering measurements. These results imply a low internal concentration compared to that of usual globular proteins. This is also confirmed by the ratio R(H)/R(g) = 0. 75 +/- 0.02 consistent with a Gaussian chain, whereas R(H)/R(g) = 1. 3 for spherical shaped molecules. However, adding a denaturing agent (urea 8 M) increases R(g) by a factor 2. This means that fibronectin "native" chain is not either completely unfolded. The average shape of fibronectin conformation was also probed by small-angle neutron scattering performed for reverse scattering vector q(-)(1) smaller than R(g) (0.2 < q(-)(1) < 15 nm). The measured form factor is in complete agreement with the form factor of a random string of 56 beads of 5 nm diameter. It rules out the possibility of unfolded chain as well as globular structures. These results have structural and biological implications as far as ECM organization is concerned.