Investigations on unconventional hydrogen bonds in RNA binding proteins: The role of CH⋯OC interactions

We have investigated the roles played by CH⋯OC interactions in RNA binding proteins. There was an average of 78 CH⋯OC interactions per protein and also there was an average of one significant CH⋯OC interaction for every 6 residues in the 59 RNA binding proteins studied. Main chain–Main chain (MM) CH⋯OC interactions are the predominant type of interactions in RNA binding proteins. The donor atom contribution to CH⋯OC interactions was mainly from aliphatic residues. The acceptor atom contribution for MM CH⋯OC interactions was mainly from Val, Phe, Leu, Ile, Arg and Ala. The secondary structure preference analysis of CH⋯OC interacting residues showed that, Arg, Gln, Glu and Tyr preferred to be in helix, while Ala, Asp, Cys, Gly, Ile, Leu, Lys, Met, Phe, Trp and Val preferred to be in strand conformation. Most of the CH⋯OC interacting polar amino acid residues were solvent exposed while, majority of the CH⋯OC interacting non polar residues were excluded from the solvent. Long and medium-range CH⋯OC interactions are the predominant type of interactions in RNA binding proteins. More than 50% of CH⋯OC interacting residues had a higher conservation score. Significant percentage of CH⋯OC interacting residues had one or more stabilization centers. Sixty-six percent of the theoretically predicted stabilizing residues were also involved in CH⋯OC interactions and hence these residues may also contribute additional stability to RNA binding proteins.     

The origin of the generalized anomeric effect: possibility of CH/n and CH/π hydrogen bonds

Ab initio MO calculations were carried out at the MP4/6-311++G(3df,3pd)//MP2/6-311++G(3df,3pd) level to investigate the conformational Gibbs energy of a series of methyl ethers CH₃O-CH₂-X (X = OH, OCH₃, F, Cl, Br, CN, CCH, C₆H₅, CHO). It was found that the Gibbs energy of the gauche conformers is lower in every case than that of the corresponding anti conformers. In the more stable gauche conformers, the interatomic distance between X and the hydrogen atom was shorter than the sum of the van der Waals radii. The natural bonding orbital (NBO) charges of group X were more negative in the gauche conformers than in the anti conformers. We suggest that the CH/n and CH/π hydrogen bonds play an important role in stabilizing the gauche conformation of these compounds.     

The role of membrane thickness in charged protein–lipid interactions

Charged amino acids are known to be important in controlling the actions of integral and peripheral membrane proteins and cell disrupting peptides. Atomistic molecular dynamics studies have shed much light on the mechanisms of membrane binding and translocation of charged protein groups, yet the impact of the full diversity of membrane physico-chemical properties and topologies has yet to be explored. Here we have performed a systematic study of an arginine (Arg) side chain analog moving across saturated phosphatidylcholine (PC) bilayers of variable hydrocarbon tail length from 10 to 18 carbons. For all bilayers we observe similar ion-induced defects, where Arg draws water molecules and lipid head groups into the bilayers to avoid large dehydration energy costs. The free energy profiles all exhibit sharp climbs with increasing penetration into the hydrocarbon core, with predictable shifts between bilayers of different thickness, leading to barrier reduction from 26 kcal/mol for 18 carbons to 6 kcal/mol for 10 carbons. For lipids of 10 and 12 carbons we observe narrow transmembrane pores and corresponding plateaus in the free energy profiles. Allowing for movements of the protein and side chain snorkeling, we argue that the energetic cost for burying Arg inside a thin bilayer will be small, consistent with recent experiments, also leading to a dramatic reduction in pK_a shifts for Arg. We provide evidence that Arg translocation occurs via an ion-induced defect mechanism, except in thick bilayers (of at least 18 carbons) where solubility-diffusion becomes energetically favored. Our findings shed light on the mechanisms of ion movement through membranes of varying composition, with implications for a range of charged protein–lipid interactions and the actions of cell-perturbing peptides. This article is part of a Special Issue entitled: Membrane protein structure and function.     

