Comparative thermodynamics for monomer and dimer sequence-dependent binding of a heterocyclic dication in the DNA minor groove

Phenylamidine cationic groups linked by a furan ring (furamidine) and related symmetric diamidine compounds bind as monomers in the minor groove of AT sequences of DNA. DB293, an unsymmetric derivative with one of the phenyl rings of furamidine replaced with a benzimidazole, can bind to AT sequences as a monomer but binds more strongly to GC-containing minor-groove DNA sites as a stacked dimer. The dimer-binding mode has high affinity, is highly cooperative and sequence selective. In order to develop a better understanding of the correlation between structural and thermodynamic aspects of DNA molecular recognition, DB293 was used as a model to compare the binding of minor-groove agents with AT and mixed sequence DNA sites. Isothermal titration calorimetry and surface plasmon resonance results clearly show that the binding of DB293 and other related compounds into the minor groove of AT sequences is largely entropy-driven while the binding of DB293 as a dimer into the minor groove of GC-containing sequences is largely enthalpy-driven. At 25 degrees C, for example, the AT binding has DeltaG degrees, DeltaH degrees and TDeltaS degrees values of -9.6, -3.6 and 6.0 kcal/mol while the values for dimer binding to a GC-containing site are -9.0, -10.9 and -1.9 kcal/mol (per mol of bound compound), respectively. These results show that the thermodynamic components for binding of compounds of this type to DNA are very dependent on the structure, solvation and sequence of the DNA binding site.     

Interpretation of 15N NMR relaxation data of globular proteins using hydrodynamic calculations with HYDRONMR

HYDRONMR is an implementation of state of the art hydrodynamic modeling to calculate the spectral density functions for NH or C-H vectors in a rigid protein structure starting from an atomic level representation. Thus HYDRONMR can be used to predict NMR relaxation times from a rigid model and to compare them with the experimental results. HYDRONMR contains a single adjustable parameter, the atomic element radius. A protocol to determine the value that gives the best agreement between calculated and experimental T₁/T₂values is described. For most proteins, the value of the atomic element radius ranges between 2.8 Å and 3.8 Å with a distribution centered at 3.3 Å. Deviations from the usual range towards larger values are associated to aggregation in several proteins. Deviations to lower values may be related to large-scale motions or inappropriate model structures.If the average structure is correct, deviations between experimental T₁/T₂values and those calculated with HYDRONMR can be used to distinguish residues affected by anisotropic motion from those that are involved in chemical exchange.     

Protein molecular dynamics with electrostatic force entirely determined by a single Poisson-Boltzmann calculation

The electrostatic force including the intramolecular Coulombic interactions and the electrostatic contribution of solvation effect were entirely calculated by using the finite difference Poisson-Boltzmann method (FDPB), which was incorporated into the GROMOS96 force field to complete a new finite difference stochastic dynamics procedure (FDSD). Simulations were performed on an insulin dimer. Different relative dielectric constants were successively assigned to the protein interior; a value of 17 was selected as optimal for our system. The simulation data were analyzed and compared with those obtained from 500-ps molecular dynamics (MD) simulation with explicit water and a 500-ps conventional stochastic dynamics (SD) simulation without the mean solvent force. The results indicate that the FDSD method with GROMOS96 force field is suitable to study the dynamics and structure of proteins in solution if used with the optimal protein dielectric constant.     

Propionic acid side chain hydrogen bonding in the malaria pigment beta-hematin

Malaria pigment, or beta-hematin, the insoluble heme detoxification product resulting from the intraerythrocitic digestion of hemoglobin by young malaria trophozoites has been structurally characterized by X-ray powder diffraction and shown to contain chains of propionic acid linked dimers. Although there is considerable spectroscopic evidence for a monodentate propionate-iron interaction in this crystalline material, the spectroscopic characterization of the propionic acid dimer is limited. Herein we demonstrate the presence of the propionic acid dimer unit by H/D isotope substitution in carboxylic acid dimer. In the Raman spectrum of the deuterium substituted compound there is a circa 12 cm(-1) shift, H: 1629 cm(-1) vs. D: 1617 cm(-1) in the symmetric ring breathing mode for the propionic acid dimer. On the other hand, the IR active asymmetric stretch has a very small shift, <3 cm(-1), upon deuteration. These, and other vibrational data, are consistent with the presence of a planar carboxylic acid dimer in the structure of beta-hematin.     

Dimer structure of magainin 2 bound to phospholipid vesicles

