Statistical and molecular dynamics studies of buried waters in globular proteins

Buried solvent molecules are common in the core of globular proteins and contribute to structural stability. Folding necessitates the burial of polar backbone atoms in the protein core, whose hydrogen-bonding capacities should be satisfied on average. Whereas the residues in alpha-helices and beta-sheets form systematic main-chain hydrogen bonds, the residues in turns, coils and loops often contain polar atoms that fail to form intramolecular hydrogen bonds. The statistical analysis of 842 high resolution protein structures shows that well-resolved, internal water molecules preferentially reside near residues without alpha-helical and beta-sheet secondary structures. These buried waters most often form primary hydrogen bonds to main-chain atoms not involved in intramolecular hydrogen bonds, providing strong evidence that hydrating main-chain atoms is a key structural role of buried water molecules. Additionally, the average B-factor of protein atoms hydrogen-bonded to waters is smaller than that of protein atoms forming intramolecular hydrogen bonds, and the average B-factor of water molecules involved in primary hydrogen bonds with main-chain atoms is smaller than the average B-factor of water molecules involved in secondary hydrogen bonds to protein atoms that form concurrent intramolecular hydrogen bonds. To study the structural coupling between internal waters and buried polar atoms in detail we simulated the dynamics of wild-type FKBP12, in which a buried water, Wat137, forms one side-chain and multiple main-chain hydrogen bonds. We mutated E60, whose side-chain hydrogen bonds with Wat137, to Q, N, S or A, to modulate the multiplicity and geometry of hydrogen bonds to the water. Mutating E60 to a residue that is unable to form a hydrogen bond with Wat137 results in reorientation of the water molecule and leads to a structural readjustment of residues that are both near and distant to the water. We predict that the E60A mutation will result in a significantly reduced affinity of FKBP12 for its ligand FK506. The propensity of internal waters to hydrogen bond to buried polar atoms suggests that ordered water molecules may constitute fundamental structural components of proteins, particularly in regions where alpha-helical or beta-sheet secondary structure is not present.     

The interrelationships of side-chain and main-chain conformations in proteins

The accurate determination of a large number of protein structures by X-ray crystallography makes it possible to conduct a reliable statistical analysis of the distribution of the main-chain and side-chain conformational angles, how these are dependent on residue type, adjacent residue in the sequence, secondary structure, residue-residue interactions and location at the polypeptide chain termini. The interrelationship between the main-chain (phi, psi) and side-chain (chi 1) torsion angles leads to a classification of amino acid residues that simplify the folding alphabet considerably and can be a guide to the design of new proteins or mutational studies. Analyses of residues occurring with disallowed main-chain conformation or with multiple conformations shed some light on why some residues are less favoured in thermophiles.     

Investigation of the solution conformation of cytochrome c-551 from Pseudomonas stutzeri.

1H NMR spectroscopy and solution structure computations have been used to examine ferrocytochrome c-551 from Pseudomonas stutzeri (ATCC 17588). Resonance assignments are proposed for all main-chain and most side-chain protons. Distance constraints were determined on the basis of nuclear Overhauser enhancements between pairs of protons. Dihedral angle constraints were determined from estimates of scaler coupling constants. Twenty-four structures were calculated by distance geometry and refined by energy minimization and simulated annealing on the basis of 1033 interproton distance and 57 torsion angle constraints. Both the main-chain and side-chain atoms are well defined except for a loop region around residues 34-40, the first two residues at the N-terminus and the last two at the C-terminus, and some side chains located on the molecular surface. The average root mean squared deviation in position for equivalent atoms between the 24 individual structures and the mean structure obtained by averaging their coordinates is 0.54 +/- 0.08 A for the main-chain atoms and 0.97 +/- 0.09 A for all non-hydrogen atoms of residues 3-80 plus the heme group. These structures were compared to the X-ray crystallographic structure of an analogous protein, cytochrome c-551 from Pseudomonas aeruginosa [Matsuura, Takano, & Dickerson (1982) J. Mol. Biol. 156, 389-409). The main-chain folding patterns are very consistent, but there are some differences. The largest difference is in a surface loop segment from residues 34 to 40.     

