Domain closure mechanism in transferrins: new viewpoints about the hinge structure and motion as deduced from high resolution crystal structures of ovotransferrin N-lobe

The crystal structure of holo hen ovotransferrin N-lobe refined at 1.65 A resolution has been obtained. The final model gave an R-factor of 0.173 in the resolution range between 10.0 and 1.65 A. The comparison of the structure with previous high-resolution apo and Fe(3+)-loaded, domain-opened intermediate structures provides new viewpoints on the domain closure mechanism upon Fe(3+) uptake in ovotransferrin N-lobe. Overall, conformational transition follows the common mechanism that has been first demonstrated for lactoferrin N-lobe; the domains 1 and 2 rotate 49.7 degrees as rigid bodies with a translation of 2.1 A around a screw-axis that passes through the two interdomain beta-strands (89-94 and 244-249). It is generally believed that the two strands display a hinge-like motion. Here, the latter strand indeed displays an ideal hinge nature: the segments 244-246 and 248-249 behave as a part of the rigid body of domain 2 and that of domain 1, respectively, and a sharp bend upon the domain closure is largely accounted for by the changes in the torsion angles phi and psi of Val247. We find, however, that the mode of the conformational change in the first beta-strand is much more complex. Two of the five inter beta-strand hydrogen bonds undergo crucial exchanges: from Ser91-N...Val247-O and Thr89-O...Ala249-N in the open apo and intermediate structures into Tyr92-N...Val247-O and Thr90-O...Ala249-N in the closed holo structure. These exchanges, which may be triggered in the intermediate state by modulation in the topological relation between the Fe(3+)-ligated hinge residue Tyr92-OH and the anion anchor residues of helix 5, are accompanied by a large conformational change and extensive hydrogen bonding rearrangements in a long stretch of segment of Glu82 to Tyr92. Such structural transition would work as a driving force for the domain closure, which highlights a "door closer"-like role, in addition to the canonical-hinge role, for the interdomain polypeptide segment pair. As an alternative hinge that secures the correct domain motion by being placed on a significant distance from the beta-strand hinge, we point out the participation of the van der Waals contacts formed between domain 1 residue of Met331 and domain 2 residues of Trp125, Ile129 and Trp140.     

The Effect of End Constraints on Protein Loop Kinematics

Despite the prevalent involvement of loops in function little is known about how the constraining of end groups influences their kinematics. Using a linear inverse-kinematics approach and assuming fixed bond lengths, bond angles, and peptide bond torsions, as well as ignoring molecular interactions to assess the effect of the end-constraint only, it is shown that the constraint creates a closed surface in torsion angle space. For pentapeptides, the constraint gives rise to inaccessible regions in a Ramachandran plot. This complex and tightly curved surface produces interesting effects that may play a functional role. For example, a small change in one torsion angle can radically change the behavior of the whole loop. The constraint also produces long-range correlations, and structures exist where the correlation coefficient is 1.0 or −1.0 between rotations about bonds separated by >30 Å. Another application allows some torsion angles to be targeted to specified values while others are constrained. When this application was used on key torsions in lactate dehydrogenase, it was found that the functional loop first folds forward and then moves sideways. For horse liver alcohol dehydrogenase, it was confirmed that the functional loop's Pro-Pro motif creates a rigid arm in an NAD-activated switch for domain closure.     

Efficient search on energy minima for structure prediction of nucleic acid motifs

