ACTH binding to phospholipid bilayers: A 1H-NMR study of the 5-14, 1-14, and 1-24 derivatives

The interaction of the 5-14, 1-14, and 1-24 fragments of ACTH with sonicated phospholipid bilayers containing egg yolk phosphatidylcholine (EPC) either pure or mixed with 10 mole % phosphatidic acid (EPA), was investigated by proton nuclear magnetic resonance (¹H-nmr). The effects observed with zwitterionic EPC vesicles were small, indicating a low binding of the ACTH derivatives. The N-terminal aromatic resonances of the ACTH peptides were markedly broadened in the presence of negatively charged vesicles (EPC/EPA 9:1 M/M), while those of the C-terminal end were barely affected, showing that ACTH interacts with its N-terminal fragment. The choline resonance of the EPC molecules of the outer monolayer was shifted and broadened upon ACTH binding to the lipid vesicles, while that of the inner layer was not affected, suggesting that the peptide molecules interact only with the external leaflet of the lipid bilayer. The C²H and C⁴H resonances of the histidine-6 side chain were both shifted downfield upon peptide binding to the negatively charged lipid interface. In the case of the 1–24 derivative, these resonances were also split into two signals reflecting two different species of membrane-bound ACTH 1–24. Analysis of the line width and chemical shift variations of the ACTH and lipid resonances observed upon peptide binding shows that the membrane-binding potency of the shorter 5–14⁺¹ fragment, which presents a +1 net charge, is roughly similar to that of the highly cationic 1–24⁺⁶ (net charge +6) derivative, implying that the 15–24⁺⁵ segment is not essential for membrane binding. The nmr measurements at a fixed lipid-to-peptide ratio in the presence of increasing amounts of spin-labeled lipids demonstrate that the N-terminal fragment of ACTH does not penetrate the hydrophobic core of the bilayer, and should lie parallel to the membrane surface. © 1997 John Wiley & Sons, Inc. Biopoly 42: 731–744, 1997     

A computational investigation of the interaction between As3+ and deoxynucleotides

A successful aptamer biosensor for As³⁺ in aqueous solution indicated that deoxynucleotides could be a potential ligand for arsenic absorption. However, the interaction mechanism between aptamer and As³⁺ is still not clear. In this paper, a density functional method was applied to investigate the interactions between As³⁺ and four kinds of deoxynucleotides in different charge states. According to the binding energy, dGMP was the best monomer with four different binding modes, and dGMP^3–As³⁺ with the O atom in the phosphate group was calculated to be the most stable one. Moreover, the interaction energy increased when the charge of deoxynucleotide was increased, suggesting that a negative charge was advantageous for deoxynucleotides to absorb As³⁺. The atoms in molecules (AIM), localised orbital locator (LOL) and Charge decomposition analyses (CDA) confirmed that the bond between N atoms and As³⁺ were covalent interactions and that between O atoms and As³⁺ were dominated by electrostatic interactions. As the charge grew, the bond strength of both interactions increased. At the same time, O atom in the vicinity of nucleobases began to participate in the interaction with As³⁺, which was originally bonded to N atoms.     

Vanadyl(IV) conalbumin. II. Mixed metal and anion binding studies

Electron paramagnetic resonance (epr) studies demonstrate that at low levels of conalbumin (CA) saturation with Fe³⁺ or VO²⁺, a ph-dependent preference of the metal exists for different protein binding-site configurations,A, B, and C. The vanadyl ion epr spectra of mixed VO²⁺, Fe³⁺-conalbumin in which Fe³⁺ is preferentially bound to the N- or C-terminal binding site are consistent with all three configurations being formed at both metal sites. At high pH the spectra suggest interaction between binding sites. In the absence of HCO₃^?, VO²⁺ is bound almost exclusively in B configuration; a full binding capacity of 2 VO²⁺ per CA is retained. Stoichiometric amounts of HCO₃^? convert the epr spectrum from B to an A, B, C type. Addition of oxalate to bicarbonate-free preparations converts the B spectrum to an A′, B, C′ type where the B resonances have lost intensity to the A′ and C′ resonances but have not changed position. The data suggest that configuration B is anion independent and that only one equivalent of binding sites at pH 9 responds to the presence of HCO₃^1? or oxalate by changing configuration but not metal binding capability. The form of the bound anion may be HCO₃^? rather than CO₃^2?. The formation rate of the colored ferric conalbumin complex by oxidizing Fe²⁺ to Fe³⁺ in limited HCO₃^? at pH 9 is also consistent with one equivalent of sites having different anion requirements than the remaining sites. Increased NaCl or NaClO₄ concentration or substitution of D₂O for water as solvent affect the environment of bound VO²⁺, but the mechanisms of action are unknown.     

