Molecular modelling studies of binding of DACD derivatives into G-Quadruplex DNA: comparison of force field and quantum polarized ligand docking methods

The DNA G-qadruplexes are one of the targets being actively explored for anti-cancer therapy by inhibiting them through small molecules. This computational study was conducted to predict the binding strengths and orientations of a set of novel dimethyl-amino-ethyl-acridine (DACA) analogues that are designed and synthesized in our laboratory, but did not diffract in Synchrotron light.Thecrystal structure of DNA G-Quadruplex(TGGGGT)4(PDB: 1O0K) was used as target for their binding properties in our studies.We used both the force field (FF) and QM/MM derived atomic charge schemes simultaneously for comparing the predictions of drug binding modes and their energetics. This study evaluates the comparative performance of fixed point charge based Glide XP docking and the quantum polarized ligand docking schemes. These results will provide insights on the effects of including or ignoring the drug-receptor interfacial polarization events in molecular docking simulations, which in turn, will aid the rational selection of computational methods at different levels of theory in future drug design programs. Plenty of molecular modelling tools and methods currently exist for modelling drug-receptor or protein-protein, or DNA-protein interactionssat different levels of complexities.Yet, the capasity of such tools to describevarious physico-chemical propertiesmore accuratelyis the next step ahead in currentresearch.Especially, the usage of most accurate methods in quantum mechanics(QM) is severely restricted by theirtedious nature. Though the usage of massively parallel super computing environments resulted in a tremendous improvement in molecular mechanics (MM) calculations like molecular dynamics,they are still capable of dealing with only a couple of tens to hundreds of atoms for QM methods. One such efficient strategy that utilizes thepowers of both MM and QM are the QM/MM hybrid methods. Lately, attempts have been directed towards the goal of deploying several different QM methods for betterment of force field based simulations, but with practical restrictions in place. One of such methods utilizes the inclusion of charge polarization events at the drug-receptor interface, that is not explicitly present in the MM FF.

Viscous coupling of shear-free turbulence across nearly flat fluid interfaces

