Intrinsic compressibility and volume compression in solvated proteins by molecular dynamics simulation at high pressure.

Constant pressure and temperature molecular dynamics techniques have been employed to investigate the changes in structure and volumes of two globular proteins, superoxide dismutase and lysozyme, under pressure. Compression (the relative changes in the proteins' volumes), computed with the Voronoi technique, is closely related with the so-called protein intrinsic compressibility, estimated by sound velocity measurements. In particular, compression computed with Voronoi volumes predicts, in agreement with experimental estimates, a negative bound water contribution to the apparent protein compression. While the use of van der Waals and molecular volumes underestimates the intrinsic compressibilities of proteins, Voronoi volumes produce results closer to experimental estimates. Remarkably, for two globular proteins of very different secondary structures, we compute identical (within statistical error) protein intrinsic compressions, as predicted by recent experimental studies. Changes in the protein interatomic distances under compression are also investigated. It is found that, on average, short distances compress less than longer ones. This nonuniform contraction underlines the peculiar nature of the structural changes due to pressure in contrast with temperature effects, which instead produce spatially uniform changes in proteins. The structural effects observed in the simulations at high pressure can explain protein compressibility measurements carried out by fluorimetric and hole burning techniques. Finally, the calculation of the proteins static structure factor shows significant shifts in the peaks at short wavenumber as pressure changes. These effects might provide an alternative way to obtain information concerning compressibilities of selected protein regions.

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.

Direct measurement of hydrogen bonding in DNA nucleotide bases by atomic force microscopy.

We have used self-assembled purines and pyrimidines on planar gold surfaces and on gold-coated atomic force microscope (AFM) tips to directly probe intermolecular hydrogen bonds. Electron spectroscopy for chemical analysis (ESCA) and thermal programmed desorption (TPD) measurements of the molecular layers suggested monolayer coverage and a desorption energy of about 25 kcal/mol. Experiments were performed under water, with all four DNA bases immobilized on AFM tips and flat surfaces. Directional hydrogen-bonding interaction between the tip molecules and the surface molecules could be measured only when opposite base-pair coatings were used. The directional interactions were inhibited by excess nucleotide base in solution. Nondirectional van der Waals forces were present in all other cases. Forces as low as two interacting base pairs have been measured. With coated AFM tips, surface chemistry-sensitive recognition atomic force microscopy can be performed.

Análise Conformacional e estudo das interações eletrônicas de algumas α-fenilseleno-α-dietóxifosforilacetofenonas para-substituídas

A presente Dissertação relata a síntese e o estudo conformacional das α-fenilseleno-α-dietóxifosforilacetofenonas para-substituídas p-X-Φ-C(O)CH[SeΦ][P(O)(OEt2] (X=OMe 1, Me 2, H 3, F 4, Cl 5, Br 6 e NO2 7) através da banda de estiramento da carbonila no infravermelho, em solventes de polaridade crescente apoiado por cálculos ab initio HF/6-3IG**. A comparação entre a freqüência e a intensidade relativa dos componentes do dubleto, para os derivados 6 e 7, e do singleto para os derivados 1-5, no solvente apolar tetracloreto de carbono, e dos componentes do dubleto, nos solventes de polaridade crescente (clorofórmio, diclorometano e acetonitrila), para os derivados 1-7, com os dados do cálculo ab initio de 3 (composto de referência), indicou que ambas as conformações estáveis (g1 e g2) apresentam a ligação C-Se na geometria anti-clinal (gauche) em relação à carbonila (C=O), enquanto que a ligação C-P assume uma geometria sin-periplanar (cis) em relação à carbonila. A análise dos contatos interatômicos de átomos relevante em comparação com a soma de seus raios de van der Waals, indicou que ambas as conformações g1 e g2 são fortemente estabilizadas pelo sinergismo das interações orbitalares e eletrostáticas π*(CO) / nSe e Oδ-[CO].....Pδ+[PO]. Analogamente, as interações mais fracas Oδ-[OR]..... Cδ+[CO], 0-Hδ+[SeΦ]....Oδ-[PO] e 0-Hδ+[ΦC(O)]....Oδ-[CO] estabilizam as conformações g1 e g2, aproximadamente na mesma extensão. No entanto, somente a conformação g1 é estabilizada pela interação eletrostática (ligação de hidrogênio) Hδ+[α-CH].....Oδ-[OR], enquanto que sómente a conformação g2 é desestabilizada pelo Efeito de Campo Repulsivo entre os dipolos Cδ+=.Oδ- e Pδ+-ORδ- Assim sendo, pode-se concluir que no dubleto de VCO no IV, o componente de maior freqüência e de menor intensidade corresponde à conformação menos estável g2 (do cálculo) enquanto que o componente de menor freqüência e mais intenso corresponde à conformação mais estável g1 (do cálculo). Estes dados estão de pleno acordo com os deslocamentos de freqüência mais negativos da carbonila (ΔVCO) do confôrmero mais estável g1 em relação ao menos estável g2.

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.

Electrostatically mediated specific adsorption of small molecules in metallo-organic frameworks

We investigate the interaction of ethylene and ethane with a Cu-tricarboxylate complex and show that at low loadings the lighter molecule has a higher binding energy as a result of an increased interaction with the framework Cu and stronger hydrogen bonding with the basic framework oxygens. This leads to selective adsorption of ethylene by a factor of about 2 at low pressure, which is overcome by the stronger van der Waals interaction of ethane at high loadings, explaining recent literature data. The results suggest the possibility of separation of light hydrocarbons at low pressures or in trace amounts.

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.