Conservation and diversification in homeodomain-DNA interactions: a comparative genetic analysis.

Nearly all metazoan homeodomains (HDs) possess DNA binding targets that are related by the presence of a TAAT sequence. We use an in vitro genetic DNA binding site selection assay to refine our understanding of the amino acid determinants for the recognition of the TAAT site. Superimposed upon the conserved ability of metazoan HDs to recognize a TAAT core is a difference in their preference for the bases that lie immediately 3' to it. Amino acid position 50 of the HD has been shown to discriminate among these base pairs, and structural studies have suggested that water-mediated hydrogen bonds and van der Waals contacts underlie for this ability. Here, we show that each of six amino acids tested at position 50 can confer a distinct DNA binding specificity.

Local densities orthogonal to beta-sheet amide planes: patterns of packing in globular proteins.

We have investigated the efficiency of packing by calculating intramolecular packing density above and below peptide planes of internal beta-pleated sheet residues in five globular proteins. The orientation of interest was chosen to allow study of regions that are approximately perpendicular to the faces of beta-pleated sheets. In these locations, nonbonded van der Waals packing interactions predominate over hydrogen bonding and solvent interactions. We observed considerable variability in packing densities within these regions, confirming that the interior packing of a protein does not result in uniform occupation of the available space. Patterns of fluctuation in packing density suggest that the regular backbone-to-backbone network of hydrogen bonds is not likely to be interrupted to maximize van der Waals interactions. However, high-density packing tends to occur toward the ends of beta-structure strands where hydrogen bonds are more likely to involve nonpolar side-chain groups or solvent molecules. These features result in internal protein folding with a central low-density core surrounded by a higher-density subsurface shell, consistent with our previous calculations regarding overall protein packing density.

Least activation path for protein folding: investigation of staphylococcal nuclease folding by stopped-flow circular dichroism.

