Solution H NMR characterization of substrate-free C. diphtheriae heme oxygenase: Pertinence for determining magnetic axes in paramagnetic substrate complexes

Proton 2D NMR was used to confirm in solution a highly conserved portion of the molecular structure upon substrate loss for the heme oxygenase from the pathogenic bacterium Corynebacterium diphtheriae, HmuO. The chemical shifts for the conserved portion of the structure are assessed as references for the dipolar shifts needed to determine the orientation of the paramagnetic susceptibility tensor, χ, in paramagnetic substrate complexes of HmuO. It is shown that the chemical shifts for the structurally conserved portion of substrate-free HmuO serve as excellent references for residues with only small to moderate sized dipolar shifts in the cyanide-inhibited substrate complex of HmuO, yielding an orientation of χ that is essentially the same as conventionally obtained from large dipolar shifts based on empirical estimates of the diamagnetic reference. The implications of these diamagnetic chemical shifts for characterizing the hydrogen bonding in the physiologically relevant, resting-state, high-spin aquo complex are discussed. The pattern of labile proton exchange in the distal H-bond network of substrate-free HmuO allowed comparison of changes in dynamic stability of tertiary contacts in the substrate-free and substrate-bound HmuO and with the same complexes of human heme oxygenase.     

Quantum-chemical analysis of C-H...O and C-H...N interactions in RNA base pairs--H-bond versus anti-H-bond pattern

Geometries and interaction energies of unusual UU and AA base pairs with one standard hydrogen bond (H-bond) and additional C-H...O or C-H...N contacts have been determined by quantum-chemical methods taking into account electron correlation. Whereas the C-H bond length in the UU C-H...O contact increases upon complex formation (H-bond pattern), the C-H bond of the AA C-H....N interaction is shortened (anti-H-bond pattern). The same properties are found for model complexes between U or A and formaldehyde that have intermolecular C-H...acceptor contacts but no standard H-bonds. Both the C-H...acceptor H-bond and anti-H-bond interactions are attractive. A possible influence of the donor CH group charge distribution on the interaction pattern is discussed.     

Quantum-Chemical Analysis of C-H…O and C-H…N Interactions in RNA Base Pairs—H-Bond Versus Anti-H-Bond Pattern

Abstract

Geometries and interaction energies of unusual UU and AA base pairs with one standard hydrogen bond (H-bond) and additional C-H…O or C-H…N contacts have been determined by quantum-chemical methods taking into account electron correlation. Whereas the C-H bond length in the UU C-H…O contact increases upon complex formation (H-bond pattern), the C-H bond of the AA C-H….N interaction is shortened (anti-H-bond pattern). The same properties are found for model complexes between U or A and formaldehyde that have intermolecular C-H…acceptor contacts but no standard H-bonds. Both the C-H…acceptor H-bond and anti-H-bond interactions are attractive. A possible influence of the donor CH group charge distribution on the interaction pattern is discussed.     

Solution 1H NMR investigation of the active site molecular and electronic structures of substrate-bound,cyanide-inhibited HmuO,a bacterial heme oxygenase from Corynebacterium diphtheriae

The molecular structure and dynamic properties of the active site environment of HmuO, a heme oxygenase (HO) from the pathogenic bacterium Corynebacterium diphtheriae, have been investigated by (1)H NMR spectroscopy using the human HO (hHO) complex as a homology model. It is demonstrated that not only the spatial contacts among residues and between residues and heme, but the magnetic axes that can be related to the direction and magnitude of the steric tilt of the FeCN unit are strongly conserved in the two HO complexes. The results indicate that very similar contributions of steric blockage of several meso positions and steric tilt of the attacking ligand are operative. A distal H-bond network that involves numerous very strong H-bonds and immobilized water molecules is identified in HmuO that is analogous to that previously identified in hHO (Li, Y., Syvitski, R. T., Auclair, K., Wilks, A., Ortiz de Montellano, P. R., and La Mar, G. N. (2002) J. Biol. Chem. 277, 33018-33031). The NMR results are completely consistent with the very recent crystal structure of the HmuO.substrate complex. The H-bond network/ordered water molecules are proposed to orient the distal water molecule near the catalytically key Asp(136) (Asp(140) in hHO) that stabilizes the hydroperoxy intermediate. The dynamic stability of this H-bond network in HmuO is significantly greater than in hHO and may account for the slower catalytic rate in bacterial HO compared with mammalian HO.     

