Computational Study of Peptide Plane Stacking with Polar and Ionizable Amino Acid Side Chains

Urea-aromatic interactions in biology

Noncovalent interactions are key determinants in both chemical and biological processes. Among such processes, the hydrophobic interactions play an eminent role in folding of proteins, nucleic acids, formation of membranes, protein-ligand recognition, etc.. Though this interaction is mediated through the aqueous solvent, the stability of the above biomolecules can be highly sensitive to any small external perturbations, such as temperature, pressure, pH, or even cosolvent additives, like, urea-a highly soluble small organic molecule utilized by various living organisms to regulate osmotic pressure. A plethora of detailed studies exist covering both experimental and theoretical regimes, to understand how urea modulates the stability of biological macromolecules. While experimentalists have been primarily focusing on the thermodynamic and kinetic aspects, theoretical modeling predominantly involves mechanistic information at the molecular level, calculating atomistic details applying the force field approach to the high level electronic details using the quantum mechanical methods. The review focuses mainly on examples with biological relevance, such as (1) urea-assisted protein unfolding, (2) urea-assisted RNA unfolding, (3) urea lesion interaction within damaged DNA, (4) urea conduction through membrane proteins, and (5) protein-ligand interactions those explicitly address the vitality of hydrophobic interactions involving exclusively the urea-aromatic moiety.

General Recognition of U-G, U-A, and C-G Pairs by Double-Stranded RNA-Binding PNAs Incorporated with an Artificial Nucleobase

Chemically modified peptide nucleic acids (PNAs) show great promise in the recognition of RNA duplexes by major-groove PNA·RNA-RNA triplex formation. Triplex formation is favored for RNA duplexes with a purine tract within one of the RNA duplex strands, and is severely destabilized if the purine tract is interrupted by pyrimidine residues. Here, we report the synthesis of a PNA monomer incorporated with an artificial nucleobase S, followed by the binding studies of a series of S-modified PNAs. Our data suggest that an S residue incorporated into short 8-mer dsRNA-binding PNAs (dbPNAs) can recognize internal Watson-Crick C-G and U-A, and wobble U-G base pairs (but not G-C, A-U, and G-U pairs) in RNA duplexes. The short S-modified PNAs show no appreciable binding to DNA duplexes or single-stranded RNAs. Interestingly, replacement of the C residue in an S·C-G triple with a 5-methyl C results in the disruption of the triplex, probably due to a steric clash between S and 5-methyl C. Previously reported PNA E base shows recognition of U-A and A-U pairs, but not a U-G pair. Thus, S-modified dbPNAs may be uniquely useful for the general recognition of RNA U-G, U-A, and C-G pairs. Shortening the succinyl linker of our PNA S monomer by one carbon atom to have a malonyl linker causes a severe destabilization of triplex formation. Our experimental and modeling data indicate that part of the succinyl moiety in a PNA S monomer may serve to expand the S base forming stacking interactions with adjacent PNA bases.

Arginine side chain stacking with peptide plane stabilizes the protein helix conformation in a cooperative way

A combined experimental and computational study is performed for arginine side chain stacking with the protein α-helix. Theremostability measurements of Aristaless homeodomain, a helical protein, suggest that mutating the arginine residue R106, R137 or R141, which has the guanidino side chain stacking with the peptide plane, to alanine, destabilizes the protein. The R-PP stacking has an energy of ∼0.2-0.4 kcal/mol. This stacking interaction mainly comes from dispersion and electrostatics, based on MP2 calculations with the energy decomposition analysis. The calculations also suggest that the stacking stabilizes 2 backbone-backbone h-bonds (i→i-4 and i-3→i-7) in a cooperative way. Desolvation and electrostatic polarization are responsible for cooperativity with the i→i-4 and i-3→i-7 h-bonds, respectively. This cooperativity is supported by a protein α-helices h-bond survey in the pdb databank where stacking shortens the corresponding h-bond distances.

Unravelling hydrogen bonding interactions of tryptamine–water dimer from neutral to cation

The physical origin for the three intermolecular hydrogen bonds in the neutral and cationic forms of the tryptamine–water dimer is explored.

An Evolving Insight into Chiral H-Bond Catalyzed Aza-Henry Reactions: A Cooperative Role for Noncovalent Attractive Interactions Unveiled by Density Functional Theory

The role of cooperative effects arising from noncovalent attractive interactions as a vital factor governing stereoinduction in chiral H-bond catalyzed aza-Henry reactions is reported. Supporting this finding were density functional theory (DFT) calculations which revealed a shape and size dependency existed between the catalyst and substrates that when matched lead to high enantioselectivity, as reflected by favorable activation parameters. Associated with optimal catalyst and substrate pairing were a closed catalytic binding pocket and a synclinal orientation of the substrates that reinforced favorable stereoelectronic effects and dispersive type forces. Meanwhile, unfavorable steric interactions were found to be a dominant effect controlling diastereoselection.