Is the pathway of protein folding determined by the relative stability of folding intermediates, or by the relative height of the activation barriers leading to these intermediates? This is a fundamental question for resolving the Levinthal paradox, which stated that protein folding by a random search mechanism would require a time too long to be plausible. To answer this question, we have studied the guanidinium chloride (GdmCl)-induced folding/unfolding of staphylococcal nuclease [(SNase, formerly EC 3.1.4.7; now called microbial nuclease or endonuclease, EC 3.1.31.1] by stopped-flow circular dichroism (CD) and differential scanning microcalorimetry (DSC). The data show that while the equilibrium transition is a quasi-two-state process, kinetics in the 2-ms to 500-s time range are triphasic. Data support the sequential mechanism for SNase folding: U3 <--> U2 <--> U1 <--> N0, where U1, U2, and U3 are substates of the unfolded protein and N0 is the native state. Analysis of the relative population of the U1, U2, and U3 species in 2.0 M GdmCl gives delta-G values for the U3 --> U2 reaction of +0.1 kcal/mol and for the U2 --> U1 reaction of -0.49 kcal/mol. The delta-G value for the U1 --> N0 reaction is calculated to be -4.5 kcal/mol from DSC data. The activation energy, enthalpy, and entropy for each kinetic step are also determined. These results allow us to make the following four conclusions. (i) Although the U1, U2, and U3 states are nearly isoenergetic, no random walk occurs among them during the folding. The pathway of folding is unique and sequential. In other words, the relative stability of the folding intermediates does not dictate the folding pathway. Instead, the folding is a descent toward the global free-energy minimum of the native state via the least activation path in the vast energy landscape. Barrier avoidance leads the way, and barrier height limits the rate. Thus, the Levinthal paradox is not applicable to the protein-folding problem. (ii) The main folding reaction (U1 --> N0), in which the peptide chain acquires most of its free energy (via van der Waals' contacts, hydrogen bonding, and electrostatic interactions), is a highly concerted process. These energy-acquiring events take place in a single kinetic phase. (iii) U1 appears to be a compact unfolded species; the rate of conversion of U2 to U1 depends on the viscosity of solution. (iv) All four relaxation times reported here depend on GdmCl concentrations: it is likely that none involve the cis/trans isomerization of prolines. Finally, a mechanism is presented in which formation of sheet-like chain conformations and a hydrophobic condensation event precede the main-chain folding reaction.

Molecular mechanism of protein-retinal coupling in bacteriorhodopsin.

Bacteriorhodopsin is a membrane protein that functions as a light-driven proton pump. Each cycle of proton transport is initiated by the light-induced isomerization of retinal from the all-trans to 13-cis configuration and is completed by the protein-driven reisomerization of retinal to the all-trans configuration. Previous studies have shown that replacement of Leu-93, a residue in close proximity to the 13-methyl group of retinal, by alanine, resulted in a 250-fold increase in the time required to complete each photocycle. Here, we show that the kinetic defect in the photocycle of the Leu-93-->Ala mutant occurs at a stage after the completion of proton transport and can be overcome in the presence of strong background illumination. Time-resolved retinal-extraction experiments demonstrate the continued presence of a 13-cis intermediate in the photocycle of the Leu-93-->Ala mutant well after the completion of proton release and uptake. These results indicate that retinal reisomerization is kinetically the rate-limiting step in the photocycle of this mutant and that the slow thermal reisomerization can be bypassed by the absorption of a second photon. The effects observed for the Leu-93-->Ala mutant are not observed upon replacement of any other residue in van der Waals contact with retinal or upon replacement of Leu-93 by valine. We conclude that the contact between Leu-93 and the 13-methyl group of retinal plays a key role in controlling the rate of protein conformational changes associated with retinal reisomerization and return of the protein to the initial state.

Calicheamicin-DNA complexes: warhead alignment and saccharide recognition of the minor groove.

The solution structures of calicheamicin gamma 1I, its cycloaromatized analog (calicheamicin epsilon), and its aryl tetrasaccharide complexed to a common DNA hairpin duplex have been determined by NMR and distance-refined molecular dynamics computations. Sequence specificity is associated with carbohydrate-DNA recognition that places the aryl tetrasaccharide component of all three ligands in similar orientations in the minor groove at the d(T-C-C-T).d(A-G-G-A) segment. The complementary fit of the ligands and the DNA minor groove binding site creates numerous van der Waals contacts as well as hydrogen bonding interactions. Notable are the iodine and sulfur atoms of calicheamicin that hydrogen bond with the exposed amino proton of the 5'- and 3'-guanines, respectively, of the d(A-G-G-A) segment. The sequence-specific carbohydrate binding orients the enediyne aglycone of calicheamicin gamma 1I such that its C3 and C6 proradical centers are adjacent to the cleavage sites. While the enediyne aglycone of calicheamicin gamma 1I is tilted relative to the helix axis and spans the minor groove, the cycloaromatized aglycone is aligned approximately parallel to the helix axis in the respective complexes. Specific localized conformational perturbations in the DNA have been identified from imino proton complexation shifts and changes in specific sugar pucker patterns on complex formation. The helical parameters for the carbohydrate binding site are comparable with corresponding values in B-DNA fibers while a widening of the groove is observed at the adjacent aglycone binding site.

The old problems of glass and the glass transition, and the many new twists.

In this paper I review the ways in which the glassy state is obtained both in nature and in materials science and highlight a "new twist"--the recent recognition of polymorphism within the glassy state. The formation of glass by continuous cooling (viscous slowdown) is then examined, the strong/fragile liquids classification is reviewed, and a new twist-the possibility that the slowdown is a result of an avoided critical point-is noted. The three canonical characteristics of relaxing liquids are correlated through the fragility. As a further new twist, the conversion of strong liquids to fragile liquids by pressure-induced coordination number increases is demonstrated. It is then shown that, for comparable systems, it is possible to have the same conversion accomplished via a first-order transition within the liquid state during quenching. This occurs in the systems in which "polyamorphism" (polymorphism in the glassy state) is observed, and the whole phenomenology is accounted for by Poole's bond-modified van der Waals model. The sudden loss of some liquid degrees of freedom through such weak first-order transitions is then related to the polyamorphic transition between native and denatured hydrated proteins, since the latter are also glass-forming systems--water-plasticized, hydrogen bond-cross-linked chain polymers (and single molecule glass formers). The circle is closed with a final new twist by noting that a short time scale phenomenon much studied by protein physicists-namely, the onset of a sharp change in d/dT ( is the Debye-Waller factor)--is general for glass-forming liquids, including computer-simulated strong and fragile ionic liquids, and is closely correlated with the experimental glass transition temperature. The latter thus originates in strong anharmonicity in certain components of the vibrational density of states, which permits the system to access the multiple minima of its configuration space. The connection between the anharmonicity in these modes, vibrational localization, the Kauzmann temperature, and the fragility of the liquid is proposed as the key problem in glass science.

Nanoparticles of tungsten as low-cost monometallic catalyst for selective hydrogenation of 3-hexyne

Low-cost tungsten monometallic catalysts containing variable amounts of metal (4.5, 7.1 and 8.5%W) were prepared by impregnating alumina with ammonium metatungstate as an inexpensive precursor. The catalysts were characterized using ICP, XPS, XRD, TPR and hydrogen chemisorption. These techniques revealed mainly WO3-Al2O3 (W6+) species on the surface. The effects of the content of W nanoparticles and reaction temperature on activity and selectivity for the partial hydrogenation of 3-hexyne, a non-terminal alkyne, were assessed under moderate conditions of temperature and pressure. The monometallic catalysts prepared were found to be active and stereoselective for the production of (Z )-3-hexene, had the following order: 7.1WN/A > 8.5 WN/A ≥ 4.5 WN/A. Additionally, the performance of the synthesized xWN/A catalysts exhibited high sensitivity to temperature variation. In all cases, the maximum 3-hexyne total conversion and selectivity was achieved at 323 K. The performance of the catalysts was considered to be a consequence of two phenomena: a) the electronic effects, related to the high charge of W (+6), causing an intensive dipole moment in the hydrogen molecule (van der Waals forces) and leading to heterolytic bond rupture; the H+ and H- species generated approach a 3-hexyne adsorbate molecule and cause heterolytic rupture of the C≡C bond into C- = C+; and b) steric effects related to the high concentration of WO3 on 8.5WN/A that block the Al2O3 support. Catalyst deactivation was detected, starting at about 50 min of reaction time. Electrodeficient W6+ species are responsible for the formation of green oil at the surface level, blocking pores and active sites of the catalyst, particularly at low reaction temperatures (293 and 303 K). The resulting best catalyst, 7.1WN/A, has low fabrication cost and high selectivity for (Z )-3-hexene (94%) at 323 K. This selectivity is comparable to that of the classical and more expensive industrial Lindlar catalyst (5 wt% Pd). The alumina supported tungsten catalysts are low-cost potential replacements for the Lindlar industrial catalyst. These catalysts could also be used for preparing bimetallic W-Pd catalysts for selective hydrogenation of terminal and non-terminal alkynes.