1H NMR characterization of the solution active site structure of substrate-bound, cyanide-inhibited heme oxygenase from Neisseria meningitidis: comparison to crystal structures

Heme oxygenase, HO, from the pathogenic bacterium Neisseria meningitidis catabolizes heme for the iron necessary for infection. The enzyme, labeled HemO, exhibits less sequence homology to mammalian HO than another studied HO from Corynebacterium diphtheriae. Solution 1H NMR has been utilized to define the active site molecular and electronic structure of the cyanide-inhibited, substrate-bound complex for comparison with those provided by several crystal structures. Extensive assignments by solely 1H NMR 2D methods reveal a structure that is very strongly conserved with respect to the crystal structure, although 1H/2H exchange indicates dynamically much more stable distal and proximal helices than those for other HOs. Several residues found with alternate orientations in crystal structures of water- and NO-ligated complexes were shown to occupy positions found solely in the NO complex, confirming that there are structural accommodations in response to ligating the substrate complex with a diatomic, H-bond acceptor ligand. The observed dipolar shifts allow the determination of the magnetic axes that show that the Fe-CN unit is tilted approximately 10 degrees toward the alpha-meso position, thereby facilitating the alpha-stereoselectivity of the enzyme. Numerous labile protons with larger than usual low-field bias are identified and, in common with the other HO complexes, shown to participate in an extended, distal side H-bond network. This H-bond network orders several water molecules, most, but not all, of which have been detected crystallographically. A series of three C-terminal residues, His207-Arg208-His209, are not detected in crystal structures. However, 1H NMR finds two residues, His207 and likely Arg208 in contact with pyrrole D, which in crystal structures is exposed to solvent. The nature of the NOEs leads us to propose a H-bond between the proximally oriented His207 ring and the carboxylate of Asp27 and a salt-bridge between the terminus of Arg208 and the reoriented 7-propionyl carboxylate. While numerous ordered water molecules are found near both propionates in the crystal structure, we find much larger water NOEs to the 6- than 7-propionate, suggesting that water molecules near the 7-propionate have been expelled from the cavity by the insertion of Arg208 into the distal pocket. The conversion of the 7-propionate link from the N-terminal region (Lys16) to the C-terminal region (Arg208) in the ligated substrate complex both closes the heme cavity more tightly and may facilitate product exit, the rate-limiting step in the enzyme activity.     

Hydrogen-bonding interactions in adrenaline–water complexes: DFT and QTAIM studies of structures, properties, and topologies

ωB97XD/6-311++G(d,p) calculations were carried out to investigate the hydrogen-bonding interactions between adrenaline (Ad) and water. Six Ad-H(2)O complexes possessing various types of hydrogen bonds (H-bonds) were characterized in terms of their geometries, energies, vibrational frequencies, and electron-density topology. Natural bond orbital (NBO) and quantum theory of atoms in molecules (QTAIM) analyses were performed to elucidate the nature of the hydrogen-bonding interactions in these complexes. The intramolecular H-bond between the amino and carboxyl oxygen atom of Ad was retained in most of the complexes, and cooperativity between the intra- and intermolecular H-bonds was present in some of the complexes. H-bonds in which hydroxyls of Ad/water acted as proton donors were stronger than other H-bonds. Both hydrogen-bonding interactions and structural deformation play important roles in the relative stabilities of the complexes. The intramolecular H-bond was broken during the formation of the most stable complex, which indicates that Ad tends to break the intramolecular H-bond and form two new intermolecular H-bonds with the first water molecule.     

