Backbone makes a significant contribution to the electrostatics of alpha/beta-barrel proteins.

The electrostatic properties of seven alpha/beta-barrel enzymes selected from different evolutionary families were studied: triose phosphate isomerase, fructose-1,6-bisphosphate aldolase, pyruvate kinase, mandelate racemase, trimethylamine dehydrogenase, glycolate oxidase, and narbonin, a protein without any known enzymatic activity. The backbone of the alpha/beta-barrel has a distinct electrostatic field pattern, which is dipolar along the barrel axis. When the side chains are included in the calculations the general effect is to modulate the electrostatic pattern so that the electrostatic field is generally enhanced and is focused into a specific area near the active site. We use the electrostatic flux through a square surface near the active site to gauge the functionally relevant magnitude of the electrostatic field. The calculations reveal that in six out of the seven cases the backbone itself contributes greater than 45% of the total flux. The substantial electrostatic contribution of the backbone correlates with the known preference of alpha/beta-barrel enzymes for negatively charged substrates.     

AhpF and other NADH:peroxiredoxin oxidoreductases, homologues of low Mr thioredoxin reductase.

A group of bacterial flavoproteins related to thioredoxin reductase contain an additional approximately 200-amino-acid domain including a redox-active disulfide center at their N-termini. These flavoproteins, designated NADH:peroxiredoxin oxidoreductases, catalyze the pyridine-nucleotide-dependent reduction of cysteine-based peroxidases (e.g. Salmonella typhimurium AhpC, a member of the peroxiredoxin family) which in turn reduce H2O2 or organic hydroperoxides. These enzymes catalyze rapid electron transfer (kcat > 165 s-1) through one tightly bound FAD and two redox-active disulfide centers, with the N-terminal-most disulfide center acting as a redox mediator between the thioredoxin-reductase-like part of these proteins and the peroxiredoxin substrates. A chimeric protein with the first 207 amino acids of S. typhimurium AhpF attached to the N-terminus of Escherichia coli thioredoxin reductase exhibits very high NADPH:peroxiredoxin oxidoreductase and thioredoxin reductase activities. Catalytic turnover by NADH:peroxiredoxin oxidoreductases may involve major domain rotations, analogous to those proposed for bacterial thioredoxin reductase, and cycling of these enzymes between two electron-reduced (EH2) and four electron-reduced (EH4) redox states.     

Two-dimensional NMR and photo-CIDNP studies of the insulin monomer: assignment of aromatic resonances with application to protein folding, structure, and dynamics

The aromatic 1H NMR resonances of the insulin monomer are assigned at 500 MHz by comparative studies of chemically modified and genetically altered variants, including a mutant insulin (PheB25----Leu) associated with diabetes mellitus. The two histidines, three phenylalanines, and four tyrosines are observed to be in distinct local environments; their assignment provides sensitive markers for studies of tertiary structure, protein dynamics, and protein folding. The environments of the tyrosine residues have also been investigated by photochemically induced dynamic nuclear polarization (photo-CIDNP) and analyzed in relation to packing constraints in the crystal structures of insulin. Dimerization involving specific B-chain interactions is observed with increasing protein concentration and is shown to depend on temperature, pH, and solvent composition. In the monomer large variations are observed in the line widths of amide resonances, suggesting intermediate exchange among conformational substates; such substates may relate to conformational changes observed in different crystal states and proposed to occur in the hormone-receptor complex. Additional evidence for multiple conformations in solution is provided by comparative studies of an insulin analogue containing a peptide bond between residues B29 and A1 (mini-proinsulin). This analogue forms dimers and higher-order oligomers under conditions in which native insulin is monomeric, suggesting that the B29-A1 peptide bond stabilizes a conformational substate favorable for dimerization. Such stabilization is not observed in corresponding studies of native proinsulin, in which a 35-residue connecting peptide joins residues B30 and A1; this extended tether is presumably too flexible to constrain the conformation of the B-chain. The differences between proinsulin and mini-proinsulin suggest a structural mechanism for the observation that the fully reduced B29-A1 analogue folds more efficiently than proinsulin to form the correct pattern of disulfide bonds. These results are discussed in relation to molecular mechanics calculations of insulin based on the available crystal structures.     

Molecular dynamics simulations of native and substrate-bound lysozyme. A study of the average structures and atomic fluctuations