Structure prediction of non-canonical motifs such as mismatches, extra unmatched nucleotides or internal and hairpin loop structures in nucleic acids is of great importance for understanding the function and design of nucleic acid structures. Systematic conformational analysis of such motifs typically involves the generation of many possible combinations of backbone dihedral torsion angles for a given motif and subsequent energy minimization (EM) and evaluation. Such approach is limited due to the number of dihedral angle combinations that grows very rapidly with the size of the motif. Two conformational search approaches have been developed that allow both an effective crossing of barriers during conformational searches and the computational demand grows much less with system size then search methods that explore all combinations of backbone dihedral torsion angles. In the first search protocol single torsion angles are flipped into favorable states using constraint EM and subsequent relaxation without constraints. The approach is repeated in an iterative manner along the backbone of the structural motif until no further energy improvement is obtained. In case of two test systems, a DNA-trinucleotide loop (sequence: GCA) and a RNA tetraloop (sequence: UUCG), the approach successfully identified low energy states close to experiment for two out of five start structures. In the second method randomly selected combinations of up to six backbone torsion angles are simultaneously flipped into preset ranges by a short constraint EM followed by unconstraint EM and acceptance according to a Metropolis acceptance criterion. This combined stochastic/EM search was even more effective than the single torsion flip approach and selected low energy states for the two test cases in between two and four cases out of five start structures.     

Comparison of simulated and experimentally determined dynamics for a variant of the Lacl DNA-binding domain, Nlac-P.

Recent advances in the experimentally determined structures and dynamics of the domains within LacI provide a rare context for evaluating dynamics calculations. A 1500-ps trajectory was simulated for a variant of the LacI DNA-binding domain, which consists of the first three helices in LacI and the hinge helix of the homologous PurR. Order parameters derived from dynamics simulations are compared to those obtained for the LacI DNA-binding domain with 15N relaxation NMR spectroscopy (Slijper et al., 1997. Biochemistry. 36:249-254). The MD simulations suggest that the unstructured loop between helices II and III does not exist in a discrete state under the conditions of no salt and neutral pH, but occupies a continuum of states between the DNA-bound and free structures. Simulations also indicate that the unstructured region between helix III and the hinge helix is very mobile, rendering motions of the hinge helix essentially independent of the rest of the protein. Finally, the alpha-helical hydrogen bonds in the hinge helix are broken after 1250 ps, perhaps as a prelude to helix unfolding.     

Structural basis for a change in substrate specificity: crystal structure of S113E isocitrate dehydrogenase in a complex with isopropylmalate, Mg2+, and NADP

Isocitrate dehydrogenase (IDH) catalyzes the oxidative decarboxylation of isocitrate and has negligible activity toward other (R)-malate-type substrates. The S113E mutant of IDH significantly improves its ability to utilize isopropylmalate as a substrate and switches the substrate specificity (k(cat)/K(M)) from isocitrate to isopropylmalate. To understand the structural basis for this switch in substrate specificity, we have determined the crystal structure of IDH S113E in a complex with isopropylmalate, NADP, and Mg(2+) to 2.0 A resolution. On the basis of a comparison with previously determined structures, we identify distinct changes caused by the amino acid substitution and by the binding of substrates. The S113E complex exhibits alterations in global and active site conformations compared with other IDH structures that include loop and helix conformational changes near the active site. In addition, the angle of the hinge that relates the two domains was altered in this structure, which suggests that the S113E substitution and the binding of substrates act together to promote catalysis of isopropylmalate. Ligand binding results in reorientation of the active site helix that contains residues 113 through 116. E113 exhibits new interactions, including van der Waals contacts with the isopropyl group of isopropylmalate and a hydrogen bond with N115, which in turn forms a hydrogen bond with NADP. In addition, the loop and helix regions that bind NADP are altered, as is the loop that connects the NADP binding region to the active site helix, changing the relationship between substrates and enzyme. In combination, these interactions appear to provide the basis for the switch in substrate specificity.     

Efficient Search on Energy Minima for Structure Prediction of Nucleic Acid Motifs

Abstract