Resolution enhanced homonuclear carbon decoupled triple resonance experiments for unambiguous RNA structural characterization

Large RNAs (>30 nucleotides) suffer from extensive resonance overlap that can seriously hamper unambiguous structural characterization. Here we present a set of 3D multinuclear NMR experiments with improved and optimized resolution and sensitivity for aiding with the assignment of RNA molecules. In all these experiments strong base and ribose carbon–carbon couplings are eliminated by homonuclear band-selective decoupling, leading to improved signal to noise and resolution of the C5, C6, and C1′ carbon resonances. This decoupling scheme is applied to base-type selective ¹³C-edited NOESY, ¹³C-edited TOCSY (HCCH, CCH), HCCNH, and ribose H1C1C2 experiments. The 3D implementation of the HCCNH experiment with both carbon and nitrogen evolution enables direct correlation of ¹³C and ¹⁵N resonances at different proton resonant frequencies. The advantages of the new experiments are demonstrated on a 36 nucleotides hairpin RNA from domain 5 (D5) of the group II intron Pylaiella littoralis using an abbreviated assignment strategy. These four experiments provided additional separation for regions of the RNA that have overlapped chemical shift resonances, and enabled the assignment of critical D5 bulge nucleotides that could not be assigned using current experimental schemes.Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1007/s10858-005-5093-6     

Investigation of the thermal unfolding of secondary and tertiary structure in E. coli tRNAfMet by high-resolution nmr

The high-resolution (300 MHz) proton nmr spectrum of E. coli tRNA^fMet has been examined in 0.17M NaCl, with and without Mg²⁺, and at various temperatures. In light of recent studies of other E. coli tRNA and fragments of tRNA^fMet, some low field (11–15 ppm) resonances previously assigned to secondary structure base pairs are reassigned to a tertiary structure A₁₄–S⁴U₈ base pair and a protected uridine residue in the anticodon loop. These two resonances and other low field resonances which are assigned to secondary structure base pairs are used to monitor the thermal unfolding of the molecule. In the absence of Mg²⁺ the tertiary structure base pair is present only to ～45°C, but in the presence of Mg²⁺ it remains until at least 70°C. Analysis of the temperature dependence of other low field resonances indicates that the melting of the dihydrouridine stem occurs more or less simultaneously with the loss of tertiary structure. The observation of the resonance from the A₁₄–S⁴U₈ base pair proves that tertiary structure is present in this molecule below 40°C, even in the absence of Mg²⁺.     

Mg2+ Dependence of the Structure and Thermodynamics of Wheat Germ and Lupin Seeds 5S rRNA

Abstract

The formation and stability of structural elements in two 5S rRNA molecules from wheat germ (WG) and lupin seeds (LS) as a function of Mg²⁺ concentration in solution was determined using the adiabatic differential scanning microcalorimetry (DSC). The experimentally determined thermodynamic parameters are compared with calculations using thermodynamic databases used for prediction of RNA structure. The 5S rRNA molecules which show minor differences in the nucleotide sequence display very different thermal unfolding profiles (DSC profiles). Numerical deconvolution of DSC profiles provided information about structural transformations that take place in both 5S rRNA molecules. A comparative analysis of DSC data and the theoretical thermodynamic models of the structure was used to establish a relationship between the constituting transitions found in the melting profiles and the unfolding of structural domains of the 5S rRNA and stability of its particular helical elements.