Molecular dynamics simulations of hen egg-white lysozyme in the free and substrate-bound states are reported and the nature of the average structures and atomic fluctuations are analyzed. Crystallographic water molecules of structural importance, as determined by hydrogen-bonding, were included in the simulations. Comparisons are made between the dynamics and the X-ray results for the atomic positions, the main-chain and side-chain dihedral angles, and the hydrogen-bonding geometry. Improvements over earlier simulations in the potential energy function and methodology resulted in stable trajectories with the C alpha co-ordinates within 1.5 A of the starting X-ray structure. Structural features analyzed in the simulations agreed well with the X-ray results except for some surface residues. The Asx chi 2 dihedral distribution and the geometry of hydrogen bonding at reverse turns show differences; possible causes are discussed. The relation between the magnitudes and time-scales of the residue fluctuations and secondary structural features, such as helices beta-sheets and coiled loops, is examined. Significant differences in the residue mobilities between the simulations of the free and substrate-bound states were found in a region of the enzyme that is in direct contact with the substrate and in a region that is distant from the active-site cleft. The dynamic behavior of the structural water molecules is analyzed by examining the correlation between the fluctuations of the water oxygens and the lysozyme heavy-atoms to which they are hydrogen-bonded.     

Application of molecular dynamics with interproton distance restraints to three-dimensional protein structure determination. A model study of crambin

The applicability of restrained molecular dynamics for the determination of three-dimensional protein structures on the basis of short interproton distances (less than 4 A) that can be realistically determined from nuclear magnetic resonance measurements in solution is assessed. The model system used is the 1.2 A resolution crystal structure of the 46 residue protein crambin, from which a set of 240 approximate distance restraints, divided into three ranges (2.5 +/- 0.5, 3.0+0.5(-1.0) and 4 +/- 1 A), is derived. This interproton distance set comprises 159 short-range ([i-j] less than or equal to 5) and 56 ([i-j] greater than 5) long-range inter-residue distances and 25 intra-residue distances. Restrained molecular dynamics are carried out using a number of different protocols starting from two initial structures: a completely extended beta-strand; and an extended structure with two alpha-helices in the same positions as in the crystal structure (residues 7 to 19, and 23 to 30) and all other residues in the form of extended beta-strands. The root-mean-square (r.m.s.) atomic differences between these two initial structures and the crystal structure are 43 A and 23 A, respectively. It is shown that, provided protocols are used that permit the secondary structure elements to form at least partially prior to folding into a tertiary structure, convergence to the correct final structure, both globally and locally, is achieved. The r.m.s. atomic differences between the converged restrained dynamics structures and the crystal structure range from 1.5 to 2.2 A for the backbone atoms and from 2.0 to 2.8 A for all atoms. The r.m.s. atomic difference between the X-ray structure and the structure obtained by first averaging the co-ordinates of the converged restrained dynamics structures is even smaller: 1.0 A for the backbone atoms and 1.6 A for all atoms. These results provide a measure with which to judge future experimental results on proteins whose crystal structures are unknown. In addition, from an examination of the dynamics trajectories, it is shown that the convergence pathways followed by the various simulations are different.     

The three-dimensional structure of alpha1-purothionin in solution: combined use of nuclear magnetic resonance, distance geometry and restrained molecular dynamics

The determination of the three-dimensional solution structure of α1-purothionin using a combination of metric matrix distance geometry and restrained molecular dynamics calculations based on n.m.r. data is presented. The experimental data comprise complete sequence-specific proton resonance assignments, a set of 310 approximate interproton distance restraints derived from nuclear Overhauser effects, 27 Ø backbone torsion angle restraints derived from vicinal coupling constants, 4 distance restraints from hydrogen bonds and 12 distance restraints from disulphide bridges. The average atomic rms difference between the final nine converged structures and the mean structure obtained by averaging their coordinates is 1.5 ± 0.1 å for the backbone atoms and 2.0 ± 0.1 å for all atoms. The overall shape of α1-purothionin is that of the capital letter L, similar to that of crambin, with the longer arm comprising two approximately parallel α-helices and the shorter arm a strand and a mini anti-parallel β sheet.     

Hemoglobin tertiary structural change on ligand binding its role in the co-operative mechanism