Interfacial interactions and structure of polyurethane intercalated nanocomposite

Understanding the interfacial interactions between the nanofiller and polymer matrix is important to improve the design and manufacture of polymer nanocomposites. This paper reports a molecular dynamic Study on the interfacial interactions and structure of a clay-based polyurethane intercalated nanocomposite. The results show that the intercalation of surfactant (i.e. dioctadecyldlmethyl ammonium) and polyurethane (PU) into the nanoconfined gallery of clay leads to the multilayer structure for both surfactant and PU, and the absence of phase separation for PU chains. Such structural characteristics are attributed to the result of competitive interactions among the surfactant, PU and the clay surface, including van der Waals, electrostatic and hydrogen bonding.

Molecular dynamics characterization of n-octyl-beta-D-glucopyranoside micelle structure in aqueous solution

n-Octyl-beta-D-glueopyranoside (OG) is a non-ionic glycolipid, which is used widely in biotechnical and biochemical applications. All-atom molecular dynamics simulations from two different initial coordinates and velocities in explicit solvent have been performed to characterize the structural behaviour of an OG aggregate at equilibrium conditions. Geometric packing properties determined from the simulations and small angle neutron scattering experiment state that OG micelles are more likely to exist in a non-spherical shape, even at the concentration range near to the critical micelle concentration (0.025 M). Despite few large deviations in the principal moment of inertia ratios, the average micelle shape calculated from both simulations is a prolate ellipsoid. The deviations at these time scales are presumably the temporary shape change of a micelle. However, the size of the micelle and the accessible surface areas were constant during the simulations with the micelle surface being rough and partially elongated. Radial distribution functions computed for the hydroxyl oxygen atoms of an OG show sharper peaks at a minimum van der Waals contact distance than the acetal oxygen, ring oxygen, and anomeric carbon atoms. This result indicates that these atoms are pointed outwards at the hydrophilic/hydrophobic interface, form hydrogen bonds with the water molecules, and thus hydrate the micelle surface effectively. (c) 2005 Elsevier Inc. All rights reserved.