Increased concentration of Mg²⁺ ions induces additional internal interactions stabilising 5S rRNA structures found at low Na⁺ concentrations. Observed conformational transitions suggest a structural model in which the extension of helical region E dominates over the postulated tertiary interaction between hairpin loops. We propose that helix E is stabilised by a sequence of non-standard pairings extending this helix by the formation of tetra loop e and an almost total reduction of loop d between helices E and D. Two hairpin structures in both 5S rRNA molecules: the extended C-C' and the extended E-E'-E” hairpins appear as the most stable elements of the structure. The cooperativity of the unfolding of helixes in these 5S rRNA molecules changes already at 2 mM Mg²⁺.     

Pseudo-4D triple resonance experiments to resolve HN overlap in the backbone assignment of unfolded proteins

The solution NMR resonance assignment of the protein backbone is most commonly carried out using triple resonance experiments that involve ¹⁵N and ¹HN resonances. The assignment becomes problematic when there is resonance overlap of ¹⁵N–¹HN cross peaks. For such residues, one cannot unambiguously link the “left” side of the NH root to the “right” side, and the residues associated with such overlapping HN resonances remain often unassigned. Here we present a solution to this problem: a hybrid (4d,3d) reduced-dimensionality HN(CO)CA(CON)CA sequence. In this experiment, the Ca(i) resonance is modulated with the frequency of the Ca(i−1) resonance, which helps in resolving the ambiguity involved in connecting the Ca(i) and Ca(i−1) resonances for overlapping NH roots. The experiment has limited sensitivity, and is only suited for small or unfolded proteins. In a companion experiment, (4d,3d) reduced-dimensionality HNCO(N)CA, the Ca(i) resonance is modulated with the frequency of the CO(i−1) resonance, hence resolving the ambiguity existent in pairing up the Ca(i) and CO(i−1) resonances for overlapping NH roots.     

Structure and electronic absorption spectra of nematogenic alkoxycinnamic acids - a comparative study based on semiempirical and DFT methods

Structure of nematogenic p-n-Alkoxy cinnamic acids (nOCAC) with various alkyl chain carbon atoms (n = 2, 4, 6, 8) has been optimized using density functional B3LYP with 6-31+G (d) basis set using crystallographic geometry as input. Using the optimized geometry, electronic structure of the molecules has been evaluated using the semiempirical methods and DFT calculations. Molecular charge distribution and phase stability of these systems have been analyzed based on Mulliken and L?wdin population analysis. The electronic absorption spectra of nOCAC molecules have been simulated by employing DFT method, semiempirical CNDO/S and INDO/S parameterizations. Two types of calculations have been performed for model systems containing single and double molecules of nOCAC. UV-Visible spectra have been calculated for all single molecules. The UV stability of the molecules has been discussed in light of the electronic transition oscillator strength (f). The dimer complexes of higher homologues (n = 6, 8) have also been reported to enable the comparison between single and double molecules.     

Amino acid type-selective backbone 1H-15N-correlations for Arg and Lys

Four novel amino acid type-selective triple resonance experiments to identify the backbone amino proton and nitrogen resonances of Arg and Lys and of their sequential neighbors in (¹³C,¹⁵N)-labeled proteins are presented: the R(i+1)-HSQC and R(i,i+1)-HSQC select signals originating from Arg side chains, the K(i+1)-HSQC and K(i,i+1)-HSQC select signals originating from Lys side chains. The selection is based on exploiting the characteristic chemical shifts of a pair of carbon atoms in Arg and Lys side chains using selective 90° pulses. The new experiments are recorded as two-dimensional ¹H-¹⁵N-correlations and their performance is demonstrated with the application to a protein domain of 83 amino acids.     

Location and motion of free cholesterol molecules in high density lipoprotein

Human high density lipoprotein (HDL₃) was reconstituted with the free cholesterol molecules replaced with 4-[¹³C]-cholesterol. 90 MHz [¹³C]-NMR spectra were obtained and two cholesterol resonances at chemical shifts of 41.73 and 42.20 ppm could be resolved. The former signal arises from the C-4 atom of cholesterol molecules associated with phospholipids and located in the surface of the HDL₃ particle while the latter resonance is due to cholesterol molecules associated with cholesterol ester and triglyceride molecules in the core. HDL₃ reconstituted without any cholesterol ester or triglyceride gave a single resonance at 41.73 ppm indicating that all the free cholesterol molecules are in the surface. 60% of the free cholesterol molecules present in normal HDL₃ are in the phospholipid monolayer around the surface where they undergo relatively restricted motion compared to the remaining 40% situated in the liquid core. The free cholesterol molecules can equilibrate between the two pools in the timescale 10ms–700s.     