Explicit-solvent molecular dynamics simulations of the polysaccharide schizophyllan in water

Schizophyllan is a beta(1-->3)-D-glucan polysaccharide with beta(1-->6)-branched lateral glucose residues that presents a very stiff triple-helical structure under most experimental conditions. Despite the remarkable stability of this structure (which persists up to 120 degrees C in aqueous solution), schizophyllan undergoes a major change of state around 7 degrees C in water that has been hypothesized to result from an order-disorder transition in the lateral residues. This hypothesis is only supported by indirect experimental evidence and detailed knowledge (at the atomic level) concerning hydrogen-bonding networks, interactions with the solvent molecules, orientational freedom of the lateral residues, and orientational correlations among them is still lacking. In this study explicit-solvent molecular dynamics simulations of a schizophyllan fragment (complemented by simulations of its tetrasaccharide monomer) are performed at three different temperatures (273 K, 350 K, and 450 K) and with two different types of boundary conditions (finite nonperiodic or infinite periodic fragment) as an attempt to provide detailed structural and dynamical information about the triple-helical conformation in solution and the mechanism of the low-temperature transition. These simulations suggest that three important driving forces for the high stability of the triple helix are i), the limited conformational work involved in its formation; ii), the formation of a dense hydrogen-bonding network at its center; and iii), the formation of interchain hydrogen bonds between main-chain and lateral glucose residues. However, these simulations evidence a moderate and continuous variation of the simulated observables upon increasing the temperature, rather than a sharp transition between the two lowest temperatures (that could be associated with the state transition). Although water-mediated hydrogen-bonded association of neighboring lateral residues is observed, this interaction is not strong enough to promote the formation of an ordered state (correlated motions of the lateral residues), even at the lowest temperature considered.     

Sequence-specific NMR assignments of the trp repressor from Escherichia coli using three-dimensional 15N/1H heteronuclear techniques.

Sequence-specific 15N and 1H assignments for the trp holorepressor from Escherichia coli are reported. The trp repressor consists of two identical 107-residue subunits which are highly helical in the crystal state [Schevitz, R., Otwinowski, Z., Joachimiak, A., Lawson, C. L. & Sigler, P. B. (1985) Nature 317, 782-786]. The high helical content and the relatively large size of the protein (Mr = 25,000) make it difficult to assign even the main-chain resonances by conventional homonuclear two-dimensional NMR methods. However, we have now assigned the main-chain resonances of 94% of the residues by using three-dimensional 15N/1H heteronuclear experiments on a sample of protein uniformly labelled with 15N. The additional resolution obtained by spreading out the signals into three dimensions proved indispensable in making these assignments. In particular, we have been able to resolve signals from residues in the N-terminal region of the A helix for the first time in solution. The observed NOE results confirm that the repressor is highly helical in solution, and contains no extended chain conformations.     

Electronic and vibrational circular dichroism of aromatic amino acids by density functional theory

Structures of model compounds mimicking aromatic amino acid residues in proteins are optimized by density functional theory (DFT), assuming that the main-chain conformation was a random coil. Excitation energies and dipole and rotational strengths for the optimized structures were calculated based on time-dependent DFT (TD-DFT). The electronic circular dichroism (ECD) bands of the models were significantly affected by side-chain conformations. Hydration models of the aromatic residues were also subjected to TD-DFT calculations, and the ECD bands of these models were found to be highly perturbed by the hydration of the main-chain amide groups. In addition to calculating the random-coil conformation, we also performed TD-DFT calculations of the aromatic residue models, assuming that the main-chain conformation was an alpha-helix or beta-strand. As expected, the overall feature of the ECD bands was also perturbed by the main-chain conformations. Moreover, vibrational circular dichroism (VCD) spectra of the hydration models in a random-coil structure were simulated by DFT, which showed that the VCD spectra are more sensitive to the side-chain conformations than the ECD spectra. The present results show that analyses combining ECD and VCD spectroscopy and using DFT calculations can elucidate the main- and side-chain conformations of aromatic residues in proteins.     

Thiol proteases. Comparative studies based on the high-resolution structures of papain and actinidin, and on amino acid sequence information for cathepsins B and H, and stem bromelain