Structure prediction of non-canonical motifs such as mismatches, extra unmatched nucleotides or internal and hairpin loop structures in nucleic acids is of great importance for understanding the function and design of nucleic acid structures. Systematic conformational analysis of such motifs typically involves the generation of many possible combinations of backbone dihedral torsion angles for a given motif and subsequent energy minimization (EM) and evaluation. Such approach is limited due to the number of dihedral angle combinations that grows very rapidly with the size of the motif. Two conformational search approaches have been developed that allow both an effective crossing of barriers during con-formational searches and the computational demand grows much less with system size then search methods that explore all combinations of backbone dihedral torsion angles. In the first search protocol single torsion angles are flipped into favorable states using constraint EM and subsequent relaxation without constraints. The approach is repeated in an iterative manner along the backbone of the structural motif until no further energy improvement is obtained. In case of two test systems, a DNA-trinucleotide loop (sequence: GCA) and a RNA tetraloop (sequence: UUCG), the approach successfully identified low energy states close to experiment for two out of five start structures. In the second method randomly selected combinations of up to six backbone torsion angles are simultaneously flipped into preset ranges by a short constraint EM followed by unconstraint EM and acceptance according to a Metropolis acceptance criterion. This combined stochastic/EM search was even more effective than the single torsion flip approach and selected low energy states for the two test cases in between two and four cases out of five start structures.     

Conformational features of a hexapeptide model Ac-TGAAKA-NH2 corresponding to a hydrated alpha helical segment from glyceraldehyde 3-phosphate dehydrogenase: implications for the role of turns in helix folding

Recent analysis of alpha helices in protein crystal structures, available in literature, revealed hydrated alpha helical segments in which, water molecule breaks open helix 5-->1 hydrogen bond by inserting itself, hydrogen bonds to both C=O and NH groups of helix hydrogen bond without disrupting the helix hydrogen bond, and hydrogen bonds to either C=O or NH of helix hydrogen bond. These hydrated segments display a variety of turn conformations and are thought to be 'folding intermediates' trapped during folding-unfolding of alpha helices. A role for reverse turns is implicated in the folding of alpha helices. We considered a hexapeptide model Ac-1TGAAKA6-NH2 from glyceraldehyde 3-phosphate dehydrogenase, corresponding to a hydrated helical segment to assess its role in helix folding. The sequence is a site for two 'folding intermediates'. The conformational features of the model peptide have been investigated by 1H 2D NMR techniques and quantum mechanical perturbative configuration interaction over localized orbitals (PCILO) method. Theoretical modeling largely correlates with experimental observations. Based upon the amide proton temperature coefficients, the observed d alpha n(i, i + 1), d alpha n(i, i + 2), dnn(i, i + 1), d beta n(i, i + 1) NOEs and the results from theoretical modeling, we conclude that the residues of the peptide sample alpha helical and neck regions of the Ramachandran phi, psi map with reduced conformational entropy and there is a potential for turn conformations at N and C terminal ends of the peptide. The role of reduced conformational entropy and turn potential in helix formation have been discussed. We conclude that the peptide sequence can serve as a 'folding intermediate' in the helix folding of glyceraldehyde 3-phosphate dehydrogenase.     

Protein flexibility predictions using graph theory

Techniques from graph theory are applied to analyze the bond networks in proteins and identify the flexible and rigid regions. The bond network consists of distance constraints defined by the covalent and hydrogen bonds and salt bridges in the protein, identified by geometric and energetic criteria. We use an algorithm that counts the degrees of freedom within this constraint network and that identifies all the rigid and flexible substructures in the protein, including overconstrained regions (with more crosslinking bonds than are needed to rigidify the region) and underconstrained or flexible regions, in which dihedral bond rotations can occur. The number of extra constraints or remaining degrees of bond-rotational freedom within a substructure quantifies its relative rigidity/flexibility and provides a flexibility index for each bond in the structure. This novel computational procedure, first used in the analysis of glassy materials, is approximately a million times faster than molecular dynamics simulations and captures the essential conformational flexibility of the protein main and side-chains from analysis of a single, static three-dimensional structure. This approach is demonstrated by comparison with experimental measures of flexibility for three proteins in which hinge and loop motion are essential for biological function: HIV protease, adenylate kinase, and dihydrofolate reductase.     

Crystal structure of the complex of plasminogen activator inhibitor 2 with a peptide mimicking the reactive center loop