Ribosomal protein L1 recognizes the same specific structural motif in its target sites on the autoregulatory mRNA and 23S rRNA

The RNA-binding ability of ribosomal protein L1 is of profound interest since the protein has a dual function as a ribosomal protein binding rRNA and as a translational repressor binding its mRNA. Here, we report the crystal structure of ribosomal protein L1 in complex with a specific fragment of its mRNA and compare it with the structure of L1 in complex with a specific fragment of 23S rRNA determined earlier. In both complexes, a strongly conserved RNA structural motif is involved in L1 binding through a conserved network of RNA–protein H-bonds inaccessible to the solvent. These interactions should be responsible for specific recognition between the protein and RNA. A large number of additional non-conserved RNA–protein H-bonds stabilizes both complexes. The added contribution of these non-conserved H-bonds makes the ribosomal complex much more stable than the regulatory one.     

CH...O hydrogen bonds at protein-protein interfaces

For the first time, a statistical potential has been developed to quantitatively describe the CH.O hydrogen bonding interaction at the protein-protein interface. The calculated energies of the CH.O pair interaction show a favorable valley at approximately 3.3 A, exhibiting a feature typical of an H-bond and similar to the ab initio quantum calculation result (Scheiner, S., Kar, T., and Gu, Y. (2001) J. Biol. Chem. 276, 9832-9837). The potentials have been applied to a set of 469 protein-protein complexes to calculate the contribution of different types of interactions to each protein complex: the average energy contribution of a conventional H-bond is approximately 30%; that of a CH.O H-bond is 17%; and that of a hydrophobic interaction is 50%. In some protein-protein complexes, the contribution of the CH.O H-bond can reach as high as approximately 40-50%, indicating the importance of the CH.O H-bond at the protein interface. At the interfaces of these complexes, C(alpha)H.O H-bonds frequently occur between adjacent strands in both parallel and antiparallel orientations, having the obvious structural motif of bifurcated H-bonds. Our study suggests that the weak CH.O H-bond makes an important contribution to the association and stability of protein complexes and needs more attention in protein-protein interaction studies.     

Hydrogen bonding interactions in noradrenaline-DMSO complexes: DFT and QTAIM studies of structure,properties and topology

The hydrogen bonding interactions between noradrenaline (NA) and DMSO were studied with density functional theory (DFT) regarding their geometries, energies, vibrational frequencies, and topological features of the electron density. The quantum theory of atoms in molecules (QTAIM) and the natural bond orbital (NBO) analyses were employed to elucidate the hydrogen bonding interaction characteristics in noradrenaline-DMSO complexes. The H-bonds involving the hydroxyls hydrogen in NA and the O atom in DMSO are dominant intermolecular H-bonds and are stronger than other H-bonds involving the methyl hydrogen of DMSO as a H-donor. The weak H-bonds also include a π H-bond which involves the benzene ring as a H-donor or H-acceptor. QTAIM identified the weak H-bonds formed between the methyl hydrogen of DMSO and the N atom in NA in some complexes (AB5, AB6 and AB7), which cannot be further confirmed by NBO and other methods, so there are probably no interactions between hydrogen and nitrogen atoms among these complexes. A good linear relationship between logarithmic electron density (lnρ _b ) at the bond critical point (BCP) and structural parameter (δR _H···Y) was found. The formations of new H-bonds in some complexes are helpful to strengthen the original intramolecular H-bond, this is attributed to the cooperativity of H-bonds in complexes and can be learned from the structure results and the NBO and QTAIM analyses. Analysis of various physically meaningful contributions arising from the energy decomposition procedures show that the orbital interactions of H-bond is predominant during the formation of the complex, moreover, both the hydrogen bonding interaction and the structural deformation are responsible for the stability of the complexes.     

Non-classical hydrogen bonds in interleukins: the role of C-H...O interactions