An accurate three-dimensional structure is known for papain (1.65 A resolution) and actinidin (1.7 A). A detailed comparison of these two structures was performed to determine the effect of amino acid changes on the conformation. It appeared that, despite only 48% identity in their amino acid sequence, different crystallization conditions and different X-ray data collection techniques, their structures are surprisingly similar with a root-mean-square difference of 0.40 A between 76% of the main-chain atoms (differences less than 3 sigma). Insertions and deletions cause larger differences but they alter the conformation over a very limited range of two to three residues only. Conformations of identical side-chains are generally retained to the same extent as the main-chain conformation. If they do change, this is due to a modified local environment. Several examples are described. Spatial positions of hydrogen bonds are conserved to a greater extent than are the specific groups involved. The greatest structural similarity is found for the active site residues of papain and actinidin, for the internal water molecules and for the main-chain conformation of residues in alpha-helices and anti-parallel beta-sheet structure. This was reflected also in the similarity of the temperature factors. It suggests that the secondary structural elements form the skeleton of the molecule and that their interaction is the main factor in directing the fold of the polypeptide chain. Therefore, substitution of residues in the skeleton will, in general, have the most drastic effect on the conformation of the protein molecule. In papain and actinidin, some main-chain-side-chain hydrogen bonds are also strongly conserved and these may determine the folding of non-repetitive parts of the structure. Furthermore, we included primary structure information for three homologous thiol proteases: stem bromelain, and the cathepsins B and H. By combining the three-dimensional structural information for papain and actinidin with sequence homologies and identities, we conclude that the overall folding pattern of the polypeptide chain is grossly the same in all five proteases, and that they utilize the same catalytic mechanism.     

Properties of polyproline II, a secondary structure element implicated in protein-protein interactions

The polyproline II (PPII) conformation of protein backbone is an important secondary structure type. It is unusual in that, due to steric constraints, its main-chain hydrogen-bond donors and acceptors cannot easily be satisfied. It is unable to make local hydrogen bonds, in a manner similar to that of alpha-helices, and it cannot easily satisfy the hydrogen-bonding potential of neighboring residues in polyproline conformation in a manner analogous to beta-strands. Here we describe an analysis of polyproline conformations using the HOMSTRAD database of structurally aligned proteins. This allows us not only to determine amino acid propensities from a much larger database than previously but also to investigate conservation of amino acids in polyproline conformations, and the conservation of the conformation itself. Although proline is common in polyproline helices, helices without proline represent 46% of the total. No other amino acid appears to be greatly preferred; glycine and aromatic amino acids have low propensities for PPII. Accordingly, the hydrogen-bonding potential of PPII main-chain is mainly satisfied by water molecules and by other parts of the main-chain. Side-chain to main-chain interactions are mostly nonlocal. Interestingly, the increased number of nonsatisfied H-bond donors and acceptors (as compared with alpha-helices and beta-strands) makes PPII conformers well suited to take part in protein-protein interactions.     

Crystal structure of an Fab fragment in complex with a meningococcal serosubtype antigen and a protein G domain.

Many pathogens present highly variable surface proteins to their host as a means of evading immune responses. The structure of a peptide antigen corresponding to the subtype P1.7 variant of the porin PorA from the human pathogen Neisseria meningitidis was determined by solution of the X-ray crystal structure of the ternary complex of the peptide (ANGGASGQVK) in complex with a Fab fragment and a domain from streptococcal protein G to 1.95 A resolution. The peptide adopted a beta-hairpin structure with a type I beta-turn between residues Gly4P and Gly7P, the conformation of the peptide being further stabilised by a pair of hydrogen bonds from the side-chain of Asn2P to main-chain atoms in Val9P. The antigen binding site within the Fab formed a distinct crevice lined by a high proportion of apolar amino acids. Recognition was supplemented by hydrogen bonds from heavy chain residues Thr50H, Asp95H, Leu97H and Tyr100H to main-chain and side-chain atoms in the peptide. Complementarity-determining region (CDR) 3 of the heavy chain was responsible for approximately 50 % of the buried surface area formed by peptide-Fab binding, with the remainder made up from CDRs 1 and 3 of the light chain and CDRs 1 and 2 of the heavy chain. Knowledge of the structures of variable surface antigens such as PorA is an essential prerequisite to a molecular understanding of antigenic variation and its implications for vaccine design.     