The structure of the serpin, plasminogen activator inhibitor type-2 (PAI-2), in a complex with a peptide mimicking its reactive center loop (RCL) has been determined at 1.6-A resolution. The structure shows the relaxed state serpin structure with a prominent six-stranded beta-sheet. Clear electron density is seen for all residues in the peptide. The P1 residue of the peptide binds to a well defined pocket at the base of PAI-2 that may be important in determining the specificity of protease inhibition. The stressed-to-relaxed state (S --> R) transition in PAI-2 can be modeled as the relative motion between a quasirigid core domain and a smaller segment comprising helix hF and beta-strands s1A, s2A, and s3A. A comparison of the Ramachandran plots of the stressed and relaxed state PAI-2 structures reveals the location of several hinge regions connecting these two domains. The hinge regions cluster in three locations on the structure, ensuring a cooperative S --> R transition. We hypothesize that the hinge formed by the conserved Gly(206) on beta-strand s3A in the breach region of PAI-2 effects the S --> R transition by altering its backbone torsion angles. This torsional change is due to the binding of the P14 threonine of the RCL to the open breach region of PAI-2.     

An Allosteric Signaling Pathway of Human 3-Phosphoglycerate Kinase from Force Distribution Analysis

3-Phosphogycerate kinase (PGK) is a two domain enzyme, which transfers a phosphate group between its two substrates, 1,3-bisphosphoglycerate bound to the N-domain and ADP bound to the C-domain. Indispensable for the phosphoryl transfer reaction is a large conformational change from an inactive open to an active closed conformation via a hinge motion that should bring substrates into close proximity. The allosteric pathway resulting in the active closed conformation has only been partially uncovered. Using Molecular Dynamics simulations combined with Force Distribution Analysis (FDA), we describe an allosteric pathway, which connects the substrate binding sites to the interdomain hinge region. Glu192 of alpha-helix 7 and Gly394 of loop L14 act as hinge points, at which these two secondary structure elements straighten, thereby moving the substrate-binding domains towards each other. The long-range allosteric pathway regulating hPGK catalytic activity, which is partially validated and can be further tested by mutagenesis, highlights the virtue of monitoring internal forces to reveal signal propagation, even if only minor conformational distortions, such as helix bending, initiate the large functional rearrangement of the macromolecule.     

The structural coupling between ATPase activation and recovery stroke in the myosin II motor

Before the myosin motor head can perform the next power stroke, it undergoes a large conformational transition in which the converter domain, bearing the lever arm, rotates approximately 65 degrees . Simultaneous with this "recovery stroke," myosin activates its ATPase function by closing the Switch-2 loop over the bound ATP. This coupling between the motions of the converter domain and of the 40 A-distant Switch-2 loop is essential to avoid unproductive ATP hydrolysis. The coupling mechanism is determined here by finding a series of optimized intermediates between crystallographic end structures of the recovery stroke (Dictyostelium discoideum), yielding movies of the transition at atomic detail. The successive formation of two hydrogen bonds by the Switch-2 loop is correlated with the successive see-saw motions of the relay and SH1 helices that hold the converter domain. SH1 helix and Switch-2 loop communicate via a highly conserved loop that wedges against the SH1-helix upon Switch-2 closing.     

NMR solution structure of the non-RGD disintegrin obtustatin