The role of classical hydrogen bonds in the structural stability of biological macro-molecules is well understood. In the present study, we explore the influence of C-H...O interactions in relation to other environmental preferences in interleukins. Main chain-main chain interactions are predominant. Pro residues might stabilize helices and strands by C-H...O H-bonds in interleukins. Majority of the C-H...O interacting residues were solvent exposed. 62% of C-H...O interactions was long-range interactions. The results presented in this study might be useful for structural stability studies in interleukins.     

Dissecting water binding sites at protein–protein interfaces: a lesson from the atomic structures in the Protein Data Bank

We dissect the protein–protein interfaces into water preservation (WP), water hydration (WH) and water dehydration (WD) sites by comparing the water-mediated hydrogen bonds (H-bond) in the bound and unbound states of the interacting subunits. Upon subunit complexation, if a H-bond between an interface water and a protein polar group is retained, we assign it as WP site; if it is lost, we assign it as WD site and if a new H-bond is created, we assign it as WH site. We find that the density of WD sites is highest followed by WH and WP sites except in antigen and (or) antibody complexes, where the density of WH sites is highest followed by WD and WP sites. Furthermore, we find that WP sites are the most conserved followed by WD and WH sites in all class of complexes except in antigen and (or) antibody complexes, where WD sites are the most conserved followed by WH and WP sites. A significant number of WP and WH sites are involved in water bridges that stabilize the subunit interactions. At WH sites, the residues involved in water bridges are significantly better conserved than the other residues. However, no such difference is observed at WP sites. Interestingly, WD sites are generally replaced with direct H-bonds upon subunit complexation. Significantly, we observe many water-mediated H-bonds remain preserved in spite of large conformational changes upon subunit complexation. These findings have implications in predicting and engineering water binding sites at protein–protein interfaces.     

Repeated parallel evolution of minimal rRNAs revealed from detailed comparative analysis

The concept of a minimal ribosomal RNA-containing ribosome, a structure with a minimal set of elements capable of providing protein biosynthesis, is essential for understanding this fundamental cellular process. Nematodes and trypanosomes have minimal mitochondrial rRNAs and detailed reconstructions of their secondary structures indicate that certain conserved helices have been lost in these taxa. In contrast, several recent studies on acariform mites have argued that minimal rRNAs may evolve via shortening of secondary structure elements but not the loss of these elements as shown for trypanosomes and nematodes. Based on extensive structural analysis of chelicerate arthropods, we demonstrate that extremely short rRNAs of acariform mites share certain structural modifications with nematodes and trypanosomes: loss of helices of the GTPase region and divergence in the evolutionarily conserved connecting loop between helices H1648 and H1764 of the large subunit rRNA. These highly concerted parallel modifications indicate that minimal rRNAs were generated under the strong selection that favored or tolerated reductions of helices in particular locations while maintaining the functionality of the rRNA molecules throughout evolution. We also discuss potential evolution of minimal rRNAs and atypical transfer RNAs.     

Structure of mammalian cytochrome P450 2B4 complexed with 4-(4-chlorophenyl)imidazole at 1.9-A resolution: insight into the range of P450 conformations and the coordination of redox partner binding

A 1.9-A molecular structure of the microsomal cytochrome P450 2B4 with the specific inhibitor 4-(4-chlorophenyl)imidazole (CPI) in the active site was determined by x-ray crystallography. In contrast to the previous experimentally determined 2B4 structure, this complex adopted a closed conformation similar to that observed for the mammalian 2C enzymes. The differences between the open and closed structures of 2B4 were primarily limited to the lid domain of helices F through G, helices B' and C, the N terminus of helix I, and the beta(4) region. These large-scale conformational changes were generally due to the relocation of conserved structural elements toward each other with remarkably little remodeling at the secondary structure level. For example, the F' and G' helices were maintained with a sharp turn between them but are placed to form the exterior ceiling of the active site in the CPI complex. CPI was closely surrounded by residues from substrate recognition sites 1, 4, 5, and 6 to form a small, isolated hydrophobic cavity. The switch from open to closed conformation dramatically relocated helix C to a more proximal position. As a result, heme binding interactions were altered, and the putative NADPH-cytochrome P450 reductase binding site was reformed. This suggests a structural mechanism whereby ligand-induced conformational changes may coordinate catalytic activity. Comparison of the 2B4/CPI complex with the open 2B4 structure yields insights into the dynamics involved in substrate access, tight inhibitor binding, and coordination of substrate and redox partner binding.     