Electronic and molecular structure of M-DNA fragments

To study M-DNA molecular structure (such DNA with transition metal ions placed between the nucleic bases is able to conduct the electric current) and its conductivity mechanisms, we carried out ab initio quantum-mechanical calculations of electronic and spatial structures, thermodynamic characteristics of adenine-thymine (АТ) and guanine-cytosine (GC) base pair complexes with Zn²⁺ and Ni²⁺. To take into account the influence of the alkaline environment, calculations for these complexes were also carried out with hydroxyl and two water molecules. Computations were performed at MP2 level of theory using 6–31+G* basis set. Analogous calculations were carried out for (AC)(TG) stacking dimer of nucleic acid base pairs with two Zn²⁺. The calculation of the interaction energy in complexes has shown the preference of locating the metal ion (instead of the imino proton) between bases in M-DNA. The electronic transition energy calculation has revealed the reduction of the first singlet transition energy in АТ and GC complexes with Ni²⁺ from 4.5 eV to 0.4 - 0.6 eV. Ni²⁺ orbitals take part in the formation of HOMO and LUMO on the complexes investigated. It was shown that charges of metal ions incorporated into complexes with nucleic bases and in dimer decrease significantly.     

Theoretical studies on interaction of anticancer drugs (dacarbazine,procarbazine and triethylenemelamine) with normal (AT and GC) and mismatch (GG,CC, AA and TT) base pairs

This study is an attempt to gain a better understanding of the physicochemical interaction between novel anticancer drugs and DNA bases. We have employed quantum chemical tools to explore the interaction of a few anticancer drugs [namely procarbazine (PR), dacarbazine (DC) and triethylenemelamine (TR)] with isolated normal (GC and AT) and mismatch (AA, CC, GG and TT) base pairs. The molecular geometries, electronic structural stability, vibrational energies, chemical reactivity and other electronic properties were studied using MP2/6-311+G^**, B3LYP/6-311+G^** and M05-2X/6-311+G^** methods. The optimised geometries of the usual and mismatch base pairs are almost planar whereas the geometries of drug-interacting complexes deviate from planarity. The presence of steric hindrance and π-bond overlaps between C–C bonds in the complexes has distorted the planarity of the four- and five-member rings in the base pairs. Among the three drugs chosen, DC and PR bond well with normal and mismatch base pairs with large interaction energy. The electron density (ED) difference maps of the most stable GG–DC, GG–PR and GG–TR drug-interacting complexes show the information about sharing of ED and gain or loss of ED within the interacting molecules. The stabilisation energy of the charge transfer interaction between the relevant donor–acceptor orbital of GG–DC and GC–DC complexes has been found to be around 16 kcal/mol and GG–PR and GC–PR complexes has been found to be around 12 kcal/mol. But, for the GG–TR and GC–TR complexes, the stabilisation energy is found to be less than 6 kcal/mol. Moreover, the topological analysis of hydrogen bond network of DC and PR drug-interacting complexes have high electron and Laplacian density with structural stability at the bond critical points (BCPs), while compared TR drug-interacting complexes by atoms in molecules and natural bond orbital analysis. Finally, we may conclude that the drugs DC and PR are highly efficient drugs to target normal and mismatch base pair for control and inhibition of DNA replication.     

Is the decrease of the total electron energy density a covalence indicator in hydrogen and halogen bonds?

In this work, halogen bonding (XB) and hydrogen bonding (HB) complexes were studied with the aim of analyzing the variation of the total electronic energy density H(r _b ) with the interaction strengthening. The calculations were performed at the MP2/6?311++G(2d,2p) level of approximation. To explain the nature of such interactions, the atoms in molecules theory (AIM) in conjunction with reduced variational space self-consistent field (RVS) energy decomposition analysis were carried out. Based on the local virial theorem, an equation to decompose the total electronic energy density H(r _b ) in two energy densities, (?G(r _b )) and 1/4?²ρ(r _b ), was derived. These energy densities were linked with the RVS interaction energy components. Through the connection between both decomposition schemes, it was possible to conclude that the decrease in H(r _b ) with the interaction strengthening observed in the HB as well as the XB complexes, is mainly due to the increase in the attractive electrostatic part of the interaction energy and in lesser extent to the increase in its covalent character, as is commonly considered.     