Molecular dynamics studies of side chain effect on the beta-1,3-D-glucan triple helix in aqueous solution

beta-1,3-D-glucans have been isolated from fungi as right-handed 6(1) triple helices. They are categorized by the side chains bound to the main triple helix through beta-(1-->6)-D-glycosyl linkage. Indeed, since a glucose-based side chain is water soluble, the presence and frequency of glucose-based side chains give rise to significant variation in the physical properties of the glucan family. Curdlan has no side chains and self-assembles to form an water-insoluble triple helical structure, while schizophyllan, which has a 1,6-D-glucose side chain on every third glucose unit along the main chain, is completely water soluble. A thermal fluctuation in the optical rotatory dispersion is observed for the side chain, indicating probable co-operative interaction between the side chains and water molecules. This paper documents molecular dynamics simulations in aqueous solution for three models of the beta-1,3-D-glucan series: curdlan (no side chain), schizophyllan (a beta-(1-->6)-D-glycosyl side-chain at every third position), and a hypothetical triple helix with a side chain at every sixth main-chain glucose unit. A decrease was observed in the helical pitch as the population of the side chain increased. Two types of hydrogen bonding via water molecules, the side chain/main chain and the side chain/side chain hydrogen bonding, play an important role in determination of the triple helix conformation. The formation of a one-dimensional cavity of diameter about 3.5 A was observed in the schizophyllan triple helix, while curdlan showed no such cavity. The side chain/side chain hydrogen bonding in schizophyllan and the hypothetical beta-1,3-D-glucan triple helix could cause the tilt of the main-chain glucose residues to the helix.     

A possibility of a protein-bound water molecule as the ionizable group responsible for pKe at the alkaline side in human matrix metalloproteinase 7 activity

Human matrix metalloproteinase 7 (MMP-7) activity exhibits broad bell-shaped pH profile with the acidic and alkaline pK(a) (pK(e1) and pK(e2)) values of about 4 and 10. The ionizable group for pK(e2) was assigned to Lys or Arg by thermodynamic analysis; however, no such residues are present in the active site. Hence, based on the crystal structure, we hypothesized that a water molecule bound to the main-chain nitrogen of Ala162 (W1) or the main-chain carbonyl oxygen of Pro217 (W2) is a candidate for the ionizable group for pK(e2) [Takeharu, H. et al. (2011) Biochim. Biophys. Acta 1814, 1940-1946]. In this study, we inspected this hypothesis. In the hydrolysis of (7-methoxycoumarin-4-yl)acetyl-L-Pro-L-Leu-Gly-L-Leu-[N(3)-(2,4-dinitrophenyl)-L-2,3-diaminopropionyl]-L-Ala-L-Arg-NH(2), all 19 variants, in which one of all Lys and Arg residues was replaced by Ala, retained activity, indicating that neither Lys nor Arg is the ionizable group. pK(e2) values of A162S, A162V and A162G were 9.6 ± 0.1, 9.5 ± 0.1 and 10.4 ± 0.2, respectively, different from that of wild-type MMP-7 (WT) (9.9 ± 0.1) by 0.3-0.5 pH unit, and those of P217S, P217V and P217G were 10.1 ± 0.1, 9.8 ± 0.1 and 9.7 ± 0.1, respectively, different from that of WT by 0.1-0.2 pH unit. These results suggest a possibility of W1 or W2 as the ionizable group for pK(e2).     

Arabinogalactan from Western larch, Part III: alkaline degradation revisited, with novel conclusions on molecular structure

