Structural dependencies of h3JNC' scalar coupling in protein H-bond chains

Study on the influence of oxidative stress on the fibrillization of fibrinogen

Human fibrinogen is an important coagulation factor as well as an independent predictor of coronary heart disease and stroke. Analysis of dysfibrinogens may provide useful information and help us to understand the molecular defects in fibrin polymerization. In the present study, we investigated the influence of oxidative stress of fibrinogen induced by H2O2 on the polymerization state of fibrin. UV absorbance spectroscopy, circular dichroism, ζ-potential, dynamic light scattering and steady shear viscosity were all employed to study the influence of oxidative stress on the molecular structure, the surface charges, and the size and shape of fibrinogen molecules. The fibrin morphology obtained was imaged and investigated using atomic force microscopy. The results demonstrated that the cross-linking, branching and height distribution of formed fibrin will be influenced by the oxidative stress of fibrinogen. This study presents new insights into the aggregation behaviour of fibrinogen and will be helpful to understand the formation mechanism of thrombosis under oxidative stress.

Explicit Hydrogen-Bond Potentials and Their Application to NMR Scalar Couplings in Proteins

Hydrogen bonds (H bonds) are fundamental for the stability, structure, and dynamics of chemically and biologically relevant systems. One of the direct means to detect H bonds in proteins is NMR spectroscopy. As H bonds are dynamic in nature, atomistic simulations offer a meaningful way to characterize and analyze properties of hydrogen bonds, provided a sufficiently accurate interaction potential is available. Here, we use explicit H-bond potentials to investigate scalar coupling constants (h3)JNC' and characterize the conformational ensemble for increasingly accurate intermolecular potentials. By considering a range of proteins with different overall topology a general procedure to improve the hydrogen-bonding potential ("morphing potentials") based on experimental information is derived. The robustness of this approach is established through explicit simulations in full solvation and comparison with experimental results. The H-bond potentials used here lead to more directional H bonds than conventional electrostatic representations employed in molecular mechanics potentials. It is found that the optimized potentials lead to H-bond geometries in remarkable agreement with previous ab initio and knowledge-based approaches to H bonds in model systems and in proteins. This suggests that, by combining theory, computation, and experimental data, H-bonding potentials can be improved and are potentially useful to better study coupling, energy transfer, and allosteric communication in proteins.

CHARMM36 all-atom additive protein force field: validation based on comparison to NMR data

Protein structure and dynamics can be characterized on the atomistic level with both nuclear magnetic resonance (NMR) experiments and molecular dynamics (MD) simulations. Here, we quantify the ability of the recently presented CHARMM36 (C36) force field (FF) to reproduce various NMR observables using MD simulations. The studied NMR properties include backbone scalar couplings across hydrogen bonds, residual dipolar couplings (RDCs) and relaxation order parameter, as well as scalar couplings, RDCs, and order parameters for side-chain amino- and methyl-containing groups. It is shown that the C36 FF leads to better correlation with experimental data compared to the CHARMM22/CMAP FF and suggest using C36 in protein simulations. Although both CHARMM FFs contains the same nonbond parameters, our results show how the changes in the internal parameters associated with the peptide backbone via CMAP and the χ1 and χ2 dihedral parameters leads to improved treatment of the analyzed nonbond interactions. This highlights the importance of proper treatment of the internal covalent components in modeling nonbond interactions with molecular mechanics FFs.

Segmental labeling to study multidomain proteins

This chapter contains a review of methodologies and recent applications of segmental labeling for NMR structural studies of proteins and protein complexes. Segmental labeling is used to specifically label a segment of protein structure with NMR active nuclei, thus reducing NMR spectral complexity and greatly facilitating structural NMR studies of large multi-domain proteins. It can also be used to introduce a synthetic fragment into a protein structure to study post-translationally modified proteins. Detailed protocols describing segmental labeling techniques are also included.

The effect of polarization on multiple hydrogen-bond formation in models of self-assembling materials