Solvent-induced differentiation of protein backbone hydrogen bonds in calmodulin

In apo and holoCaM, almost half of the hydrogen bonds (H-bonds) at the protein backbone expected from the corresponding NMR or X-ray structures were not detected by h3JNC' couplings. The paucity of the h3JNC' couplings was considered in terms of dynamic features of these structures. We examined a set of seven proteins and found that protein-backbone H-bonds form two groups according to the h3JNC' couplings measured in solution. H-bonds that have h3JNC' couplings above the threshold of 0.2 Hz show distance/angle correlation among the H-bond geometrical parameters, and appear to be supported by the backbone dynamics in solution. The other H-bonds have no such correlation; they populate the water-exposed and flexible regions of proteins, including many of the CaM helices. The observed differentiation in a dynamical behavior of backbone H-bonds in apo and holoCaM appears to be related to protein functions.     

Long-term molecular dynamics simulation of copper plastocyanin in water

A long molecular dynamics simulation (1.1 ns) of fully hydrated plastocyanin has been performed and analysed to relate protein dynamics to structural elements and functional properties. The solvated structure is described in detail by the analysis of H-bond network. During all the simulation, the crystal H-bond network is maintained in the beta-sheet regions, while several H-bonds are broken or formed on the external surface of the protein. To evaluate whether such changes could be due to conformational rearrangements or to solvent competition, we have examined the average number of H-bonds between protein atoms and water molecules, and the root mean square deviations from crystal structure as a function of protein residues. Protein mobility and flexibility have been examined by positional and dihedral angle rms fluctuations. Finally, cross-correlation maps have revealed the existence of correlated motions among residues connected by hydrogen bonds.     

New light on allostery: dynamic resonance Raman spectroscopy of hemoglobin kempsey

On the basis of static and time-resolved resonance Raman spectroscopy of HbA and of a mutant, HbK (Dalpha99N), a specific reaction coordinate is proposed for the allosteric transition in human hemoglobin. The heme is held between proximal (F) and distal (E) helices, whose orientation is responsive to forces generated by ligation and deligation. The E and F helices are in turn tethered via H-bonds to the A and H helices. These outer helices follow the E-F motion, thereby repositioning the N- and C-termini, which form the intersubunit salt bridges in the T quaternary structure. When the T state interface is weakened by Asp --> Asn substitution at a quaternary H-bond (HbK), the Fe-His bond is relaxed and becomes responsive to allosteric effectors. The same E-F motion is observed in HbK, but the A-H following motion is delayed, relative to HbA, as is the Asn H-bond formation.     

Change in entropy of the free polypeptide chain during formation of hydrogen bonds]

Formation probabilities of different hydrogen bonds between carbonyl oxygen and amide hydrogen were determined by Monte Carlo simulations using a computer model in the space of sterically allowable conformations of alanine and glycine oligopeptides, and the corresponding entropy losses for the peptide backbone, T delta S, were calculated. The model was studied at different criteria of steric interactions. Comparison with the data of other authors showed the values of T delta S to be mainly determined by overall extent and type of the state space and to be only slightly dependent on its energy profile. Both short-range and long-range steric interactions were shown to prevent hydrogen bonding, especially in alanine peptides. In the model studied, the initiation of alpha(R)-helices is associated with T delta S = 8-10 kT, and prior formation of a 3/10-turn or one three-center H-bond does not appreciably decrease this entropy barrier. Elongation of the alpha(R)-helix by one residue leads to T delta S = 3.0-3.7 kT, the helices begin to stabilize after at least three sequential H-bonds are formed. The difference in the probability of insertion of Ala and Gly into the helix is lower than it follows from comparison of their mobility. The results could be explained assuming that factors different from helical H-bonds take part in the stabilization of the helices. One may suppose upon modeling of folding that even three sequential H-bonds are unable to fix the structure of a flexible peptide loop, while the elongation of alpha(R)-helices in the supersecondary helix-loop-helix structure is favorable as long as the loop conformation remains nearly optimal.     