X-ray and model-building studies on the specificity of the active site of proteinase K

Proteinase K, the extracellular serine endopeptidase (E.C. 3.4.21.14) from the fungus Tritirachium album limber, is homologous to the bacterial subtilisin proteases. The binding geometry of the synthetic inhibitor carbobenzoxy-Ala-Phechloromethyl Ketone to the active site of proteinase K was the first determined from a Fourier synthesis based on synchrotron X-ray diffraction data between 1.8 Å and 5.0 Å resolution. The protein inhibitor complexes was refined by restrained least-squares minimization with the data between 10.0 and 1.8 Å. The final R factor was 19.1% and the model contained 2,018 protein atoms, 28 inhibitors atoms, 125 water molecules, and two Ca²⁺ ions. The peptides portion of the inhibitor is bound to the active center of proteinase K by means of a three-stranded antiparallel pleated sheet, with the side chain of the phenylalanine located in the P1 site. Model building studies, with lysine replacing phenylalanine in the inhibitor, explain the relatively unspecific catalytic activity of the enzyme.     

A B3LYP and MP2(full) theoretical investigation into explosive sensitivity upon the formation of the molecule-cation interaction between the nitro group of 3,4-dinitropyrazole and H+, Li+, Na+, Be2+ or Mg2+

The explosive sensitivity upon the formation of molecule-cation interaction between the nitro group of 3,4-dinitropyrazole (DNP) and H⁺, Li⁺, Na⁺, Be²⁺ or Mg²⁺ has been investigated using the B3LYP and MP2(full) methods with the 6-311++G** and 6-311++G(2df,2p) basis sets. The bond dissociation energy (BDE) of the C3–N7 trigger bond has also been discussed for the DNP monomer and the corresponding complex. The interaction between the oxygen atom of nitro group and H⁺ in DNP…H⁺ is partly covalent in nature. The molecule-cation interaction and bond dissociation energy of the C3–N7 trigger bond follow the order of DNP…Be²⁺ > DNP…Mg²⁺ > DNP…Li⁺ > DNP…Na⁺. Except for DNP…H⁺, the increment of the trigger bond dissociation energy in comparison with the DNP monomer correlates well with the molecule-cation interaction energy, natural charge of the nitro group, electron density ρ _BCP(C3–N7), delocalization energy E ⁽²⁾ and NBO charge transfer. The analyses of atoms in molecules (AIM), natural bond orbital (NBO) and electron density shifts have shown that the electron density of the nitro group shifts toward the C3–N7 trigger bond upon the formation of the molecule-cation interaction. Thus, the trigger bond is strengthened and the sensitivity of DNP is reduced.     

Calmodulin structure and function: Implication of arginine in the compaction related to ligand binding

Calmodulin (CAM) is a modulatory protein that regulates cellular activity by binding to a large number of proteins. Key elements in the Ca²⁺-dependent mechanism of interaction between CAM and the proteins it activates are the selectivity for Ca²⁺ ions and the requirement for Ca²⁺-dependent conformational changes. We report on results from a series of molecular dynamics simulations that identified discrete steps in the mechanism of structural rearrangement of CAM. The findings implicate the side chains of arginine residues in the bending of the central alpha helix. Structural and energetic considerations point to a dynamic hydrogen bonding pattern around the arginine residues as a ratcheting-type mechanism, causing the kinking of the central helix in consecutive steps stabilized by each new pattern of hydrogen bonds. Initial model building studies to locate potential binding sites of ligands such as trifluoperazine (TFP) indicate that the compaction of CAM results in several structural changes, that explain the selective binding of molecules such as TFP in the N-terminal domain. The present studies identify specific residues involved in the process of compaction and point to specific CAM residues involved in the binding of the ligand. These insights lead directly to propositions for experimental engineering of the molecular structure of CAM in order to probe the hypotheses and their consequences for the function of this important protein.     

Stereochemistry of the Brivaracetam Diastereoisomers