The solution structure of obtustatin, a novel non-RGD disintegrin of 41 residues isolated from Vipera lebetina obtusa venom, and a potent and selective inhibitor of the adhesion of integrin alpha(1)beta(1) to collagen IV, has been determined by two-dimensional nuclear magnetic resonance. Almost the whole set of chemical shifts for 1H, 13C and 15N were assigned at natural abundance from 2D homonuclear and heteronuclear 500 MHz, 600 MHz and 800 MHz spectra at pH 3.0 recorded at 298 K and 303 K. Final structural constraints consisted of 302 non-redundant NOE (95 long-range, 60 medium, 91 sequential and 56 intra-residue), four disulfide bond distances, five chi1 dihedral angles and four hydrogen bonds. The 20 conformers with lowest total energy had no NOE violations greater than 0.35A or dihedral angle violations greater than 12 degrees. The average root-mean-square deviation (RMSD) for backbone atoms of all residues among the 20 conformers was 1.1A and 0.6A for the 29 best-defined residues. Obtustatin lacks any secondary structure. Compared to all known disintegrin structures in which the RGD motif is located at the apex of an 11 residue hairpin loop, the active KTS tripeptide of obtustatin is oriented towards a side of its nine residue integrin-binding loop. The C-terminal tail is near to the active loop, and these two structural elements display the largest atomic displacements due to local conformational disorder. Double cross-peaks for W20, Y28 and H27 in the aromatic region of TOCSY spectra, local RMSD values for these residues, and positive cross-peaks in a ROESY spectrum (600 MHz, 100 ms mixing time), suggest that these residues act as a hinge allowing for the overall flexibility of the entire integrin-binding loop. These distinct structural features, along with its different electrostatic surface potential in relation to other known disintegrins, may confer to obtustatin its reported alpha(1)beta(1) integrin inhibitory selectivity.     

Hsc70-interacting HPD loop of the J domain of polyomavirus T antigens fluctuates in ps to ns and micros to ms

The backbone dynamics of the J domain from polyomavirus T antigens have been investigated using 15N NMR relaxation and molecular dynamics simulation. Model-free relaxation analysis revealed picosecond to nanosecond motions in the N terminus, the I-II loop, the C-terminal end of helix II through the HPD loop to the beginning of helix III, and the C-terminal end of helix III to the C terminus. The backbone dynamics of the HPD loop and termini are dominated by motions with moderately large amplitudes and correlation times of the order of a nanosecond or longer. Conformational exchange on the microsecond to millisecond timescale was identified in the HPD loop, the N and C termini, and the I-II loop. A 9.7ns MD trajectory manifested concerted swings of the HPD loop. Transitions between major and minor conformations of the HPD loop featured distinct patterns of change in backbone dihedral angles and hydrogen bonds. Fraying of the C-terminal end of helix II and the N-terminal end of helix III correlated with displacements of the HPD loop. Correlation of crankshaft motions of Gly46 and Gly47 with the collective motions of the HPD loop suggested an important role of the two glycine residues in the mobility of the loop. Fluctuations of the HPD loop correlated with relative reorientation of side-chains of Lys35 and Asp44 that interact with Hsc70.     

Visualisation of a 2'-5' parallel stranded double helix at atomic resolution: crystal structure of cytidylyl-2',5'-adenosine

X-ray crystallographic studies on 3'-5' oligomers have provided a great deal of information on the stereochemistry and conformational flexibility of nucleic acids and polynucleotides. In contrast, there is very little information available on 2'-5' polynucleotides. We have now obtained the crystal structure of Cytidylyl-2',5'-Adenosine (C2'p5'A) at atomic resolution to establish the conformational differences between these two classes of polymers. The dinucleoside phosphate crystallises in the monoclinic space group C2, with a = 33.912(4)A, b = 16.824(4)A, c = 12.898(2)A and beta = 112.35(1) with two molecules in the asymmetric unit. Spectacularly, the two independent C2'p5'A molecules in the asymmetric unit form right handed miniature parallel stranded double helices with their respective crystallographic two fold (b axis) symmetry mates. Remarkably, the two mini duplexes are almost indistinguishable. The cytosines and adenines form self-pairs with three and two hydrogen bonds respectively. The conformation of the C and A residues about the glycosyl bond is anti same as in the 3'-5' analog but contrasts the anti and syn geometry of C and A residues in A2'p5'C. The furanose ring conformation is C3' endo, C2' endo mixed puckering as in the C3'p5'A-proflavine complex. A comparison of the backbone torsion angles with other 2'-5' dinucleoside structures reveals that the major deviations occur in the torsion angles about the C3'-C2' and C4'-C3' bonds. A right-handed 2'-5' parallel stranded double helix having eight base pairs per turn and 45 degrees turn angle between them has been constructed using this dinucleoside phosphate as repeat unit. A discussion on 2'-5' parallel stranded double helix and its relevance to biological systems is presented.     

