A new numerical approach to dense polymer brushes and surface instabilities.

We present a numerical self-consistent field (SCF) method which describes freely jointed chains of spherical monomers applied to densely grafted polymer brushes. We discuss both the Flory-Huggins model and the Carnahan-Starling equation of state and show the latter being preferable within our model at polymer volume fractions above 10%. We compare the results of our numerical method with data from molecular dynamics (MD) simulations [G.-L. He, H. Merlitz, J.-U. Sommer, and C.-X. Wu, Macromolecules 40, 6721 (2007)] and analytical SCF calculations [P. M. Biesheuvel, W. M. de Vos, and V. M. Amoskov, Macromolecules 41, 6254 (2008)] and obtain close agreement between the density profiles up to high grafting densities. In contrast to prior numerical and analytical studies of densely grafted polymer brushes our method provides detailed information about chain configurations including fluctuation, depletion, and packing effects. Using our model we could study the recently discovered instability of densely grafted polymer brushes with respect to slight variations of individual chain lengths, driven by fluctuation effects [H. Merlitz, G.-L. He, C.-X. Wu, and J.-U. Sommer, Macromolecules 41, 5070 (2008)]. The obtained results are in very close agreement with corresponding MD simulations.

Receptor flexibility for large-scale in silico ligand screens: chances and challenges.

An important contribution to today's computer-aided drug design is the automated screening of large compound databases against structurally resolved protein receptors targets. The introduction of ligand flexibility has, by now, become a standardized procedure. In contrast, a general approach to treat target degrees of freedom is still to be found, a consequence of the extreme increase of computational complexity, which comes along with the relaxation of protein degrees of freedom. In this chapter, we discuss in some detail both benefits and present limitations of target flexibility for high-throughput in silico database screens. Among the benefits are an improved diversity of binding modes, which allows one to identify a wider class of drug candidates. The limitations are related to a diminishing docking accuracy and an increased number of false hits. Using the thymidine kinase receptor and ten known inhibitors as an example, we describe in detail how target flexibility was implemented and how it affected the screening performance.

Polymer brushes near the crystallization density.

In this paper, polymer brushes are studied via molecular-dynamics simulations at very high grafting densities, where the crossover between the brush regime and the polymer-crystal regime is taking place. This crossover is directly observed with the structure factor and pair-correlation function. With increasing grafting density, this crystallization is progressing from the core layer of the brush towards the surface layer. The same process is analyzed using the lateral fluctuations of the monomers as a signature of their diminishing mobility. Additionally, bond forces and the chain excess free energy indicate a transition from the brush regime to the overstretched regime, which is in agreement with predictions of a modified self-consistent field theory.

Fluctuation analysis and accuracy of a large-scale in silico screen.

Using a cascadic version of the stochastic tunneling method we perform an all-atom database screen over 186,000 flexible ligands of the NCI 3D database against the thymidine kinase receptor. By analyzing the errors in the binding energy we demonstrate how the cascadic technique is superior to conventional sequential docking techniques and how reliable results for the determination of the top-scoring ligands could be achieved. The substrate corresponding to the crystal structure used in the screen ranks in the upper 0.05% of the database, validating both docking methodology and the applicability of the scoring function to this substrate. Several high ranking ligands of the database display significant structural similarity with known substrates. A detailed analysis of the accuracy of the screening method is carried out, and its dependence on the flexibility of the ligand is quantified.

A Brownian dynamics program for the simulation of linear and circular DNA and other wormlike chain polyelectrolytes.

For the interpretation of solution structural and dynamic data of linear and circular DNA molecules in the kb range, and for the prediction of the effect of local structural changes on the global conformation of such DNAs, we have developed an efficient and easy way to set up a program based on a second-order explicit Brownian dynamics algorithm. The DNA is modeled by a chain of rigid segments interacting through harmonic spring potentials for bending, torsion, and stretching. The electrostatics are handled using precalculated energy tables for the interactions between DNA segments as a function of relative orientation and distance. Hydrodynamic interactions are treated using the Rotne-Prager tensor. While maintaining acceptable precision, the simulation can be accelerated by recalculating this tensor only once in a certain number of steps.

Looping dynamics of linear DNA molecules and the effect of DNA curvature: a study by Brownian dynamics simulation.

A Brownian dynamics (BD) model described in the accompanying paper (Klenin, K., H. Merlitz, and J. Langowski. 1998. A Brownian dynamics program for the simulation of linear and circular DNA, and other wormlike chain polyelectrolytes. Biophys. J. 74:000-000) has been used for computing the end-to-end distance distribution function, the cyclization probability, and the cyclization kinetics of linear DNA fragments between 120 and 470 basepairs with optional insertion of DNA bends. Protein-mediated DNA loop formation was modeled by varying the reaction distance for cyclization between 0 and 10 nm. The low cyclization probability of DNA fragments shorter than the Kuhn length (300 bp) is enhanced by several orders of magnitude when the cyclization is mediated by a protein bridge of 10 nm diameter, and/or when the DNA is bent. From the BD trajectories, end-to-end collision frequencies were computed. Typical rates for loop formation of linear DNAs are 1.3 x 10(3) s(-1) (235 bp) and 4.8 x 10(2) s(-1) (470 bp), while the insertion of a 120 degree bend in the center increases this rate to 3.0 x 10(4) s(-1) (235 bp) and 5.5 x 10(3) s(-1) (470 bp), respectively. The duration of each encounter is between 0.05 and 0.5 micros for these DNAs. The results are discussed in the context of the interaction of transcription activator proteins.

Salt effects on the structure and internal dynamics of superhelical DNAs studied by light scattering and Brownian dynamics.

Using laser light scattering, we have measured the static and dynamic structure factor of two different superhelical DNAs, p1868 (1868 bp) and simian virus 40 (SV40) (5243 bp), in dilute aqueous solution at salt concentrations between 1 mM and 3 M NaCl. For both DNA molecules, Brownian dynamics (BD) simulations were also performed, using a previously described model. A Fourier mode decomposition procedure was used to compute theoretical light scattering autocorrelation functions (ACFs) from the BD trajectories. Both measured and computed autocorrelation functions were then subjected to the same multiexponential decomposition procedure. Simulated and measured relaxation times as a function of scattering angle were in very good agreement. Similarly, computed and measured static structure factors and radii of gyration agreed within experimental error. One main result of this study is that the amplitudes of the fast-relaxing component in the ACF show a peak at 1 M salt concentration. This nonmonotonic behavior might be caused by an initial increase in the amplitudes of internal motions due to diminishing long-range electrostatic repulsions, followed by a decrease at higher salt concentration due to a compaction of the structure.