Alkaline degradation of larch arabinogalactan (AG) involves rapid peeling of the (1→3)-galactan main chain from the reducing end. The products are the original side chains attached to galactometasaccharinic (GalMS) acids derived from main-chain residues. These products have been separated by GPC and studied by compositional, methylation and NMR analyses. Results confirm that in a typical AG molecule most main-chain residues carry a side chain on C-6. About half of these side chains are β-(1→6)-linked Galp dimers, and about a quarter are single Galp residues. The rest contain three or more residues and include most of the Ara (arabinose) found in the polysaccharide. About one-fifth of the original Ara is consumed by the peeling reaction, and Ara groups at the non-reducing end of the main chain are proposed, predominantly as arabinobiosyl groups [β-l-Arap-(1→3)-L-Araf-(1→.. In general for the larger side chains abundance decreases with size, while branching and Ara content increase with size. Most terminal Araf residues occur in three- and four-residue side chains, while most arabinobiosyl groups are found in side chains of more than three residues. The Ara probably occurs only as these monomeric or dimeric groups, and since no Ara is found in side chains smaller than three residues, this implies that there are no Ara branches attached directly to the main chain.     

Analysis of the relationship between side-chain conformation and secondary structure in globular proteins

The relationship between the preferred side-chain dihedral angles and the secondary structure of a residue was examined. The structures of 61 proteins solved to a resolution of 2.0 A (1 A = 0.1 nm) or better were analysed using a relational database to store the information. The strongest feature observed was that the chi 1 distribution for most side-chains in an alpha-helix showed an absence of the g- conformation and a shift towards the t conformation when compared to the non-alpha/beta structures. The exceptions to this tendency were for short polar side-chains that form hydrogen bonds with the main-chain which prefer g+. Shifts in the chi 1 preferences for residues in the beta-sheet were observed. Other side-chain dihedral angles (chi 2, chi 3, chi 4) were found to be influenced by the main-chain. This paper presents more accurate distributions for the side-chain dihedral angles which were obtained from the increased number of proteins determined to high resolution. The means and standard deviations for chi 1 and chi 2 angles are presented for all residues according to the secondary structure of the main-chain. The means and standard deviations are given for the most popular conformations for side-chains in which chi 3 and chi 4 rotations affect the position of C atoms.     

Solution structure of switch Arc,a mutant with 3(10) helices replacing a wild-type beta-ribbon

Adjacent N11L and L12N mutations in the antiparallel beta-ribbon of Arc repressor result in dramatic changes in local structure in which each beta-strand is replaced by a right-handed helix. The full solution structure of this "switch" Arc mutant shows that irregular 3(10) helices compose the new secondary structure. This structural metamorphosis conserves the number of main-chain and side-chain to main-chain hydrogen bonds and the number of fully buried core residues. Apart from a slight widening of the interhelical angle between alpha-helices A and B and changes in side-chain conformation of a few core residues in Arc, no large-scale structural adjustments in the remainder of the protein are necessary to accommodate the ribbon-to-helix change. Nevertheless, some changes in hydrogen-exchange rates are observed, even in regions that have very similar structures in the two proteins. The surface of switch Arc is packed poorly compared to wild-type, leading to approximately 1000A(2) of additional solvent-accessible surface area, and the N termini of the 3(10) helices make unfavorable head-to-head electrostatic interactions. These structural features account for the positive m value and salt dependence of the ribbon-to-helix transition in Arc-N11L, a variant that can adopt either the mutant or wild-type structures. The tertiary fold is capped in different ways in switch and wild-type Arc, showing how stepwise evolutionary transformations can arise through small changes in amino acid sequence.     

Conformational behavior of poly-N5-(3-hydroxypropyl)-L-glutamine in water-methanol mixtures studied by 13C Nmr and CD spectroscopy

The helix–coil transition of poly-N⁵-(3-hydroxypropyl)-L -glutamine (PHPLG) has been studied in methanol–water by CD and cmr spectroscopy. For polydisperse PHPLG, two separate peaks arising from residues in helical and random-coil conformations are observed during the transition for both main-chain carbons. These results are discussed and compared to those observed in the case of a polymer sample obtained by racemization of PHPLG in 0.1 M NaOH and of PHPLG samples of controlled molecular weight and dispersity. The dominant influence of the molecular-weight heterogeneity on the double-peak phenomenon has been verified. The linewidths and chemical shift of the ¹³C resonances are discussed in terms of side-chain–main-chain interactions and side-chain solvation.     

Refined structure of monomeric diphtheria toxin at 2.3 A resolution.

The structure of toxic monomeric diphtheria toxin (DT) was determined at 2.3 A resolution by molecular replacement based on the domain structures in dimeric DT and refined to an R factor of 20.7%. The model consists of 2 monomers in the asymmetric unit (1,046 amino acid residues), including 2 bound adenylyl 3'-5' uridine 3' monophosphate molecules and 396 water molecules. The structures of the 3 domains are virtually identical in monomeric and dimeric DT; however, monomeric DT is compact and globular as compared to the "open" monomer within dimeric DT (Bennett MJ, Choe S, Eisenberg D, 1994b, Protein Sci 3:0000-0000). Detailed differences between monomeric and dimeric DT are described, particularly (1) changes in main-chain conformations of 8 residues acting as a hinge to "open" or "close" the receptor-binding (R) domain, and (2) a possible receptor-docking site, a beta-hairpin loop protruding from the R domain containing residues that bind the cell-surface DT receptor. Based on the monomeric and dimeric DT crystal structures we have determined and the solution studies of others, we present a 5-step structure-based mechanism of intoxication: (1) proteolysis of a disulfide-linked surface loop (residues 186-201) between the catalytic (C) and transmembrane (T) domains; (2) binding of a beta-hairpin loop protruding from the R domain to the DT receptor, leading to receptor-mediated endocytosis; (3) low pH-triggered open monomer formation and exposure of apolar surfaces in the T domain, which insert into the endosomal membrane; (4) translocation of the C domain into the cytosol; and (5) catalysis by the C domain of ADP-ribosylation of elongation factor 2.     

Solution conformation of cytochrome c-551 from Pseudomonas stutzeri ZoBell determined by NMR.

1H NMR spectroscopy and solution structure computations have been used to examine ferrocytochrome c-551 from Pseudomonas stutzeri ZoBell (ATCC 14405). Resonance assignments are proposed for all main-chain and most side-chain protons. Stereospecific assignments were also made for some of the beta-methylene protons and valine methyl protons. Distance constraints were determined based upon nuclear Overhauser enhancements between pairs of protons. Dihedral angle constraints were determined from estimates of scalar coupling constants and intra-residue NOEs. Twenty structures were calculated by distance geometry and refined by energy minimization and simulated annealing on the basis of 1012 interproton distance and 74 torsion angle constraints. Both the main-chain and side-chain atoms are well defined except for two terminal residues, and some side-chain atoms located on the molecular surface. The average root mean squared deviation in the position for equivalent atoms between the 20 individual structures and the mean structure obtained by averaging their coordinates is 0.56 +/- 0.10 A for the main-chain atoms, and 0.95 +/- 0.09 A for all nonhydrogen atoms of residue 3 to 80 plus the heme group. The average structure was compared with an analogous protein, cytochrome c-551 from pseudomonas stutzeri. The main-chain folding patterns are very consistent, but there are some differences, some of which can be attributed to the loss of normally conserved aromatic residues in the ZoBell c-551.     

Analysis of protein main-chain solvation as a function of secondary structure

We have analysed the hydration of main-chain carbonyl and amide groups in 24 high-resolution well-refined protein structures as a function of the secondary structure in which these polar groups occur. We find that main-chain atoms in beta-sheets are as hydrated as those in alpha-helices, with most interactions involving "free" amide and carbonyl groups that do not participate in secondary structure hydrogen bonds. The distributions of water molecules around these non-bonded carbonyl groups reflect specific steric interactions due to the local secondary structure. Approximately 20% and 4%, respectively of bonded carbonyl and amide groups interact with solvent. These include interactions with carbonyl groups on the exposed faces of alpha-helices that have been correlated previously with bending of the helix. Water molecules interacting with alpha-helices occur mainly at the amino and carbonyl termini of the helices, in which case the solvent sites maintain the hydrogen bonding by bridging between residues i and i-3 or i-4 at the amino terminus and between i and i+3 or i+4 at the carbonyl terminus. We also see a number of solvent-mediated Ncap and Ccap interactions. The water molecules interacting with beta-sheets occur mainly at the edges, in which case they extend the sheet structure, or at the ends of strands, in which case they extend the beta-ladder. In summary, the solvent networks appear to extend the hydrogen-bonding structure of the secondary structures. In beta-turns, which usually occur at the surface of a protein, exposed amide and carbonyl groups are often hydrated, especially close to glycine residues. Occasionally water molecules form a bridge between residues i and i+3 in the turn and this may provide extra stabilization.     

Amino acid conformational preferences and solvation of polar backbone atoms in peptides and proteins

Amino acids in peptides and proteins display distinct preferences for alpha-helical, beta-strand, and other conformational states. Various physicochemical reasons for these preferences have been suggested: conformational entropy, steric factors, hydrophobic effect, and backbone electrostatics; however, the issue remains controversial. It has been proposed recently that the side-chain-dependent solvent screening of the local and non-local backbone electrostatic interactions primarily determines the preferences not only for the alpha-helical but also for all other main-chain conformational states. Side-chains modulate the electrostatic screening of backbone interactions by excluding the solvent from the vicinity of main-chain polar atoms. The deficiency of this electrostatic screening model of amino acid preferences is that the relationships between the main-chain electrostatics and the amino acid preferences have been demonstrated for a limited set of six non-polar amino acid types in proteins only. Here, these relationships are determined for all amino acid types in tripeptides, dekapeptides, and proteins. The solvation free energies of polar backbone atoms are approximated by the electrostatic contributions calculated by the finite difference Poisson-Boltzmann and the Langevin dipoles methods. The results show that the average solvation free energy of main-chain polar atoms depends strongly on backbone conformation, shape of side-chains, and exposure to solvent. The equilibrium between the low-energy beta-strand conformation of an amino acid (anti-parallel alignment of backbone dipole moments) and the high-energy alpha conformation (parallel alignment of backbone dipole moments) is strongly influenced by the solvation of backbone polar atoms. The free energy cost of reaching the alpha conformation is by approximately 1.5 kcal/mol smaller for residues with short side-chains than it is for the large beta-branched amino acid residues. This free energy difference is comparable to those obtained experimentally by mutation studies and is thus large enough to account for the distinct preferences of amino acid residues. The screening coefficients gamma(local)(r) and gamma(non-local)(r) correlate with the solvation effects for 19 amino acid types with the coefficients between 0.698 to 0.851, depending on the type of calculation and on the set of point atomic charges used. The screening coefficients gamma(local)(r) increase with the level of burial of amino acids in proteins, converging to 1.0 for the completely buried amino acid residues. The backbone solvation free energies of amino acid residues involved in strong hydrogen bonding (for example: in the middle of an alpha-helix) are small. The hydrogen bonded backbone is thus more hydrophobic than the peptide groups in random coil. The alpha-helix forming preference of alanine is attributed to the relatively small free energy cost of reaching the high-energy alpha-helix conformation. These results confirm that the side-chain-dependent solvent screening of the backbone electrostatic interactions is the dominant factor in determining amino acid conformational preferences.     

Correlation between calculated local stability and hydrogen exchange rates in proteins

The attempt is made to find new correlations between local structural characteristics of proteins and the hydrogen exchange rates of their individual main-chain amides, and to relate such correlations to possible mechanisms of hydrogen exchange. It is found that in bovine pancreatic trypsin inhibitor (BPTI) the surface area buried by a particular residue and its neighbors correlates with the exchange rate of the main-chain amide of that residue. As the area buried by a particular fragment can be associated with the stabilization of the protein structure by this fragment, the correlation suggests a role for the energetics of the local unfolding in the mechanism of hydrogen exchange. Calculations based on the assumption that the exchange mechanism involves local unfolding lead to quantitative agreement between the calculated and experimentally measured exchange rates for 80% of the amides of BPTI that are buried or hydrogen bonded to the main-chain or to internal water molecules. The same degree of correlation is found between the calculated exchange rates and partial exchange data for ribonuclease S, hen lysozyme and cytochrome c. A similarly strong correlation is found between calculated exchange rates and the exchange rates of ribonuclease A determined by neutron diffraction in the crystal. The criteria of correlation are, however, less stringent in this case because of the experimental errors, which are larger than for solution data. It is suggested that the observed correlation be used for predictions of hydrogen exchange rates in proteins.