Uniformity,ideality, and hydrogen bonds in transmembrane alpha-helices

Protein environments substantially influence the balance of molecular interactions that generate structural stability. Transmembrane helices exist in the relatively uniform low dielectric interstices of the lipid bilayer, largely devoid of water and with a very hydrophobic distribution of amino acid residues. Here, through an analysis of bacteriorhodopsin crystal structures and the transmembrane helix structure from M2 protein of influenza A, some helices are shown to be exceptionally uniform in hydrogen bond geometry, peptide plane tilt angle, and backbone torsion angles. Evidence from both the x-ray crystal structures and solid-state NMR structure suggests that the intramolecular backbone hydrogen bonds are shorter than their counterparts in water-soluble proteins. Moreover, the geometry is consistent with a dominance of electrostatic versus covalent contributions to these bonds. A comparison of structure as a function of resolution shows that as the structures become better characterized the helices become much more uniform, suggesting that there is a possibility that many more uniform helices will be observed, even among the moderate resolution membrane protein structures that are currently in the Protein Data Bank that do not show such features.     

Ribonuclease A. Analysis of the hydrogen bond geometry, and spatial accessibility at the active site

The hydrogen bonding of bovine ribonuclease A derived from the high resolution X-ray structure has been studied in detail. Correlations have been examined for main-chain-main-chain hydrogen bond angles, torsion angles and distances, respectively. Differences are found consistently for correlations associated with alpha-helix and beta-sheet, respectively. Ten of the 124 side-chains have four or more hydrogen bond contacts; two, including Glu-101, have five or more. Three potential C = O---H, three N---X and three potential side-chain H-bonds fail to form. A search for highly inaccessible buried residues resulted in nine outstanding examples, all of which are conserved across 38 known mammalian ribonuclease A sequences, indicating the importance of these residues for structural stability. Of the two histidines in the active site, His-12 has five hydrogen bonds and His-119 three. The conformational space accessible to these two catalytically important residues studied by means of simple non-bonded contact energy calculations confirms the existence of two alternative, interchangeable locations for His-119, while His-12 is locked in a local energy minimum.     

Conformational flexibility in the allosteric regulation of human UDP-α-D-glucose 6-dehydrogenase

UDP-α-D-xylose (UDX) acts as a feedback inhibitor of human UDP-α-D-glucose 6-dehydrogenase (hUGDH) by activating an unusual allosteric switch, the Thr131 loop. UDX binding induces the Thr131 loop to translate ~5 ? through the protein core, changing packing interactions and rotating a helix (α6(136-144)) to favor the formation of an inactive hexameric complex. But how does to conformational change occur given the steric packing constraints of the protein core? To answer this question, we deleted Val132 from the Thr131 loop to approximate an intermediate state in the allosteric transition. The 2.3 ? resolution crystal structure of the deletion construct (Δ132) reveals an open conformation that relaxes steric constraints and facilitates repacking of the protein core. Sedimentation velocity studies show that the open conformation stabilizes the Δ132 construct as a hexamer with point group symmetry 32, similar to that of the active complex. In contrast, the UDX-inhibited enzyme forms a lower-symmetry, horseshoe-shaped hexameric complex. We show that the Δ132 and UDX-inhibited structures have similar hexamer-building interfaces, suggesting that the hinge-bending motion represents a path for the allosteric transition between the different hexameric states. On the basis of (i) main chain flexibility and (ii) a model of the conformational change, we propose that hinge bending can occur as a concerted motion between adjacent subunits in the high-symmetry hexamer. We combine these results in a structurally detailed model for allosteric feedback inhibition and substrate--product exchange during the catalytic cycle.     

Stereochemical analysis of the secondary structure of polypeptide chains with the aid of Courtauld three-dimensional models. II. Hydrogen and hydrophobic bonds.

Sterically permissible hydrogen bonds between side chains and the backbone, which fix a small number of angles of internal rotation have been identified. The hydrogen bonds on the ends and within the secondary structures, as well as the most important bonds in the alpha helix and in irregular structures, have been considered.