The interactions between shear-free turbulence in two regions (denoted as + and − on either side of a nearly flat horizontal interface are shown here to be controlled by several mechanisms, which depend on the magnitudes of the ratios of the densities, ρ+/ρ−, and kinematic viscosities of the fluids, μ+/μ−, and the root mean square (r.m.s.) velocities of the turbulence, u0+/u0−, above and below the interface. This study focuses on gas–liquid interfaces so that ρ+/ρ− ≪ 1 and also on where turbulence is generated either above or below the interface so that u0+/u0− is either very large or very small. It is assumed that vertical buoyancy forces across the interface are much larger than internal forces so that the interface is nearly flat, and coupling between turbulence on either side of the interface is determined by viscous stresses. A formal linearized rapid-distortion analysis with viscous effects is developed by extending the previous study by Hunt & Graham (J. Fluid Mech., vol. 84, 1978, pp. 209–235) of shear-free turbulence near rigid plane boundaries. The physical processes accounted for in our model include both the blocking effect of the interface on normal components of the turbulence and the viscous coupling of the horizontal field across thin interfacial viscous boundary layers. The horizontal divergence in the perturbation velocity field in the viscous layer drives weak inviscid irrotational velocity fluctuations outside the viscous boundary layers in a mechanism analogous to Ekman pumping. The analysis shows the following. (i) The blocking effects are similar to those near rigid boundaries on each side of the interface, but through the action of the thin viscous layers above and below the interface, the horizontal and vertical velocity components differ from those near a rigid surface and are correlated or anti-correlated respectively. (ii) Because of the growth of the viscous layers on either side of the interface, the ratio uI/u0, where uI is the r.m.s. of the interfacial velocity fluctuations and u0 the r.m.s. of the homogeneous turbulence far from the interface, does not vary with time. If the turbulence is driven in the lower layer with ρ+/ρ− ≪ 1 and u0+/u0− ≪ 1, then uI/u0− ~ 1 when Re (=u0−L−/ν−) ≫ 1 and R = (ρ−/ρ+)(v−/v+)1/2 ≫ 1. If the turbulence is driven in the upper layer with ρ+/ρ− ≪ 1 and u0+/u0− ≫ 1, then uI/u0+ ~ 1/(1 + R). (iii) Nonlinear effects become significant over periods greater than Lagrangian time scales. When turbulence is generated in the lower layer, and the Reynolds number is high enough, motions in the upper viscous layer are turbulent. The horizontal vorticity tends to decrease, and the vertical vorticity of the eddies dominates their asymptotic structure. When turbulence is generated in the upper layer, and the Reynolds number is less than about 106–107, the fluctuations in the viscous layer do not become turbulent. Nonlinear processes at the interface increase the ratio uI/u0+ for sheared or shear-free turbulence in the gas above its linear value of uI/u0+ ~ 1/(1 + R) to (ρ+/ρ−)1/2 ~ 1/30 for air–water interfaces. This estimate agrees with the direct numerical simulation results from Lombardi, De Angelis & Bannerjee (Phys. Fluids, vol. 8, no. 6, 1996, pp. 1643–1665). Because the linear viscous–inertial coupling mechanism is still significant, the eddy motions on either side of the interface have a similar horizontal structure, although their vertical structure differs.

The role of protein hydrophobicity in thionin–phospholipid interactions: a comparison of α1 and α2-purothionin adsorbed anionic phospholipid monolayers

The plant defence proteins α1- and α2-purothionin (Pth) are type 1 thionins from common wheat (Triticum aestivum). These highly homologous proteins possess characteristics common amongst antimicrobial peptides and proteins, that is, cationic charge, amphiphilicity and hydrophobicity. Both α1- and α2-Pth possess the same net charge, but differ in relative hydrophobicity as determined by C18 reversed phase HPLC. Brewster angle microscopy, X-ray and neutron reflectometry, external reflection FTIR and associated surface pressure measurements demonstrated that α1 and α2-Pth interact strongly with condensed phase 1,2-dipalmitoyl-sn-glycero-3-phospho-(1'-rac-glycerol) (DPPG) monolayers at the air/liquid interface. Both thionins disrupted the in-plane structure of the anionic phospholipid monolayer, removing lipid during this process and both penetrated the lipid monolayer in addition to adsorbing as a single protein layer to the lipid head-group. However, analysis of the interfacial structures revealed that the α2-Pth showed faster disruption of the lipid film and removed more phospholipid (12%) from the interface than α1-Pth. Correlating the protein properties and lipid binding activity suggests that hydrophobicity plays a key role in the membrane lipid removal activity of thionins.

Effect of shear rupture on aggregate scale formation in sea ice

A discrete element model is used to study shear rupture of sea ice under convergent wind stresses. The model includes compressive, tensile, and shear rupture of viscous elastic joints connecting floes that move under the action of the wind stresses. The adopted shear rupture is governed by Coulomb’s criterion. The ice pack is a 400 km long square domain consisting of 4 km size floes. In the standard case with tensile strength 10 times smaller than the compressive strength, under uniaxial compression the failure regime is mainly shear rupture with the most probable scenario corresponding to that with the minimum failure work. The orientation of cracks delineating formed aggregates is bimodal with the peaks around the angles given by the wing crack theory determining diamond-shaped blocks. The ice block (floe aggregate) size decreases as the wind stress gradient increases since the elastic strain energy grows faster leading to a higher speed of crack propagation. As the tensile strength grows, shear rupture becomes harder to attain and compressive failure becomes equally important leading to elongation of blocks perpendicular to the compression direction and the blocks grow larger. In the standard case, as the wind stress confinement ratio increases the failure mode changes at a confinement ratio within 0.2–0.4, which corresponds to the analytical critical confinement ratio of 0.32. Below this value, the cracks are bimodal delineating diamond shape aggregates, while above this value failure becomes isotropic and is determined by small-scale stress anomalies due to irregularities in floe shape.

Relaxation of surface tension in the free-surface boundary layers of simple Lennard-Jones liquids

In this paper we use molecular dynamics to answer a classical question: how does the surface tension on a liquid/gas interface appear? After defining surface tension from the first principles and performing several consistency checks, we perform a dynamic experiment with a single simple liquid nanodroplet. At time zero, we remove all molecules of the interfacial layer of molecules, creating a fresh bare interface with the bulk arrangement of molecules. After that the system evolves towards equilibrium, and the expected surface tension is re-established. We found that the system relaxation consists of three distinct stages. First, the mechanical balance is quickly re-established. During this process the notion of surface tension is meaningless. In the second stage, the surface tension equilibrates, and the density profile broadens to a value which we call “intrinsic” interfacial width. During the third stage, the density profile continues to broaden due to capillary wave excitations, which does not however affect the surface tension.We have observed this scenario for monatomic Lennard-Jones (LJ) liquid as well as for binary LJ mixtures at different temperatures, monitoring a wide range of physical observables.

Potential reduction of concrete deterioration through controlled DEF in hydrated concrete

Delayed ettringite formation (DEF) is a chemical reaction with proven damaging effects on hydrated concrete. Ettringite crystals can cause cracks and their widening due to pressure on cracked walls caused by the positive volume difference in the reaction. Concrete may show improvements in strength at early ages but further growth of cracks causes widening and spreading through the concrete structure. In this study, finely dispersed crystallization nuclei achieved by adding air-entraining agent (AEA) and short vibration of specimens is presented as the main prerequisite for reducing DEF-induced deterioration of hydrated concrete. The study presents the method and mechanism for obtaining the required nucleation. Controlling long-term DEF by providing AEA-induced crystallisation nuclei, prevented excessive and rapid initial strength improvements, and resulted in a slight increase of compressive strength of fine grained concrete with only marginally lower density.

Electrochemical sensing of 2D condensation in amyloid peptides

The interfacial behavior of the model amyloid peptide octamer YYKLVFFC (peptide 1) and two other amyloid peptides YEVHHQKLVFF (peptide 2) and KKLVFFA (peptide 3) at the metal|aqueous solution interface was studied by voltammetric and constant current chronopotentiometric stripping (CPS). All three peptides are adsorbed in a wide potential range and exhibit different interfacial organizations depending on the electrode potential. At the least negative potentials, chemisorption of peptide 1 occurs through the formation of a metal sulfur bond. This bond is broken close to −0.6 V. The peptide undergoes self-association at more negative potentials, leading to the formation of a “pit” characteristic of a 2D condensed film. Under the same conditions the other peptides do not produce such a pit. Formation of the 2D condensed layer in peptide 1 is supported by the time, potential and temperature dependences of the interfacial capacity and it is shown that presence of the 2D layer is reflected by the peptide CPS signals due to the catalytic hydrogen evolution. The ability of peptide 1 to form the potential-dependent 2D condensed layer has been reported neither for any other peptide nor for any protein molecule. This ability might be related to the well-known oligomerization and aggregation of Alzheimer amyloid peptides.

Modelling the reorientation of sea-ice faults as the wind changes direction

A discrete-element model of sea ice is used to study how a 90° change in wind direction alters the pattern of faults generated through mechanical failure of the ice. The sea-ice domain is 400km in size and consists of polygonal floes obtained through a Voronoi tessellation. Initially the floes are frozen together through viscous–elastic joints that can break under sufficient compressive, tensile and shear deformation. A constant wind-stress gradient is applied until the initially frozen ice pack is broken into roughly diamond-shaped aggregates, with crack angles determined by wing-crack formation. Then partial refreezing of the cracks delineating the aggregates is modelled through reduction of their length by a particular fraction, the ice pack deformation is neglected and the wind stress is rotated by 90°. New cracks form, delineating aggregates with a different orientation. Our results show the new crack orientation depends on the refrozen fraction of the initial faults: as this fraction increases, the new cracks gradually rotate to the new wind direction, reaching 90° for fully refrozen faults. Such reorientation is determined by a competition between new cracks forming at a preferential angle determined by the wing-crack theory and at old cracks oriented at a less favourable angle but having higher stresses due to shorter contacts across the joints

Predicting the orientation of lipid cubic phase films

Lipid cubic phase films are of increasingly widespread importance, both in the analysis of the cubic phases themselves by techniques including microscopy and X-ray scattering, and in their applications, especially as electrode coatings for electrochemical sensors and for templates for the electrodeposition of nanostructured metal. In this work we demonstrate that the crystallographic orientation adopted by these films is governed by minimization of interfacial energy. This is shown by the agreement between experimental data obtained using grazing-incidence small-angle X-ray scattering (GI-SAXS), and the predicted lowest energy orientation determined using a theoretical approach we have recently developed. GI-SAXS data show a high degree of orientation for films of both the double diamond phase and the gyroid phase, with the [111] and [110] directions respectively perpendicular to the planar substrate. In each case, this matches the lowest energy facet calculated for that particular phase.

Electronic structure of Pd multi-layers on Re(0001): the role of charge transfer

Understanding the origin of the properties of metal-supported metal thin films is important for the rational design of bimetallic catalysts and other applications, but it is generally difficult to separate effects related to strain from those arising from interface interactions. Here we use density functional (DFT) theory to examine the structure and electronic behavior of few-layer palladium films on the rhenium (0001) surface, where there is negligible interfacial strain and therefore other effects can be isolated. Our DFT calculations predict stacking sequences and interlayer separations in excellent agreement with quantitative low-energy electron diffraction experiments. By theoretically simulating the Pd core-level X-ray photoemission spectra (XPS) of the films, we are able to interpret and assign the basic features of both low-resolution and high-resolution XPS measurements. The core levels at the interface shift to more negative energies, rigidly following the shifts in the same direction of the valence d-band center. We demonstrate that the valence band shift at the interface is caused by charge transfer from Re to Pd, which occurs mainly to valence states of hybridized s-p character rather than to the Pd d-band. Since the d-band filling is roughly constant, there is a correlation between the d-band center shift and its bandwidth. The resulting effect of this charge transfer on the valence d-band is thus analogous to the application of a lateral compressive strain on the adlayers. Our analysis suggests that charge transfer should be considered when describing the origin of core and valence band shifts in other metal / metal adlayer systems.

Influence of energy and pulse repetition rate of Er : YAG laser on enamel ablation ability and morphological analysis of the laser-irradiated surface

Objective: To assess the influence of energy and pulse repetition rate of Er:YAG laser on the enamel ablation ability and substrate morphology. Methods: Fifteen crowns of molars were sectioned in four fragments, providing 60 samples, which were ground to flatten the enamel surface. The initial mass was obtained by weighing the fragments. The specimens were hydrated for I h, fixed, and a 3-mm-diameter area was delimited. Twelve groups were randomly formed according to the combination of laser energies (200, 250, 300, or 350 mJ) and pulse repetition rates (2, 3, or 4 Hz). The final mass was obtained and mass loss was calculated by the difference between the initial and final mass. The specimens were prepared for SEM. Data were submitted to ANOVA and Scheffe test. Results: The 4 Hz frequency resulted in higher mass loss and was statistically different from 2 and 3 Hz (p < 0.05). The increase of frequency produced more melted areas, cracks, and unselective and deeper ablation. The 350 mJ energy promoted greater mass loss, similar to 300 mJ. Conclusions: The pulse repetition rate influenced more intensively the mass loss and morphological alteration. Among the tested parameters, 350 mJ/3 Hz improved the ability of enamel ablation with less surface morphological alterations. (C) 2007 Wiley Periodicals, Inc. J Biomed Mater Res.

Effect of extrusion on the emulsifying properties of soybean proteins and pectin mixtures modelled by response surface methodology

The objective of this study was to apply response surface methodology to estimate the emulsifying capacity and stability of mixtures containing isolated and textured soybean proteins combined with pectin and to evaluate if the extrusion process affects these interfacial properties. A simplex-centroid design was applied to the model emulsifying activity index (EAI), average droplet size (D-[4.3]) and creaming inhibition (Cl%) of the mixtures. All models were significant and able to explain more than 86% of the variation. The high predictive capacity of the models was also confirmed. The mean values for EAI, D-[4.3] and Cl% observed in all assays were 0.173 +/- 0.015 mn, 19.2 +/- 1.0 mu m and 53.3 +/- 2.6%, respectively. No synergism was observed between the three compounds. This result can be attributed to the low soybean protein solubility at pH 6.2 (<35%). Pectin was the most important variable for improving all responses. The emulsifying capacity of the mixture increased 41% after extrusion. Our results showed that pectin could substitute or improve the emulsifying properties of the soybean proteins and that the extrusion brings additional advantage to interfacial properties of this combination. (C) 2008 Elsevier Ltd. All rights reserved.

Selection of operational parameters for the production of instant soy protein isolate by pulsed fluid bed agglomeration

The objective of this study was to select the optimal operational conditions for the production of instant soy protein isolate (SPI) by pulsed fluid bed agglomeration. The spray-dried SPI was characterized as being a cohesive powder, presenting cracks and channeling formation during its fluidization (Geldart type A). The process was carried out in a pulsed fluid bed, and aqueous maltodextrin solution was used as liquid binder. Air pulsation, at a frequency of 600 rpm, was used to fluidize the cohesive SPI particles and to allow agglomeration to occur. Seventeen tests were performed according to a central composite design. Independent variables were (i) feed flow rate (0.5-3.5 g/min), (ii) atomizing air pressure (0.5-1.5 bar) and (iii) binder concentration (10-50%). Mean particle diameter, process yield and product moisture were analyzed as responses. Surface response analysis led to the selection of optimal operational parameters, following which larger granules with low moisture content and high process yield were produced. Product transformations were also evaluated by the analysis of size distribution, flowability, cohesiveness and wettability. When compared to raw material, agglomerated particles were more porous and had a more irregular shape, presenting a wetting time decrease, free-flow improvement and cohesiveness reduction. (C) 2010 Elsevier B.V. All rights reserved.

Surface properties of eucalyptus pulp fibres as reinforcement of cement-based composites

The objective of the present work is to evaluate the effects of the surface properties of unrefined eucalyptus pulp fibres concerning their performance in cement-based composites. The influence of the fibre surface on the microstructure of fibre-cement composites was evaluated after accelerated ageing cycles, which simulate natural weathering. The surface of unbleached pulp is a thin layer that is rich in cellulose, lignin, hemicelluloses, and extractives. Such a layer acts as a physical and chemical barrier to the penetration of low molecular components of cement. The unbleached fibres are less hydrophilic than the bleached ones. Bleaching removes the amorphous lignin and extractives from the surface and renders it more permeable to liquids. Atomic force microscopy (AFM) helps in understanding the fibre-cement interface. Bleaching improved the fibre- cement interfacial bonding, whereas fibres in the unbleached pulp were less susceptible to the re-precipitation of cement hydration products into the fibre cavities (lumens). Therefore, unbleached fibres can improve the long-term performance of the fibre-cement composite owing to their delayed mineralization.

Molecular Ordering of Layer-by-Layer Polyelectrolyte Films Studied by Sum-Frequency Vibrational Spectroscopy

The molecular arrangement in organic thin films is crucial for their increasing technological applications. Here, we use vibrational spectroscopy by sum-frequency generation (SFG) to study the ordering of polyelectrolyte layers adsorbed on silica for all steps of layer-by-layer (LbL) self-assembly. In situ measurements during adsorption and rinsing showed that the adsorbed polymer has a disordered conformation and confirmed surface charge overcompensation upon polyelectrolyte adsorption by probing the interfacial electric field. In dry films, the polymer chains acquired a net orientational ordering, which was affected, however, by the adsorption of subsequent layers. Such a detailed characterization may allow the control of LbL film structure and functionality with unprecedented power.