Structure, dynamics, and energetics of siRNA-cationic vector complexation:a molecular dynamics study

The design and synthesis of safe and efficient nonviral vectors for gene delivery has attracted significant attention in recent years. Previous experiments have revealed that the charge density of a polycation (the carrier) plays a crucial role in complexation and the release of the gene from the complex in the cytosol. In this work, we adopt an atomistic molecular dynamics simulation approach to study the complexation of short strand duplex RNA with six cationic carrier systems of varying charge and surface topology. The simulations reveal detailed molecular-level pictures of the structures and dynamics of the RNA-polycation complexes. Estimates for the binding free energy indicate that electrostatic contributions are dominant followed by van der Waals interactions. The binding free energy between the 8(+)polymers and the RNA is found to be larger than that of the 4(+)polymers, in general agreement with previously published data. Because reliable binding free energies provide an effective index of the ability of the polycationic carrier to bind the nucleic acid and also carry implications for the process of gene release within the cytosol, these novel simulations have the potential to provide us with a much better understanding of key mechanistic aspects of gene-polycation complexation and thereby advance the rational design of nonviral gene delivery systems.

Solvent induced NMR chemical shifts that arise from molecular encounters

Recently Homer and Percival have postulated that intermolecular van der Waals dispersion forces can be characterized by three mechanisms. The first arises via the mean square reaction field < R1; 2> due to the transient dipole of a particular solute molecule that is considered situated in a cavity surrounded by solvent molecules; this was characterized by an extended Onsager approach. The second stems from the extra cavity mean square reaction field < R2; 2> of the near neighbour solvent molecules. The third originates from square field electric fields E2BI due to a newly characterized effect in which solute atoms are `buffeted' by the peripheral atoms of adjacent solvent molecules. The present work concerns more detailed studies of the buffeting screening, which is governed by sterically controlled parameter (2T - T)2, where and are geometric structural parameters. The original approach is used to characterise the buffeting shifts induced by large solvent molecules and the approach is found to be inadequate. Consequently, improved methods of calculating and are reported. Using the improved approach it is shown that buffeting is dependent on the nature of the solvent as well as the nature of the solute molecule. Detailed investigation of the buffeting component of the van der Waals chemical shifts of selected solutes in a range of solvents containing either H or Cl as peripheral atoms have enabled the determination of a theoretical acceptable value for the classical screening coefficient B for protons. 1H and 13C resonance studies of tetraethylmethane and 1H, 13C and 29Si resonance studies of TMS have been used to support the original contention that three (< R1; 2> , < R2; 2> and E2BI) components of intermolecular van der Waals dispersion fields are required to characterise vdW chemical shifts.

Viscoelasticity and high buckling stress of dense carbon nanotube brushes

We report on the mechanical behavior of a dense brush of small-diameter (1–3 nm) non-catalytic multiwall (2–4 walls) carbon nanotubes (CNTs), with ~10 times higher density than CNT brushes produced by other methods. Under compression with spherical indenters of different radii, these highly dense CNT brushes exhibit a higher modulus (~17–20 GPa) and orders of magnitude higher resistance to buckling than vapor phase deposited CNT brushes or carbon walls. We also demonstrate the viscoelastic behavior, caused by the increased influence of the van der Waals’ forces in these highly dense CNT brushes, showing their promise for energy-absorbing coatings.

The structural conversion from α-AgVO3 to β-AgVO3: Ag nanoparticle decorated nanowires with application as cathode materials for Li-ion batteries

The majority of electrode materials in batteries and related electrochemical energy storage devices are fashioned into slurries via the addition of a conductive additive and a binder. However, aggregation of smaller diameter nanoparticles in current generation electrode compositions can result in non-homogeneous active materials. Inconsistent slurry formulation may lead to inconsistent electrical conductivity throughout the material, local variations in electrochemical response, and the overall cell performance. Here we demonstrate the hydrothermal preparation of Ag nanoparticle (NP) decorated α-AgVO3 nanowires (NWs) and their conversion to tunnel structured β-AgVO3 NWs by annealing to form a uniform blend of intercalation materials that are well connected electrically. The synthesis of nanostructures with chemically bound conductive nanoparticles is an elegant means to overcome the intrinsic issues associated with electrode slurry production, as wire-to-wire conductive pathways are formed within the overall electrode active mass of NWs. The conversion from α-AgVO3 to β-AgVO3 is explained in detail through a comprehensive structural characterization. Meticulous EELS analysis of β-AgVO3 NWs offers insight into the true β-AgVO3 structure and how the annealing process facilitates a higher surface coverage of Ag NPs directly from ionic Ag content within the α-AgVO3 NWs. Variations in vanadium oxidation state across the surface of the nanowires indicate that the β-AgVO3 NWs have a core–shell oxidation state structure, and that the vanadium oxidation state under the Ag NP confirms a chemically bound NP from reduction of diffused ionic silver from the α-AgVO3 NWs core material. Electrochemical comparison of α-AgVO3 and β-AgVO3 NWs confirms that β-AgVO3 offers improved electrochemical performance. An ex situ structural characterization of β-AgVO3 NWs after the first galvanostatic discharge and charge offers new insight into the Li+ reaction mechanism for β-AgVO3. Ag+ between the van der Waals layers of the vanadium oxide is reduced during discharge and deposited as metallic Ag, the vacant sites are then occupied by Li+.

Understanding catalytic reactions over zeolites: A density functional theory study of selective catalytic reduction of NOX by NH3 over Cu-SAPO-34

Metal exchanged CHA-type (SAPO-34 and SSZ-13) zeolites are promising catalysts for selective catalytic reduction (SCR) of NOx by NH3. However, the understanding of the process at the molecular level is still limited, which hinders the identification of its mechanism and the design of more efficient zeolite catalysts. In this work, modelling the reaction over Cu-SAPO-34, a periodic density functional theory (DFT) study of NH3-SCR was performed using hybrid functional with the consideration of van der Waals (vdW) interactions. A mechanism with a low N–N coupling barrier is proposed to account for the activation of NO. The redox cycle of Cu2+ and Cu+, which is crucial for the SCR process, is identified with detailed analyses. Besides, the decomposition of NH2NO is shown to readily occur on the Brønsted acid site by a hydrogen push-pull mechanism, confirming the collective efforts of Brønsted acid and Lewis acid (Cu2+) sites. The special electronic and structural properties of Cu-SAPO-34 are demonstrated to play an essential role the reaction, which may have a general implication on the understanding of zeolite catalysis.

Toward controlled growth of helicity-specific carbon nanotubes

The underlying mechanisms for the nucleation of carbon nanotubes as well as their helicity, remain elusive. Here, using van der Waals dispersion force calculations implemented within density functional theory, we study the cap formation, believed to be responsible for the chirality of surface-catalyzed carbon nanotubes. We find the energetics associated with growth along different facets to be independent of the surface orientation and that the growth across an edge along the axis of the metal particle leads to a perfect honeycomb lattice in a curved geometry. The formation of defects in the graphene matrix, which bend the carbon plane, requires that two or more graphene embryos with significantly different growth axis merge. Such scenario is only possible at the front- or back-end of the metal particle where growth symmetry is broken. The graphene embryos reconstruct their hexagonal structure into pentagons, heptagons, and octagons counterpart to accommodate the tube curvature.