Magainin 2 from African clawed frog Xenopus laevis is an antimicrobial peptide with broad spectra and action mechanisms considered to permeabilize bacterial membranes. CD, vibration, and solid-state NMR spectroscopies indicate the peptide adopts an alpha-helical conformation on binding to phospholipid bilayers, and its micelle-bound conformation, being monomeric and alpha-helical, is well detailed. We showed, however, that the peptide dimerizes on binding to phospholipid bilayers. This difference in the conformation and aggregation state between micelle- and bilayer-bound states prompted us to analyze the conformation of an equipotent analog of magainin 2 (F5Y,F16W magainin 2) bound to phosphatidylcholine vesicles using transferred nuclear Overhauser enhancement (TRNOE) spectroscopy. While observed medium-range TRNOE cross peaks were characteristic of alpha-helix, many long-range cross peaks were not compatible with the peptide's monomeric state. Simulated annealing calculations generated dimer structures indicating (1) two peptide molecules have a largely helical conformation in antiparallel orientation forming a short coiled-coil structure, (2) residues 4-20 are well converged and residues 9-20 are in an alpha-helical conformation, and (3) the interface of the two peptide molecules is formed by well-defined side chains of hydrophobic residues. Finally, determined structures are compatible with numerous investigations examining magainin-phospholipid interactions.     

A partially unfolded state of equine beta-lactoglobulin at pH 8.7

The urea-induced unfolding transition of equine -lactoglobulin was studied at pH 8.7 using circular dichroism (CD), ultracentrifugation, and gel filtration chromatography. The unfolding transition curves showed that at least one intermediate accumulates at moderate concentrations of urea. Furthermore, analytical ultracentrifugation experiments indicated that the intermediate forms a dimer. Thus, the urea-induced unfolding transition was measured by CD at various protein concentrations and was analyzed by a model assuming the four conformational states (the native, intermediate, dimeric intermediate, and unfolded states). The characteristics of the intermediate are markedly different from those of the intermediate previously observed at pH 4.0 or 1.5. The intermediate at pH 8.7 does not show the intense far-ultraviolet CD suggestive of the nonnative -helix.     

Dimer formation by a "monomeric" protein

Dimeric proteins can arise by the swapping of structural domains between monomers. The prevalence of this occurrence is unknown. Ribonuclease A (RNase A) is assumed to be a monomer near physiological conditions. Here, this hypothesis is tested and found to be imprecise. The two histidine residues (His12 and His119) in the active site of RNase A arise from two domains (S-peptide and S-protein) of the protein. The H12A and H119A variants have 10(5)-fold less ribonucleolytic activity than does the wild-type enzyme. Incubating a 1:1 mixture of the H12A and H119A variants at pH 6.5 and 65 degrees C results in a 10(3)-fold increase in ribonucleolytic activity. A large quantity of active dimer can be produced by lyophilizing a 1:1 mixture of the H12A and H119A variants from acetic acid. At pH 6.5 and 65 degrees C, the ribonucleolytic activity of this dimer converges to that of the dimer formed by simply incubating the monomers, as expected for a monomer-dimer equilibrium. The equilibrium dissociation constant for the dimer is near 2 mM at both 65 and 37 degrees C. This value of Kd is only 20-fold greater than the concentration of RNase A in the cow pancreas, suggesting that RNase A dimers exist in vivo. The intrinsic ability of RNase A to form dimers under physiological conditions is consistent with a detailed model for the evolution of homodimeric proteins. Dimers of "monomeric" proteins could be more prevalent than is usually appreciated.     

Ab Initio and DFT Study of the Structures,Vibrations and Energetics of the Urea Dimer and Trimer

Stable structures and electronic properties of small urea clusters are investigated with ab initio calculations. We optimized the cluster geometries and calculated the vibrational frequencies with Hartree–Fock (HF), second-order Møller-Plesset perturbation theory (MP2), and Density Functional Theory (DFT) methods using different basis sets. The most stable dimer was found to consist of two nonplanar urea molecules which are connected by two N–-H...O bonds in a common plane, and the most stable trimer has a flat structure of complex and planar C2 form for each urea molecule, like in the crystal. The interaction energies were corrected for the basis set superposition error (BSSE) using the full Boys*ndash;Bernardi counterpoise correction scheme. The stability of different dimer and trimer structures, the features of formation of H-bonds and presented here are compared to the available experimental data.     

Difference in the structures of alanine tri‐ and tetra‐peptides with antiparallel β‐sheet assessed by X‐ray diffraction,solid‐state NMR and chemical shift calculations by GIPAW

Alanine oligomers provide a key structure for silk fibers from spider and wild silkworms.We report on structural analysis of l ‐alanyl‐l ‐alanyl‐l ‐alanyl‐l ‐alanine (Ala)₄ with anti‐parallel (AP) β‐structures using X‐ray and solid‐state NMR. All of the Ala residues in the (Ala)₄ are in equivalent positions, whereas for alanine trimer (Ala)₃ there are two alternative locations in a unit cell as reported previously (Fawcett and Camerman, Acta Cryst., 1975, 31, 658–665). (Ala)₄ with AP β‐structure is more stable than AP‐(Ala)₃ due to formation of the stronger hydrogen bonds. The intermolecular structure of (Ala)₄ is also different from polyalanine fiber structure, indicating that the interchain arrangement of AP β‐structure changes with increasing alanine sequencelength. Furthermore the precise ¹H positions, which are usually inaccesible by X‐ray diffraction method, are determined by high resolution ¹H solid state NMR combined with the chemical shift calculations by the gauge‐including projector augmented wave method. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 13–20, 2014.