Antiparallel triple helices. Structural characteristics and stabilization by 8-amino derivatives

The structural, dynamical, and recognition properties of antiparallel DNA triplexes formed by the antiparallel d(G#G.C), d(A#A.T), and d(T#A.T) motifs (the pound sign and dot mean reverse-Hoogsteen and Watson-Crick hydrogen bonds, respectively) are studied by means of "state of the art" molecular dynamics simulations. Once the characteristics of the helix are defined, molecular dynamics and thermodynamic integration calculations are used to determine the expected stabilization of the antiparallel triplex caused by the introduction of 8-aminopurines. Finally, oligonucleotides containing 8-aminopurine derivatives are synthesized and tested experimentally using several approaches in a variety of systems. A very large stabilization of the triplex is found experimentally, as predicted by simulations. These results open the possibility for the use of oligonucleotides carrying 8-aminopurines to bind single-stranded nucleic acids by formation of antiparallel triplexes.     

TreeDock: a tool for protein docking based on minimizing van der Waals energies

Predicting protein-protein and protein-ligand docking remains one of the challenging topics of structural biology. The main problems are (i) to reliably estimate the binding free energies of docked states, (ii) to enumerate possible docking orientations at a high resolution, and (iii) to consider mobility of the docking surfaces and structural rearrangements upon interaction. Here we present a novel algorithm, TreeDock, that addresses the enumeration problem in a rigid-body docking search. By representing molecules as multidimensional binary search trees and by exploring a sufficient number of docking orientations such that two chosen atoms, one from each molecule, are always in contact, TreeDock is able to explore all clash-free orientations at very fine resolution in a reasonable amount of time. Due to the speed of the program, many contact pairs can be examined to search partial or complete surface areas. The deterministic systematic search of TreeDock is in contrast to most other docking programs that use stochastic searches such as Monte Carlo or simulated annealing methods. At this point, we have used the Lennard-Jones potential as the only scoring function and show that this can predict the correct docked conformation for a number of protein-protein and protein-ligand complexes. The program is most powerful if some information is known about the location of binding faces from NMR chemical-shift perturbation studies, orientation information from residual dipolar coupling, or mutational screening. The approach has the potential to include docking-site mobility by performing molecular dynamics or other randomization methods of the docking site and docking families to families of structures. The performance of the algorithm is demonstrated by docking three complexes of immunoglobulin superfamily domains, CD2 to CD58, the V(alpha) domain of a T-cell receptor to its V(beta) domain, and a T-cell receptor to a pMHC complex as well as a small molecule inhibitor to a phosphatase.     

Diastereoselective synthesis of psi[(E)-CH=CMe]- and psi[(Z)-CH=CMe]-type dipeptide isosteres by organocopper-mediated anti-S(N)2' reaction

[reaction: see text] Acyclic psi[(E)-CH=CMe]- and psi[(Z)-CH=CMe]-type dipeptide isosteres were efficiently synthesized. In a key reaction, alpha-alkylation of gamma-mesyloxy-beta-methyl-alpha,beta-unsaturated esters with organocyanocuprates in diethyl ether or tetrahydrofuran preferentially afforded the psi[(E)-CH=CMe]- or psi[(Z)-CH=CMe]-isomer, respectively, via anti-S(N)2' mechanism.     

Free-energy component analysis of 40 protein-DNA complexes: a consensus view on the thermodynamics of binding at the molecular level

Noncovalent association of proteins to specific target sites on DNA--a process central to gene expression and regulation--has thus far proven to be idiosyncratic and elusive to generalizations on the nature of the driving forces. The spate of structural information on protein--DNA complexes sets the stage for theoretical investigations on the molecular thermodynamics of binding aimed at identifying forces responsible for specific macromolecular recognition. Computation of absolute binding free energies for systems of this complexity transiting from structural information is a stupendous task. Adopting some recent progresses in treating atomic level interactions in proteins and nucleic acids including solvent and salt effects, we have put together an energy component methodology cast in a phenomenological mode and amenable to systematic improvements and developed a computational first atlas of the free energy contributors to binding in approximately 40 protein-DNA complexes representing a variety of structural motifs and functions. Illustrating vividly the compensatory nature of the free energy components contributing to the energetics of recognition for attaining optimal binding, our results highlight unambiguously the roles played by packing, electrostatics including hydrogen bonds, ion and water release (cavitation) in protein-DNA binding. Cavitation and van der Waals contributions without exception favor complexation. The electrostatics is marginally unfavorable in a consensus view. Basic residues on the protein contribute favorably to binding despite the desolvation expense. The electrostatics arising from the acidic and neutral residues proves unfavorable to binding. An enveloping mode of binding to short stretches of DNA makes for a strong unfavorable net electrostatics but a highly favorable van der Waals and cavitation contribution. Thus, noncovalent protein-DNA association is a system-specific fine balancing act of these diverse competing forces. With the advances in computational methods as applied to macromolecular recognition, the challenge now seems to be to correlate the differential (initial vs. final) energetics to substituent effects in drug design and to move from affinity to specificity.     

Cooling rate dependence of specific heat in systems out of equilibrium

The anomaly of specific heat in systems out of equilibrium, especially the measurement procedure dependence of specific heat, is investigated by means of free energy landscape. Introducing measurement procedure which is based on experimental method, we propose a calculation method of specific heat in systems out of equilibrium and find an abrupt change in specific heat between annealed and quenched states. For longer observation time the change in specific heat occurs at lower temperature and becomes sharper. For slower cooling of a system the transition temperature becomes lower. This cooling rate dependence of the transition temperature is consistent with experiments and thus the abrupt change in specific heat can be regarded as the glass transition which is thermally identified.     

Adhesion of bacterial exopolymers to alpha-FeOOH: inner-sphere complexation of phosphodiester groups

Extracellular polymeric substances (EPS) constitute a heterogeneous mixture of polyelectrolytes that mediate biomineralization and bacterial adhesion and stabilize biofilm matrixes in natural and artificial environments. Although nucleic acids are exuded extracellularly and are purported to be required for biofilm formation, direct evidence of the active mechanism is lacking. EPS were extracted from both Bacillus subtilis (a gram-positive bacterium) and Pseudomonas aeruginosa (a gram-negative bacterium) and their interaction with the goethite (alpha-FeOOH) surface was studied using attenuated total internal reflection infrared spectroscopy. Correspondence between spectral data and quantum chemical calculations demonstrate that phosphodiester groups of nucleic acids mediate the binding of EPS to mineral surfaces. Our data indicate that these groups emerge from the EPS mixture to form monodentate complexes with Fe centers on the goethite (alpha-FeOOH) surface, providing an energetically stable bond for further EPS or cell adhesion.     

The use of solid-phase supported 1-N-piperazine-4-N-carboxaldehyde in Vilsmeier reactions

A Vilsmeier salt supported on solid phase was prepared using piperazine bound to Merrifield resin. Piperazine was selected because it contains two secondary amines: one of the amines is protected upon binding to the resin, and the second was formylated to give resin-1-N-piperazine-4-N-carboxaldehyde (9). Activation of the formamide with either bis(trichloromethyl)carbonate (BTC) or POCl(3) afforded the Vilsmeier salt 10. Several olefins were used to test the supported Vilsmeier reagent. The in-solution activation with BTC and POCl(3) of various secondary amides was also evaluated: dimethylformamide (1), N-methyformanilide (4), 4-formylmorpholine (5), and 1,4-dicarboxylpiperazine (6), which showed that amides with one additional heteroatom increase the yields in the Vilsmeier salt formation.     

Selective metal promoted hydrolysis of nitrile groups in the side chain of tetraazamacrocyclic Cu2+-complexes

The metal promoted hydrolysis of nitrile groups in the side chains of tetraazamacrocyclic Cu2+ complexes has been studied by stopped-flow techniques. It is shown that the reaction proceeds by an intramolecular attack of an axially coordinated OH- onto the nitrile group to give the corresponding amide. In alkaline solution the amide then deprotonates and binds to the axial position of the Cu2+ thus preventing further coordination of an OH-. This explains mechanistically that in the Cu2+ complexes of macrocycles carrying two nitrile functions only one is selectively hydrolysed. The nitrile hydrolysis has also been used on a preparative scale to synthesize tetraazamacrocycles with two different side chains. X-Ray diffractions of several products are presented to confirm the structures and the results from the kinetics and equilibria measurements.