The role of protein–solvent hydrogen bond dynamics in the structural relaxation of a protein in glycerol versus water

We used MD simulations to investigate the dependence of the dynamics of a soluble protein, RNase A, on temperature and solvent environment. Consistent with neutron scattering data, the simulations predict that the protein undergoes a dynamical transition in both glycerol and aqueous solutions that is absent in the dry protein. The temperature of the transition is higher, while the rate of increase with temperature of the amplitudes of motion on the 100 ps timescale is lower, in glycerol versus water. Analysis of the dynamics of hydrogen bonds revealed that the protein dynamical transition is connected to the relaxation of the protein-solvent hydrogen bond network, which, in turn, is associated with solvent translational diffusion. Thus, it appears that the role of solvent dynamics in affecting the protein dynamical transition is qualitatively similar in water and glycerol.     

Proton transfer in and polarizability of hydrogen bonds in proteins. Tyrosine-lysine and glutamic acid-lysine hydrogen bonds — IR investigations

The OH N O^– H⁺N hydrogen bonds formed between tyrosine and lysine, and between glutamic acid and lysine residues are studied by infrared spectroscopy considering the following systems: (l-lys)_n + phenol, copoly (l-lys, l-tyr)_n, (l-lys)_n + (l-tyr)_n and (l-lys)_n + (l-glu)_n. The phenol-lysine hydrogen bonds are largely symmetrical in the average if the pK_a of the protonated lysine is 2.2 units larger than that of the phenols. In the case of the hydrogen bonds between tyrosine and lysine residues in copoly (l-lys, l-tyr)_n and (l-lys)_n + (l-tyr)_n, the weight of the proton limiting structure OH N is 80–90%, and that of the polar O^– H⁺N structure 10–20%. Double minimum proton potentials occur but the proton is preferentially present at the tyrosine residues. In the (l-lys)_n + (l-glu)_n system, the protons are present at the lysine residues. Thus, these hydrogen bonds have very large dipole moments (about 10 D). With the lysine-phenole hydrogen bonds, hydration shifts the proton transfer equilibrium a little in favour of the polar proton limiting structure O^– H⁺N. These hydrogen bonds are broken to a large extent, however, when only about 3 water molecules are present per lysine residue. When less water is present, as in the copoly (l-lys, l-tyr)_n and (l-lys)_n + (l-tyr)_n systems, these hydrogen bonds are, however, formed quantitatively. Thus — as discussed in this paper — the tyrosine-lysine hydrogen bonds can participate in proton conducting hydrogen bonded systems — as, for instance, present in bacteriorhodopsin — performing the proton transport through hydrophobic regions of biological membranes.     

The role of parent–offspring interactions during and after fledging in the Black-legged Kittiwake

Most bird species endure a high mortality at fledging, and selection should favour parental behaviour diminishing these costs. Post-fledging parental care varies greatly among species and is often linked to parent–offspring recognition. In the Black-legged Kittiwake (Rissa tridactyla), fledglings need to return to the natal nest to be fed by their parents until independence. Rejections of fledglings by non-parent adults may be fairly violent, and parents are expected to recognize and help their chicks at the time of first return. However, previous cross-fostering experiments pointed out that parents are not able to recognize their chicks up to 15 days before fledging. In this paper, we study the behaviour of both parents and juveniles at fledging. We found that parents answered significantly more to their fledgling's calls than to those of others. Compared to silent juveniles, juveniles that called before landing were more likely to be accepted by their parents. No such pattern was observed with foreign juveniles, indicating that fledglings’ voice may carry individual identity. Furthermore, fledglings found their way back to the natal nest faster when parents attended the natal nest and reacted to their offspring's calls than when they were absent or inactive. Such interactions may therefore diminish juvenile mortality at fledging.     

Enhancing effects of hydrogen/halogen bonds on σ-hole interactions involving ylide