The stereochemistry of all four stereoisomers of brivaracetam was determined using vibrational circular dichroism (VCD) spectroscopy. By comparing experimentally obtained VCD spectra and computationally simulated ones, the absolute configurations can be confidently assigned without prior knowledge of their relative stereochemistry. Neither the corrected mean absolute errors analysis of the nuclear magnetic resonance (NMR) data, nor the matching of experimental and calculated infrared spectra allowed the diastereoisomers to be distinguished. VCD spectroscopy itself suffices to establish the absolute configurations of all diastereoisomers. The relative stereochemistry could also be statistically confirmed by matching experimental and computed NMR spectra using the CP3 algorithm. The combination of VCD and NMR is recommended for molecules bearing more than one chiral center, as the relative configurations obtained from NMR serve as an independent check for those established with VCD. Analysis of the calculated VCD spectra reveals that the localized NH₂ scissoring mode at around 1600 cm^‐1 is characteristic for intramolecular hydrogen bonding, while the orientation of the ethyl group is reflected by the delocalized modes between 1150 and 1050 cm^‐1. Chirality 28:215–225, 2016. © 2016 Wiley Periodicals, Inc.     

1H, 15N and 13C resonance assignments of yeast Saccharomyces cerevisiae calmodulin in the Ca2+-free state

Calmodulin (CaM) is a small Ca²⁺-binding protein, which has been found in all of eucaryotic cells examined. CaMs isolated from various species have highly conserved amino acid sequence (more than 90% identical), and show the same biological functions. CaM isolated from the baker's yeast (Saccharomyces cerevisiae) (yCaM), however, shares only 60% identity in the amino acid sequence with CaM from vertebrate, and shows quite distinct conformational and biochemical properties compared with those of CaM from other species. The conformational details of yCaM, however, have not been revealed yet. We achieved the chemical shift assignments of yCaM (146 amino acids) in the apo-state using uniformly ¹⁵N- and ¹³C-labeled protein. Consequently, the resonances of 95% atoms in the backbone amides were successfully assigned.     

Peripheral membrane molecules of leukocytes and NK cytotoxicity

Some leukocyte effector cell-surface molecules movement toward the adjoining target cells takes place during the reaction of NK cytotoxicity (NK R). The majority of the moving molecules are usually anchoredvia a divalent-ion-dependent interaction (PMM-M²⁺). The released PMM-M²⁺ can interact also with the secreted tumor necrosis factor alfa (TNF-α). In agreement with PMM-M²⁺ movement, the number of TNF-α binding sites on the target cell surface increases during NK R. In addition, antibodies against PMM-M²⁺, as well asd-mannose- or N-acetyl-d-glucosamine-terminated oligosaccharides of PMM-M²⁺ inhibit NK R. A more detailed analysis of PMM-M²⁺ with monoclonal antibodies used flow cytometry and cell-surface biotinylation. Only 3 of 31 tested CD antigens (CD2, LAK-1 and CD45) were passed through this first strongly restricted experimental screening. The EDTA-released LAK-1 antigen, but not CD2 and CD45, interact with TNF-α and cell surfacevia a mannose-inhibitable interaction dependent on the presence of Ca²⁺ ions. The mechanism of possible participation of PMM-M²⁺ in cytotoxic events is discussed in relation to Ca²⁺ influx and subsequent cytolysin secretion.     

Solid-state structure of N-o-, N-m-, and N-p-nitrophenyl-2,3,4-tri-O-acetyl-β-d-xylopyranosylamines

Comprehensive structural analyses were performed for N-o-, N-m-, and N-p-nitrophenyl-2,3,4-tri-O-acetyl-β-d-xylopyranosylamines. Single-crystal X-ray diffraction data were collected and revealed that one compound under investigation undergoes temperature-dependent polymorph transitions (crystal structures of three polymorphs were obtained). The number of molecules in the independent part of the crystal unit cells was in agreement with the number of resonances in solid-state ¹³C NMR spectra. Therefore, the compounds exist as single polymorphs at room temperature, as confirmed by powder X-ray diffraction measurements. Significant differences in ¹³C chemical shifts between solution and solid-state NMR for selected carbon atoms confirmed the existence of intra- and/or intermolecular interactions.