Analysis of the tertiary structural alterations in hemoglobin induced by ligand binding demonstrates that an allosteric core composed of the heme, histidine F8, the FG corner and part of the F-helix plays an essential role in co-operativity. This conclusion is based on structural and spectroscopic data and theoretical studies of hemoglobin chains. The methodology employed in the calculations is presented with details of the empirical energy function. Energy minimized structures of the unliganded hemoglobin chains, which serve as reference systems for the analysis, are described. To determine the structural changes induced by ligand binding, the effects of FeN bond shortening and of heme translation and tilting perturbations are examined. Energy minimization in the presence of the perturbations serves to provide information concerning the globin structural modifications produced by them. The validity of the results is supported by comparisons with the X-ray data of Anderson, Pulsinelli, Baldwin and Chothia on tertiary changes in the hemoglobin subunits.Internal to the allosteric core, there appear to be two stable positions for its elements: one of these corresponds to the liganded and the other to the unliganded species. The unliganded geometry fits without strain into the deoxy tetramer, while the liganded one fits without strain into the oxy tetramer. On ligation of a subunit in the deoxy tetramer, the structural changes within the allosteric core are in the direction of those found in going from the unliganded deoxy to the liganded oxy system, although they are reduced by the presence of constraints due to the other subunits in the deoxy tetramer. In addition, the quaternary constraints in the deoxy tetramer prevent the large overall displacement of the allosteric core that occurs in the transition to the liganded oxy tetramer. The coupling between the changes internal to the allosteric core, produced on ligation and the overall displacement of the core that accompanies the quaternary transition, is an essential element of the co-operative mechanism. As shown in previous work (Gelin & Karplus, 1977), the proximal histidine serves as the link between the position of the heme and the F-helix; the asymmetric orientation of the histidine in the deoxy structure, coupled with contributions from other heme-protein interactions, appears to initiate the tertiary structural changes induced by ligand binding. The reduced oxygen affinity of hemoglobin results not from tension on the heme in the unliganded structure (there is none) but instead from strain in the liganded subunit of the tetramer within the deoxy quaternary structure. Further, the changes in the allosteric core provide a relatively localized reaction path for transmitting information concerning ligand binding from the heme group to the surface of the subunit; particularly in the α-chain, the residue Val FG5 appears to play an important role in the reaction path.The present analysis has important implications for realistic statistical thermodynamic models of hemoglobin co-operativity. It suggests that the previously formulated model (Szabo & Karplus, 1972) should be generalized by the introduction of two different subunit tertiary structures in the deoxy and in the oxy tetramer; they would be associated with the unliganded and the liganded allosteric core, respectively, and would take account of steric constraints that reduce the ligand affinity of the deoxy tetramer.     

Helix–coil transitions in a simple polypeptide model

A simplified model of a polypeptide chain is described. Each residue is represented by a single interaction center. The energy of the chain and the force acting on each residue are given as a function of the residue coordinates. Terms to approximate the effect of solvent and the stabilization energy of helix formation are included. The model is used to study equilibrium and dynamical aspects of the helix–coil transition. The equilibrium properties examined include helix–coil equilibrium constants and their dependence on chain position. Dynamical properties are examined by a stochastic simulation of the Brownian motion of the chain in its solvent surroundings. Correlations in the motions of the residues are found to have an important influence on the helix–coil transition rates.     

Molecular dynamics of ferrocytochrome c. Magnitude and anisotropy of atomic displacements

Two computer simulations of the atomic motion in tuna ferrocytochrome c have been carried out. The average structures and the structural correlations of the magnitudes of the atomic position fluctuations are in substantial agreement with recent X-ray diffraction results, particularly for the protein interior. The simulations show, however, that the atomic displacements are quite anisotropic. The degree of anisotropy and the preferred directions of atomic displacement exhibit correlations with structural features of the protein.     

A new two-disulphide intermediate in the refolding of reduced bovine pancreatic trypsin inhibitor

The apparently complete refolding of reduced bovine pancreatic trypsin inhibitor (BPTI) is shown to produce a mixture of two species. One of these is native BPTI, but the other lacks the disulphide bond between cysteines 30 and 51. The latter species has a folded conformation very like that of native BPTI, and is oxidized by air to native BPTI on warming in aqueous solution. The two unreactive cysteine thiol groups appear to be buried in the interior of the molecule, which restricts access by reagents that can alkylate them or oxidize them to form the disulphide bond. The implications of this intermediate and its conformation for the understanding of protein folding are discussed.