PqsD mediates the conversion of anthraniloyl-coenzyme A (ACoA) to 2-heptyl-4-hydroxyquinoline (HHQ), a precursor of the Pseudomonas quinolone signal (PQS) molecule. Due to the role of the quinolone signaling pathway of Pseudomonas aeruginosa in the expression of several virulence factors and biofilm formation, PqsD is a potential target for controlling this nosocomial pathogen, which exhibits a low susceptibility to standard antibiotics. PqsD belongs to the β-ketoacyl-ACP synthase family and is similar in structure to homologous FabH enzymes in E. coli and Mycobacterium tuberculosis. Here, we used molecular dynamics simulations to obtain the structural position of the substrate ACoA in the binding pocket of PqsD, and semiempirical molecular orbital calculations to study the reaction mechanism for the catalytic cleavage of ACoA. Our findings suggest a nucleophilic attack of the deprotonated sulfur of Cys112 at the carbonyl carbon of ACoA and a switch in the protonation pattern of His257 whereby N_δ is protonated and the proton of N_ε is shifted to the sulfur of CoA during the reaction. This is in agreement with the experimentally determined decreased catalytic activity of the Cys112Ser mutant, whereas the Cys112Ala, His257Phe, and Asn287Ala mutants are all inactive. ESI mass-spectrometric measurements of the Asn287Ala mutant show that anthraniloyl remains covalently bound to Cys112, thus further supporting the inference from our computed mechanism that Asn287 does not take part in the cleavage of ACoA. Since this mutant is inactive, we suggest instead that Asn287 must play an essential role in the subsequent formation of HHQ in vitro.     

Control and role of pH in peptide–lipid interactions in oriented membrane samples

To understand the molecular mechanisms of amphiphilic membrane-active peptides, one needs to study their interactions with lipid bilayers under ambient conditions. However, it is difficult to control the pH of the sample in biophysical experiments that make use of mechanically aligned multilamellar membrane stacks on solid supports. HPLC-purified peptides tend to be acidic and can change the pH in the sample significantly. Here, we have systematically studied the influence of pH on the lipid interactions of the antimicrobial peptide PGLa embedded in oriented DMPC/DMPG bilayers. Using solid-state NMR (³¹P, ²H, ¹⁹F), both the lipid and peptide components were characterized independently, though in the same oriented samples under typical conditions of maximum hydration. The observed changes in lipid polymorphism were supported by DSC on multilamellar liposome suspensions. On this basis, we can present an optimized sample preparation protocol and discuss the challenges of performing solid-state NMR experiments under controlled pH. DMPC/DMPG bilayers show a significant up-field shift and broadening of the main lipid phase transition temperature when lowering the pH from 10.0 to 2.6. Both, strongly acidic and basic pH, cause a significant degree of lipid hydrolysis, which is exacerbated by the presence of PGLa. The characteristic re-alignment of PGLa from a surface-bound to a tilted state is not affected between pH of 7 to 4 in fluid bilayers. On the other hand, in gel-phase bilayers the peptide remains isotropically mobile under acidic conditions, displays various co-existing orientational states at pH 7, and adopts an unknown structural state at basic pH.     

The role of N-terminal heterocycles in hydrogen bonding to α-chymotrypsin

A series of dipeptide aldehydes containing different N-terminal heterocycles was prepared and assayed in vitro against α-chymotrypsin to ascertain the importance of the heterocycle in maintaining a β-strand geometry while also providing a hydrogen bond donor equivalent to the backbone amide nitrogen of the surrogate amino acid. The dipeptide containing a pyrrole constraint (10) was the most potent inhibitor, with >30-fold improved activity over dipeptides which lacked a nitrogen hydrogen bond donor (namely thiophene 11, furan 12 and pyridine 13). Molecular docking studies of 10 bound to α-chymotrypsin demonstrates a hydrogen bond between the pyrrole nitrogen donor and the backbone carbonyl of Gly₂₁₆ located in the S₃ pocket which is proposed to be critical for overall binding.     

