Modified AutoDock for accurate docking of protein kinase inhibitors

Protein kinases are an important class of enzymes controlling virtually all cellular signaling pathways. Consequently, selective inhibitors of protein kinases have attracted significant interest as potential new drugs for many diseases. Computational methods, including molecular docking, have increasingly been used in the inhibitor design process [1]. We have considered several docking packages in order to strengthen our kinase inhibitor work with computational capabilities. In our experience, AutoDock offered a reasonable combination of accuracy and speed, as opposed to methods that specialize either in fast database searches or detailed and computationally intensive calculations.However, AutoDock did not perform well in cases where extensive hydrophobic contacts were involved, such as docking of SB203580 to its target protein kinase p38. Another shortcoming was a hydrogen bonding energy function, which underestimated the attraction component and, thus, did not allow for sufficiently accurate modeling of the key hydrogen bonds in the kinase-inhibitor complexes.We have modified the parameter set used to model hydrogen bonds, which increased the accuracy of AutoDock and appeared to be generally applicable to many kinase-inhibitor pairs without customization. Binding to largely hydrophobic sites, such as the active site of p38, was significantly improved by introducing a correction factor selectively affecting only carbon and hydrogen energy grids, thus, providing an effective, although approximate, treatment of solvation.     

Effect of Coriolis force on counter-current chromatographic separation by centrifugal partition chromatography

The effect of Coriolis force on the counter-current chromatographic separation was studied using centrifugal partition chromatography (CPC) with four different two-phase solvent systems including n-hexane-acetonitrile (ACN); tert-butyl methyl ether (MtBE)-aqueous 0.1% trifluoroacetic acid (TFA) (1:1); MtBE-ACN-aqueous 0.1% TFA (2:2:3); and 12.5% (w/w) polyethylene glycol (PEG) 1000-12.5% (w/w) dibasic potassium phosphate. Each separation was performed by eluting either the upper phase in the ascending mode or the lower phase in the descending mode, each in clockwise (CW) and counterclockwise column rotation. Better partition efficiencies were attained by the CW rotation in both mobile phases in all the two-phase solvent systems examined. The mathematical analysis also revealed the Coriolis force works favorably under the CW column rotation for both mobile phases. The overall results demonstrated that the Coriolis force produces substantial effects on CPC separation in both organic-aqueous and aqueous-aqueous two-phase systems.     

Simple test system for single molecule recognition force microscopy

We have established an easy-to-use test system for detecting receptor-ligand interactions on the single molecule level using atomic force microscopy (AFM). For this, avidin-biotin, probably the best characterized receptor-ligand pair, was chosen. AFM sensors were prepared containing tethered biotin molecules at sufficiently low surface concentrations appropriate for single molecule studies. A biotin tether, consisting of a 6 nm poly(ethylene glycol) (PEG) chain and a functional succinimide group at the other end, was newly synthesized and covalently coupled to amine-functionalized AFM tips. In particular, PEG₈₀₀ diamine was glutarylated, the mono-adduct NH₂-PEG-COOH was isolated by ion exchange chromatography and reacted with biotin succinimidylester to give biotin-PEG-COOH which was then activated as N-hydroxysuccinimide (NHS) ester to give the biotin-PEG-NHS conjugate which was coupled to the aminofunctionalized AFM tip. The motional freedom provided by PEG allows for free rotation of the biotin molecule on the AFM sensor and for specific binding to avidin which had been adsorbed to mica surfaces via electrostatic interactions. Specific avidin-biotin recognition events were discriminated from nonspecific tip-mica adhesion by their typical unbinding force (∼40 pN at 1.4 nN/s loading rate), unbinding length (<13 nm), the characteristic nonlinear force-distance relation of the PEG linker, and by specific block with excess of free d-biotin. The convenience of the test system allowed to evaluate, and compare, different methods and conditions of tip aminofunctionalization with respect to specific binding and nonspecific adhesion. It is concluded that this system is well suited as calibration or start-up kit for single molecule recognition force microscopy.     

Comparative characterisation by atomic force microscopy and ellipsometry of soft and solid thin films

Ellipsometry and atomic force microscopy (AFM) were used to study the film thickness and the surface roughness of both ‘soft’ and solid thin films. ‘Soft’ polymer thin films of polystyrene and poly(styrene–ethylene/butylene–styrene) block copolymer were prepared by spin‐coating onto planar silicon wafers. Ellipsometric parameters were fitted by the Cauchy approach using a two‐layer model with planar boundaries between the layers. The smooth surfaces of the prepared polymer films were confirmed by AFM. There is good agreement between AFM and ellipsometry in the 80–130 nm thickness range. Semiconductor surfaces (Si) obtained by anisotropic chemical etching were investigated as an example of a randomly rough surface. To define roughness parameters by ellipsometry, the top rough layers were treated as thin films according to the Bruggeman effective medium approximation (BEMA). Surface roughness values measured by AFM and ellipsometry show the same tendency of increasing roughness with increased etching time, although AFM results depend on the used window size. The combined use of both methods appears to offer the most comprehensive route to quantitative surface roughness characterisation of solid films. Copyright © 2007 John Wiley & Sons, Ltd.     

Viscoelastic properties of single poly(ethylene glycol) molecules.

The viscoelastic properties of single poly(ethylene glycol) (PEG) molecules were measured by analysis of thermally and magnetically driven oscillations of an atomic force microscope (AFM) cantilever/molecule system. The molecular and monomer stiffness and friction of the PEG polymer were derived using a simple harmonic oscillator (SHO) model. Excellent agreement between the values of these two parameters obtained by the two approaches indicates the validity of the SHO model under the experimental regimes and the excellent reproducibility of the techniques. A sharp minimum in the monomeric friction is seen at around 180 pN applied force which we propose is due to a force induced change in the shape of the energy landscape describing the conformational transition of PEG from a helical to a planar state, which in turn affects the timescale of the transition and therefore modifies the measured internal friction. A knowledge of the viscoelastic response of PEG monomers is particularly important since PEG is widely used as a linker molecule for tethering groups of interest to the AFM tip in force spectroscopy experiments, and we show here that care must be exercised because of the force-dependent viscoelastic properties of these linkers.     