Truncation of subunit ND2 disrupts the threefold symmetry of the antiporter-like subunits in complex I from higher metazoans

Three of the conserved, membrane-bound subunits in NADH:ubiquinone oxidoreductase (complex I) are related to one another, and to Mrp sodium-proton antiporters. Recent structural analysis of two prokaryotic complexes I revealed that the three subunits each contain fourteen transmembrane helices that overlay in structural alignments: the translocation of three protons may be coordinated by a lateral helix connecting them together (Efremov, R.G., Baradaran, R. and Sazanov, L.A. (2010). The architecture of respiratory complex I. Nature 465, 441-447). Here, we show that in higher metazoans the threefold symmetry is broken by the loss of three helices from subunit ND2; possible implications for the mechanism of proton translocation are discussed.     

Quantitative comparison of the hydrogen bond network of A-state and native ubiquitin by hydrogen bond scalar couplings

The backbone hydrogen bond (H-bond) network of the partially folded A-state of ubiquitin (60% methanol, 40% water, pH 2) has been characterized quantitatively by (h3)J(NC)(') H-bond scalar couplings between the (15)N nuclei of amino acid H-bond donors and the (13)C carbonyl nuclei of the acceptors. Results on (h3)J(NC)(') couplings and the amide proton ((1)H(N)) chemical shifts for the A-state are compared quantitatively to the native state. The (h3)J(NC)(') correlations of the A-state show intact, nativelike H-bonds of the first beta-hairpin beta1/beta2 and the alpha-helix, albeit at lower strength, whereas the H-bonds in the C-terminal part change from a pure beta-structure to an all alpha-helical H(N)(i)-->O(i-4) connectivity pattern. A residue-specific analysis reveals that the conformations within the conserved secondary structure segments are much more homogeneous in the A-state than in the native state. Thus, the strong asymmetry of (h3)J(NC)(') couplings and (1)H(N) chemical shifts between the interior and exterior sides of the native state alpha-helix vanishes in the A-state. This indicates that the bend of this helix around the native state hydrophobic core is released in the homogeneous solvent environment of the A-state. Similarly, an irregularity in the behavior of H-bond I3-->L15 in hairpin beta1/beta2, which results from strong contacts to strand beta5 in the native state, is absent in the A-state. These findings rationalize the behavior of the (1)H(N) chemical shifts in both states and indicate that the A-state is in many aspects similar to the onset of thermal denaturation of the native state.     

Proposal of a new hydrogen-bonding form to maintain curdlan triple helix

Curdlan and other beta-1,3-D-glucans form right-handed triple helices, and it has been believed that the intermolecular H-bond is present at the center of the helix to maintain the structure. In this H-bond model, three secondary OH groups form an inequilateral hexagonal shape perpendicular to the helix axis. This hexagonal form seems to be characteristic for beta-1,3-D-glucans and is widely accepted. We carried out MOPAC and ab initio calculations for the curdlan helix, and we propose a new intermolecular H-bonding model. In our model, the H-bonds are formed between the O2-atoms on different x-y planes along the curdlan helix, hence the H-bonds are not perpendicular to the helix axis. The new H-bonds are connected along the helix, traversing three curdlan chains to make a left-handed helix. Therefore, the H-bonding array leads to a reverse helix of the main chain. According to our MOPAC calculation, this model is more stable than the previous one. We believe that the continuous H-bonding array is stabilized by cooperative phenomena in the polymeric system.