We report density functional theory calculations at the B3LYP/D95(d,p) level on several different cyclic H-bonding dimers, where the monomers of each are connected by a pair of N-H···O=C H-bonding interactions, and the H-bonding donors and acceptors on each monomer are separated by polarizable spacers. Depending on the structures, the individual H-bonds vary in strength (enthalpy) by over a factor of four, from 2.41 to 10.99 kcal/mol. We attribute most of the variation in interaction energies to differences in the extent of polarization due to each of the H-bonds, which can either combine constructively or destructively. The dipole-dipole interactions between the pair of H-bonds also contribute somewhat to the relative stabilities. The relevance of these results to the design of self-assembling materials is discussed.

A scheme for rapid prediction of cooperativity in hydrogen bond chains of formamides, acetamides, and N-methylformamides

A scheme is proposed in this article to predict the cooperativity in hydrogen bond chains of formamides, acetamides, and N-methylformamides. The parameters needed in the scheme are derived from fitting to the hydrogen bonding energies of MP2/6-31+G** with basis set superposition error (BSSE) correction of the hydrogen bond chains of formamides containing from two to eight monomeric units. The scheme is then used to calculate the individual hydrogen bonding energies in the chains of formamides containing 9 and 12 monomeric units, in the chains of acetamides containing from two to seven monomeric units, in the chains of N-methylformamides containing from two to seven monomeric units. The calculation results show that the cooperativity predicted by the scheme proposed in this paper is in good agreement with those obtained from MP2/6-31+G** calculations by including the BSSE correction, demonstrating that the scheme proposed in this article is reasonable. Based on our scheme, a cooperativity effect of almost 240% of the dimer hydrogen bonding energy in long hydrogen bond formamide chains, a cooperativity effect of almost 190% of the dimer hydrogen bonding energy in long hydrogen bond acetamide chains, and a cooperativity effect of almost 210% of the dimer hydrogen bonding energy in long hydrogen bond N-methylformamide chains are predicted. The scheme is further applied to some heterogeneous chains containing formamide, acetamide, and N-methylformamide. The individual hydrogen bonding energies in these heterogeneous chains predicted by our scheme are also in good agreement with those obtained from Møller-Plesset calculations including BSSE correction.

Theoretical study of bifurcated hydrogen bonding effects on the 1J(N,H), 1hJ(N,H), 2hJ(N,N) couplings and 1H, 15N shieldings in model pyrroles

According to the density functional theory calculations, the X...H...N (X=N, O) intramolecular bifurcated (three-centered) hydrogen bond with one hydrogen donor and two hydrogen acceptors causes a significant decrease of the (1h)J(N,H) and (2h)J(N,N) coupling constants across the N-H...N hydrogen bond and an increase of the (1)J(N,H) coupling constant across the N-H covalent bond in the 2,5-disubstituted pyrroles. This occurs due to a weakening of the N-H...N hydrogen bridge resulting in a lengthening of the N...H distance and a decrease of the hydrogen bond angle at the bifurcated hydrogen bond formation. The gauge-independent atomic orbital calculations of the shielding constants suggest that a weakening of the N-H...N hydrogen bridge in case of the three-centered hydrogen bond yields a shielding of the bridge proton and deshielding of the acceptor nitrogen atom. The atoms-in-molecules analysis shows that an attenuation of the (1h)J(N,H) and (2h)J(N,N) couplings in the compounds with bifurcated hydrogen bond is connected with a decrease of the electron density rho(H...N) at the hydrogen bond critical point and Laplacian of this electron density nabla(2)rho(H...N). The natural bond orbital analysis suggests that the additional N-H...X interaction partly inhibits the charge transfer from the nitrogen lone pair to the sigma*(N-H) antibonding orbital across hydrogen bond weakening of the (1h)J(N,H) and (2h)J(N,N) trans-hydrogen bond couplings through Fermi-contact mechanism. An increase of the nitrogen s-character percentage of the N-H bond in consequence of the bifurcated hydrogen bonding leads to an increase of the (1)J(N,H) coupling constant across the N-H covalent bond and deshielding of the hydrogen donor nitrogen atom.