A new GROMOS force field for hexopyranose-based carbohydrates

A new parameter set (referred to as 45A4) is developed for the explicit-solvent simulation of hexopyranose-based carbohydrates. This set is compatible with the most recent version of the GROMOS force field for proteins, nucleic acids, and lipids, and the SPC water model. The parametrization procedure relies on: (1) reassigning the atomic partial charges based on a fit to the quantum-mechanical electrostatic potential around a trisaccharide; (2) refining the torsional potential parameters associated with the rotations of the hydroxymethyl, hydroxyl, and anomeric alkoxy groups by fitting to corresponding quantum-mechanical profiles for hexopyranosides; (3) adapting the torsional potential parameters determining the ring conformation so as to stabilize the (experimentally predominant) (4)C(1) chair conformation. The other (van der Waals and nontorsional covalent) parameters and the rules for third and excluded neighbors are taken directly from the most recent version of the GROMOS force field (except for one additional exclusion). The new set is general enough to define parameters for any (unbranched) hexopyranose-based mono-, di-, oligo- or polysaccharide. In the present article, this force field is validated for a limited set of monosaccharides (alpha- and beta-D-glucose, alpha- and beta-D-galactose) and disaccharides (trehalose, maltose, and cellobiose) in solution, by comparing the results of simulations to available experimental data. More extensive validation will be the scope of a forthcoming article. (c) 2005 Wiley Periodicals, Inc. J Comput Chem 26: 1400-1412, 2005.     

Sublimation and Diffusion Kinetics of 2,4,6-Trinitrotoluene (TNT) Single Crystals by Atomic Force Microscopy (AFM)

In this article, we report the in-situ nanoscale experimental measurement of sublimation rates, activation energy of sublimation, and diffusion coefficients of 2,4,6-trinitrotoluene (TNT) single crystals in air using atomic force microscopy (AFM). The crystals were prepared by slow evaporation at 5 °C using acetone-dissolved TNT. The mass loss was calculated by monitoring the shrinkage of the surface area of layered islands formed on the surface of the TNT crystals due to sublimation upon isothermal heating at temperatures below the melting point. The results suggest the sublimation process occurs via two-dimensional detachment of TNT molecules from the non-prominent facets on the crystal surface which imitates the nucleation and crystal growth process. Sublimation rates are one order of magnitude smaller than previously reported values. However, the calculated activation energy (112.15 ± 3.2 kJ/mol) and temperature-dependent sublimation rates agree well with the reported values for TNT thin films and microcrystals determined by UV-vis absorbance spectroscopy and quartz crystal microscopy (QCM) (90–141 kJ/mol). The average diffusion coefficient is (4.35 × 10^–6 m²/s) which is within the range of the reported theoretical values with an average of 5.59 × 10^–6 m²/s, and about 25% less than that determined using thermogravimetric analysis for powder TNT.     

Folding Free Energy Determination of an RNA Three-Way Junction Using Fluctuation Theorems

Nonequilibrium work relations and fluctuation theorems permit us to extract equilibrium information from nonequilibrium measurements. They find application in single-molecule pulling experiments where molecular free energies can be determined from irreversible work measurements by using unidirectional (e.g., Jarzynski’s equality) and bidirectional (e.g., Crooks fluctuation theorem and Bennet’s acceptance ratio (BAR)) methods. However, irreversibility and the finite number of pulls limit their applicability: the higher the dissipation, the larger the number of pulls necessary to estimate ΔG within a few kBT. Here, we revisit pulling experiments on an RNA three-way junction (3WJ) that exhibits significant dissipation and work-distribution long tails upon mechanical unfolding. While bidirectional methods are more predictive, unidirectional methods are strongly biased. We also consider a cyclic protocol that combines the forward and reverse work values to increase the statistics of the measurements. For a fixed total experimental time, faster pulling rates permit us to efficiently sample rare events and reduce the bias, compensating for the increased dissipation. This analysis provides a more stringent test of the fluctuation theorem in the large irreversibility regime.     

双旋转附属圆柱对主圆柱尾流控制的影响

基于Fluent软件平台,采用数值模拟方法对非稳态圆柱体结构尾流流动特性进行了研究。对在Re=50～200范围内,双旋转附属圆柱的转速对主圆柱体尾流流动特性的影响进行了分析。研究结果表明:随着附属圆柱旋转速率的增加,主圆柱体表面所受阻力系数平均值与均方根值、升力系数均方根值均会减小。同时,旋转速率的变化对柱体结构表面压力分布的影响显著,压力系数在附属圆柱的位置产生了跳跃性变化。另外,当附属圆柱转速达到临界值时,尾流涡街变窄,涡脱落现象消失,并且系统的能量效率到达最佳状态。     

Ir(111)表面石墨烯中缺陷的原子结构确认

识别和解析石墨烯中缺陷的精确原子结构是研究不同类型缺陷的物化特性,实现石墨烯物性调控的前提,可以为在原子尺度研究石墨烯缺陷的构效关系提供重要的实验依据.本文结合扫描隧道显微镜(STM)和原子力显微镜(AFM)确认了在Ir(111)表面生长的石墨烯中自发形成的缺陷,以及通过离子轰击方法在石墨烯中引入的多种缺陷结构,包括单空位缺陷、非六元环拓扑结构以及石墨烯层下的基底缺陷.