Structure of pulmonary surfactant membranes and films: The role of proteins and lipid–protein interactions

The pulmonary surfactant system constitutes an excellent example of how dynamic membrane polymorphism governs some biological functions through specific lipid–lipid, lipid–protein and protein–protein interactions assembled in highly differentiated cells. Lipid–protein surfactant complexes are assembled in alveolar pneumocytes in the form of tightly packed membranes, which are stored in specialized organelles called lamellar bodies (LB). Upon secretion of LBs, surfactant develops a membrane-based network that covers rapidly and efficiently the whole respiratory surface. This membrane-based surface layer is organized in a way that permits efficient gas exchange while optimizing the encounter of many different molecules and cells at the epithelial surface, in a cross-talk essential to keep the whole organism safe from potential pathogenic invaders.The present review summarizes what is known about the structure of the different forms of surfactant, with special emphasis on current models of the molecular organization of surfactant membrane components. The architecture and the behaviour shown by surfactant structures in vivo are interpreted, to some extent, from the interactions and the properties exhibited by different surfactant models as they have been studied in vitro, particularly addressing the possible role played by surfactant proteins. However, the limitations in structural complexity and biophysical performance of surfactant preparations reconstituted in vitro will be highlighted in particular, to allow for a proper evaluation of the significance of the experimental model systems used so far to study structure–function relationships in surfactant, and to define future challenges in the design and production of more efficient clinical surfactants.     

1H-1H correlations across N–H•••N hydrogen bonds in nucleic acids

In 2HJ(NN)-COSY experiments, which correlate protons with donor/acceptor nitrogens across Nd...HNa bonds, the receptor nitrogen needs to be assigned in order to unambiguously identify the hydrogen bond. For many situations this is a non-trivial task which is further complicated by poor dispersion of (Na,Nd) resonances. To address these problems, we present pulse sequences to obtain direct, internucleotide correlations between protons in uniformly 13C/15N labeled nucleic acids containing Nd...HNa hydrogen bonds. Specifically, the pulse sequence H2(N1N3)H3 correlates H2(A,omega1):H3(U,omega2) protons across Watson-CrickA-U and mismatched G.A base pairs, the sequences H5(N3N1)H1/H6(N3N1)H1 correlate H5(C,omega1)/H6(C,omega1):H1(G,omega2) protons across Watson-Crick G-C base pairs, and the H2(N2N7)H8 sequence correlates NH2(G,A,C;omega1):H8(G,A;omega2) protons across G.G, A.A, sheared G.A and other mismatch pairs. These 1H-1H connectivities circumvent the need for independent assignment of the donor/acceptor nitrogen and related degeneracy issues associated with poorly dispersed nitrogen resonances. The methodology is demonstrated on uniformly 13C/15N labeled samples of (a) an RNA regulatory element involving the HIV-1 TAR RNA fragment, (b) a multi-stranded DNA architecture involving a G.(C-A) triad-containing G-quadruplex and (c) a peptide-RNA complex involving an evolved peptide bound to the HIV-1 Rev response element (RRE) RNA fragment.     

NMR experiments for the rapid identification of P=O···H–X type hydrogen bonds in nucleic acids

Implementation of a new algorithm, SMILE, is described for reconstruction of non-uniformly sampled two-, three- and four-dimensional NMR data, which takes advantage of the known phases of the NMR spectrum and the exponential decay of underlying time domain signals. The method is very robust with respect to the chosen sampling protocol and, in its default mode, also extends the truncated time domain signals by a modest amount of non-sampled zeros. SMILE can likewise be used to extend conventional uniformly sampled data, as an effective multidimensional alternative to linear prediction. The program is provided as a plug-in to the widely used NMRPipe software suite, and can be used with default parameters for mainstream application, or with user control over the iterative process to possibly further improve reconstruction quality and to lower the demand on computational resources. For large data sets, the method is robust and demonstrated for sparsities down to ca 1%, and final all-real spectral sizes as large as 300 Gb. Comparison between fully sampled, conventionally processed spectra and randomly selected NUS subsets of this data shows that the reconstruction quality approaches the theoretical limit in terms of peak position fidelity and intensity. SMILE essentially removes the noise-like appearance associated with the point-spread function of signals that are a default of five-fold above the noise level, but impacts the actual thermal noise in the NMR spectra only minimally. Therefore, the appearance and interpretation of SMILE-reconstructed spectra is very similar to that of fully sampled spectra generated by Fourier transformation.     