A theoretical structural analysis of the factors that affect (1)J(NH), (1h)J(NH) and (2h)J(NN) in N-H...N hydrogen-bonded complexes

Calculations of (1)J(NH), (1h)J(NH) and (2h)J(NN) spin-spin coupling constants of 27 complexes presenting N-H...N hydrogen bonds have allowed to analyze these through hydrogen-bond coupling as a function of the hybridization of both nitrogen atoms and the charge (+1, 0, - 1) of the complex. The main conclusions are that the hybridization of N atom of the hydrogen bond donor is much more important than that of the hydrogen bond acceptor. Positive and negative charges (cationic and anionic complexes) exert opposite effects while the effect of the transition states 'proton-in-the-middle' is considerable.

A review with comprehensive data on experimental indirect scalar NMR spin-spin coupling constants across hydrogen bonds

Scalar NMR spin-spin coupling constants across hydrogen bonds are fundamental in structural studies and as test grounds for theoretical calculations. Since they are scattered among many articles of different kinds, it seems useful to collect them in the most comprehensive way.

CLOPPA-IPPP analysis of cooperative effects in hydrogen-bonded molecular complexes. Application to intermolecular 2hJ(N,C) spin-spin coupling constants in linear (CNH)n complexes

The cooperative effects on NMR indirect nuclear coupling constants are analyzed by means of the IPPP-CLOPPA approach (where CLOPPA is the Contributions from Localized Orbitals within the Polarization Propagator Approach and IPPP is the Inner Projections of the Polarization Propagator). The decomposition of the J coupling allows one to classify these effects as those due to changes in the geometric structure and those that directly involve the transmission mechanisms. This latter contribution admits a further classification, taking into account its electronic origin. As an example, the cooperative effects on intermolecular 2hJ(N,C) couplings of the linear complexes (CNH)n (n = 2, 3, 4) are discussed.

Spin-spin coupling across intramolecular N-H(+)-N hydrogen bonds in models for proton sponges: an ab initio investigation

Ab initio calculations have been performed to obtain structures and coupling constants (1)J(N-H), (1h)J(H-N), and (2h)J(N-N) for models of proton sponges with symmetric and asymmetric N-H(+)-N intramolecular hydrogen bonds (IMHBs). For a given model, the asymmetric structure has a lower energy, a longer N-N distance, and a hydrogen bond which has a greater deviation from linearity. The computed values of (2h)J(N-N) for the models are significantly less than predicted values based on the distance dependence of (2h)J(N-N) for complexes with intermolecular N-H(+)-N hydrogen bonds. However, the reduced values of (2h)J(N-N) cannot be attributed solely to the distortion of the hydrogen bond in the models, but also reflect differences in s electron populations at the nitrogens in both the ground state and the excited states which couple to it through the Fermi-contact (FC) operator. Values of (2h)J(N-N) for IMHBs can be related quadratically to the N-N distances in the models, and demonstrate that there is no discrepancy between computed values of (2h)J(N-N) at the short N-N distances found in these systems and experimental data for proton sponges.

Structural dependencies of protein backbone 2JNC' couplings

Protein folding can introduce strain in peptide covalent geometry, including deviations from planarity that are difficult to detect, especially for a protein in solution. We have found dependencies in protein backbone (2)J(NC') couplings on the planarity and the relative orientation of the sequential peptide planes. These dependences were observed in experimental (2)J(NC') couplings from seven proteins, and also were supported by DFT calculations for a model tripeptide. Findings indicate that elevated (2)J(NC') couplings may serve as reporters of structural strain in the protein backbone imposed by protein folds. Such information, supplemented with the H-bond strengths derived from (h3)J(NC') couplings, provides useful insight into the overall energy profile of the protein backbone in solution.

Characterizing the Cooperativity in H-Bonded Amino Structures

Density functional theory calculations were used to examine the effect of H-bond cooperativity on the magnitude of the NMR chemical shifts and spin-spin coupling constants in a C4h-symmetric G-quartet and in structures consisting of six cyanamide monomers. These included two ring structures (a planar C6h-symmetric structure and a nonplanar S6-symmetric structure) and two linear chain structures (a fully optimized planar Cs-symmetric chain and a planar chain structure where all intra- and intermolecular parameters were constrained to be identical). The NMR parameters were computed for the G-quartet and cyanamide structures, as well as for shorter fragments derived from these assemblies without reoptimization. In the ring structures and the chain with identical monomers, the intra- and intermolecular geometries of the cyanamides were identical, thereby allowing the study of cooperative effects in the absence of geometry changes. The magnitude of the |1JNH| coupling, 1H and 15N chemical shifts of the H-bonding amino N-H group, and the |h2JNN| H-bond coupling increased, whereas the size of the |1JNH| coupling of the non-H-bonded amino N-H bonds of the first amino group in the chain, which are roughly perpendicular to the H-bonding network, decreased in magnitude when H-bonding monomers were progressively added to extending ring or chain structures. These effects are attributed to electron redistribution induced by the presence of the nearby H-bonding guanine or cyanamide molecules.

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.

Direct Calculation of trans-Hydrogen-Bond 13C-15N 3-Bond J-Couplings in Entire Polyalanine α-Helices. A Density Functional Theory Study

We present the first trans-H-bond 13C-15N 3-bond J couplings calculated from entire neutral and protonated alpha-helical polyalanines. The neutral helices considered are those of the capped peptides, acetyl(Ala)NNH2, where N = 8, 16, 17, and 18, while the protonated peptides are the uncapped (Ala)17 protonated at three different positions. The calculated J values correlate well with O...H distances and somewhat less well with N...O distances, particularly if the terminal H-bonds are eliminated from the correlation. The J values calculated using the entire helix are about 6% lower in magnitude than those recently reported for H-bonding chains whose geometries were extracted from the same helices. Aqueous solvation favors protonation of the alpha-helix on the terminal COOH. Experimental measurements of the trans-H-bond 13C-15N 3-bond J couplings in acidic solution should be interpreted with this context.

Density functional theory reveals an increase in the amino H1 chemical shift in guanine due to hydrogen bonding with water

Electronic structure calculations underestimate the chemical shift of the non H-bonded amino proton 1H22 in isolated G-quartet structures. The current work shows that this underestimation is due to the absence of a water environment in the calculations: coordination of at least two water molecules is required to obtain good agreement with experiment. The results indicate how improved agreement between calculated and experimental (solution-phase) NMR data can be obtained.

Characterization of the monovalent ion position and hydrogen-bond network in guanine quartets by DFT calculations of NMR parameters