The roles of ABA in plant–pathogen interactions

Defence against abiotic and biotic stresses is crucial for the fitness and survival of plants under adverse or suboptimal growth conditions. The phytohormone abscisic acid (ABA) is not only important for mediating abiotic stress responses, but also plays a multifaceted and pivotal role in plant immunity. This review presents examples demonstrating the importance of crosstalk between ABA and the key biotic stress phytohormone salicylic acid in determining the outcome of plant–pathogen interactions. We then provide an overview of how ABA influences plant defence responses against various phytopathogens with particular emphasis on the Arabidopsis–Pseudomonas syringae model pathosystem. Lastly, we discuss future directions for studies of ABA in plant immunity with emphasis on, its role in the crosstalk between biotic and abiotic stress responses, the importance of distinguishing direct and indirect effects of ABA, as well as the prospect of utilizing the recently elucidated core ABA signaling network to gain further insights into the roles of ABA in plant immunity.     

The expanding role of miR-302–367 in pluripotency and reprogramming

MicroRNA (miRNA) has been shown to be essential for regulating cell fate and pluripotency; however, our knowledge of miRNA function in stem cells is incomplete due to experimental limitations and difficulties in identifying their physiological targets. Recent studies implicated hESC-expressed miRNAs (miR?302–367 and miR?371–373 clusters) in regulating BMP signaling and promoting pluripotency, suggesting that low levels of BMP signaling may promote pluripotency by preventing neural induction. A comprehensive list of miR?302–367 targets recently identified by genome-wide approaches suggests a number of additional cellular processes and signaling pathways whose regulation by miR?302–367 may promote pluripotency and reprogramming, such as cell cycle, epigenetic changes, metabolism and vesicular transfer.     

Intramolecular CH···O hydrogen bonds in the AI and BI DNA-like conformers of canonical nucleosides and their Watson-Crick pairs. Quantum chemical and AIM analysis

The aim of this work is to cast some light on the H-bonds in double-stranded DNA in its AI and BI forms. For this purpose, we have performed the MP2 and DFT quantum chemical calculations of the canonical nucleoside conformers, relative to the AI and BI DNA forms, and their Watson-Crick pairs, which were regarded as the simplest models of the double-stranded DNA. Based on the atoms-in-molecules analysis (AIM), five types of the CH···O hydrogen bonds, involving bases and sugar, were detected numerically from 1 to 3 per a conformer: C2'H···O5', C1'H···O2, C6H···O5', C8H···O5', and C6H···O4'. The energy values of H-bonds occupy the range of 2.3-5.6 kcal/mol, surely exceeding the kT value (0.62 kcal/mol). The nucleoside CH···O hydrogen bonds appeared to "survive" turns of bases against the sugar, sometimes in rather large ranges of the angle values, pertinent to certain conformations, which points out to the source of the DNA lability, necessary for the conformational adaptation in processes of its functioning. The calculation of the interactions in the dA·T nucleoside pair gives evidence, that additionally to the N6H···O4 and N1···N3H canonical H-bonds, between the bases adenine and thymine the third one (C2H···O2) is formed, which, though being rather weak (about 1 kcal/mol), satisfies the AIM criteria of H-bonding and may be classified as a true H-bond. The total energy of all the CH···O nontraditional intramolecular H-bonds in DNA nucleoside pairs appeared to be commensurable with the energy of H-bonds between the bases in Watson-Crick pairs, which implies their possible important role in the DNA shaping.