Conformational stability of G-quartets found in telomeric DNA quadruplex structures requires the coordination of monovalent ions. Here, an extensive Hartree-Fock and density functional theory analysis of the energetically favored position of Li+, Na+, and K+ ions is presented. The calculations show that at quartet-quartet distances observed in DNA quadruplex structures (3.3 A), the Li+ and Na+ ions favor positions of 0.55 and 0.95 A outside the plane of the G-quartet, respectively. The larger K+ ion prefers a central position between successive G-quartets. The energy barrier separating the minima in the quartet-ion-quartet model are much smaller for the Li+ and Na+ ions compared with the K+ ion; this suggests that K+ ions will not move as freely through the central channel of the DNA quadruplex. Spin-spin coupling constants and isotropic chemical shifts in G-quartets extracted from crystal structures of K+- and Na+-coordinated DNA quadruplexes were calculated with B3LYP/6-311G(d). The results show that the sizes of the trans-hydrogen-bond couplings are influenced primarily by the hydrogen bond geometry and only slightly by the presence of the ion. The calculations show that the R(N2N7) distance of the N2-H2...N7 hydrogen bond is characterized by strong correlations to both the chemical shifts of the donor group atoms and the (h2)J(N2N7) couplings. In contrast, weaker correlations between the (h3)J(N1C6') couplings and single geometric factors related to the N1-H1...O6=C6 hydrogen bond are observed. As such, deriving geometric information on the hydrogen bond through the use of trans-hydrogen-bond couplings and chemical shifts is more complex for the N1-H1...O6=C6 hydrogen bond than for the N2-H2...N7 moiety. The computed trans-hydrogen-bond couplings are shown to correlate with the experimentally determined couplings. However, the experimental values do not show such strong geometric dependencies.

Three-dimensional structure of the water-insoluble protein crambin in dodecylphosphocholine micelles and its minimal solvent-exposed surface

We chose crambin, a hydrophobic and water-insoluble protein originally isolated from the seeds of the plant Crambe abyssinica, as a model for NMR investigations of membrane-associated proteins. We produced isotopically labeled crambin(P22,L25) (variant of crambin containing Pro22 and Leu25) as a cleavable fusion with staphylococcal nuclease and refolded the protein by an approach that has proved successful for the production of proteins with multiple disulfide bonds. We used NMR spectroscopy to determine the three-dimensional structure of the protein in two membrane-mimetic environments: in a mixed aqueous-organic solvent (75%/25%, acetone/water) and in DPC micelles. With the sample in the mixed solvent, it was possible to determine (>NH...OC<) hydrogen bonds directly by the detection of (h3)J(NC)' couplings. H-bonds determined in this manner were utilized in the refinement of the NMR-derived protein structures. With the protein in DPC (dodecylphosphocholine) micelles, we used manganous ion as an aqueous paramagnetic probe to determine the surface of crambin that is shielded by the detergent. With the exception of the aqueous solvent exposed loop containing residues 20 and 21, the protein surface was protected by DPC. This suggests that the protein may be similarly embedded in physiological membranes. The strategy described here for the expression and structure determination of crambin should be applicable to structural and functional studies of membrane active toxins and small membrane proteins.

Calcium-binding proteins afford calibration of dihedral-angle dependence of 3J(NC(gamma)) coupling constant in aspartate and asparagine residues

Calibration of the 3J(NC(gamma)) couplings across the N-C(alpha)-C(beta)-C(gamma) fragment of aspartate and asparagine residues is afforded by two interactions that produce fixed conformations of the side chains in solution. One is the binding of these side chains to calcium ions; the other is the H-bond interaction of these side chains with a backbone amide.

Calculation of trans-hydrogen-bond 13C-15N three-bond and other scalar J-couplings in cooperative peptide models. A density functional theory study

We report B3LYP DFT calculations on peptide models that consider the effects of cooperative interactions with proximate H-bonds and local geometry at the H-bonding site upon trans-H-bond (13)C-(15)N three-bond scalar J-couplings. The calculations predict that cooperative interactions with other H-bonds within a H-bonding chain can significantly increase the magnitude of these couplings. Such increases are due to a combination of the presence of the neighboring H-bonds and the slight increase in C=O distances expected for peptide H-bonds near the centers of H-bonding chains. The energies of H-bonds inferred from H-bonding distances, alone, could be significantly in error if the effects of neighboring H-bonds are ignored.

Peptaibol zervamicin IIb structure and dynamics refinement from transhydrogen bond J couplings

Zervamicin IIB (Zrv-IIB) is a channel-forming peptaibol antibiotic of fungal origin. The measured transhydrogen bond (3h)J(NC') couplings in methanol solution heaving average value of -0.41 Hz indicate that the stability of the Zrv-IIB helix in this milieu is comparable to the stability of helices in globular proteins. The N-terminus of the peptide forms an alpha-helix, whereas 3(10)-helical hydrogen bonds stabilize the C-terminus. However, two weak transhydrogen bond peaks are observed in a long-range HNCO spectrum for HN Aib(12). Energy calculations using the Empirical Conformation Energy Program for Peptides (ECEPP)/2 force field and the implicit solvent model show that the middle of the peptide helix accommodates a bifurcated hydrogen bond that is simultaneously formed between HN Aib(12) and CO Leu(8) and CO Aib(9). Several lowered (3h)J(NC') on a polar face of the helix correlate with the conformational exchange process observed earlier and imply dynamic distortions of a hydrogen bond pattern with the predominant population of a properly folded helical structure. The refined structure of Zrv-IIB on the basis of the observed hydrogen bond pattern has a small ( approximately 20 degrees ) angle of helix bending that is virtually identical to the angle of bending in dodecylphosphocholine (DPC) micelles, indicating the stability of a hinge region in different environments. NMR parameters ((1)HN chemical shifts and transpeptide bond (1)J(NC') couplings) sensitive to hydrogen bonding along with the solvent accessible surface area of carbonyl oxygens indicate a large polar patch on the convex side of the helix formed by three exposed backbone carbonyls of Aib(7), Aib(9), and Hyp(10) and polar side chains of Hyp(10), Gln(11), and Hyp(13). The unique structural features, high helix stability and the enhanced polar patch, set apart Zrv-IIB from other peptaibols (for example, alamethicin) and possibly underlie its biological and physiological properties.

Semisynthesis of a segmental isotopically labeled protein splicing precursor: NMR evidence for an unusual peptide bond at the N-extein-intein junction

Protein splicing is a posttranslational autocatalytic process in which an intervening sequence, termed an intein, is removed from a host protein, the extein. Although we have a reasonable picture of the basic chemical steps in protein splicing, our knowledge of how these are catalyzed and regulated is less well developed. In the current study, a combination of NMR spectroscopy and segmental isotopic labeling has been used to study the structure of an active protein splicing precursor, corresponding to an N-extein fusion of the Mxe GyrA intein. The (1)J(NC') coupling constant for the (-1) scissile peptide bond at the N-extein-intein junction was found to be approximately 12 Hz, which indicates that this amide is highly polarized, perhaps because of nonplanarity. Additional mutagenesis and NMR studies indicate that conserved box B histidine residue is essential for catalysis of the first step of splicing and for maintaining the (-1) scissile bond in its unusual conformation. Overall, these studies support the "ground-state destabilization" model as part of the mechanism of catalysis.

Investigation of the NMR Spin−Spin Coupling Constants across the Hydrogen Bonds in Ubiquitin: The Nature of the Hydrogen Bond as Reflected by the Coupling Mechanism

The indirect scalar NMR spin-spin coupling constants across the H-bonds of the protein ubiquitin were calculated, including the Fermi contact, the diamagnetic spin-orbit, the paramagnetic spin-orbit, and the spin dipole term, employing coupled perturbed density functional theory in combination with the B3LYP functional and different basis sets: (9s,5p,1d/5s,1p)[6s,4p,1d/3s,1p] and (11s,7p,2d/5s,1p)[7s,6p,2d/4s,2p]. Four different models based on either the crystal or the aqueous solution structure of ubiquitin were used to describe H-bonding for selected residue pairs of ubiquitin. Calculated and measured 3hJ(NC') coupling constants differ depending on the model used, which is due to the fact that the geometry of ubiquitin is different in the solid state and in aqueous solution. Also, conformational averaging leads to a decrease of the magnitude of the measured 3hJ(NC') constants, which varies locally (larger for -sheets, smaller for -helix). Two different spin-spin coupling mechanisms were identified. While mechanism I transmits spin polarization via an electric field effect, mechanism II involves also electron delocalization from the lone pair of the carbonyl oxygen to the antibonding orbital of the N-H bond. Mechanism I is more important in the crystal structure of ubiquitin, while in aqueous solution, mechanism II plays a larger role. It is possible to set up simple relationships between the spin-spin coupling constants associated with the H bond in proteins and the geometrical features of these bonds. The importance of the 3hJ(NC') and 1J(N-H) constants as descriptors for the H-bond is emphasized.

H-bonding mediates polarization of peptide groups in folded proteins

The carbon-nitrogen J-couplings in the hydrogen bonding chains of proteins show that H-bonding mediates peptide-group polarization, which results in the general reduction of peptide-group polarity of folded proteins in solution. The net effect is to make large regions of protein secondary structure, especially beta-sheets, intrinsically more hydrophobic, contributing thereby to overall stability of the tertiary structure.