Theoretical and Experimental Study of Surface Plasmon Radiation Force on Micrometer-Sized Spheres

We report theoretical and experimental study of surface plasmon (SP) radiation force on micrometer-sized polystyrene latex spheres with different radii. Two theoretical methods are introduced to numerically investigate the SP sphere system: the arbitrary beam theory for electromagnetic field distribution and the generalized Lorenz–Mie scattering theory for radiation force calculation. Transport velocity of the spheres in water is tested experimentally and compared with numerical results. Our work also shows effective coupling from SP field to micrometer-sized spheres, with high consistency between experiment and calculation, which has significant potential in future applications.     

The Effect of Electrical Deformation Forces on the Electropermeabilization of Erythrocyte Membranes in Low- and High-Conductivity Media

Electrical breakdown of erythrocytes induces hemoglobin release which increases markedly with decreasing conductivity of the pulse medium. This effect presumably results from the transient, conductivity-dependent deformation forces (elongation or compression) on the cell caused by Maxwell stress. The deformation force is exerted on the plasma membrane of the cell, which can be viewed as a transient dipole induced by an applied DC electric field pulse. The induced dipole arises from the free charges that accumulate at the cell interfaces via the Maxwell-Wagner polarization mechanism. The polarization response of erythrocytes to a DC field pulse was estimated from the experimental data obtained by using two complementary frequency-domain techniques. The response is very rapid, due to the highly conductive cytosol. Measurements of the electrorotation and electrodeformation spectra over a wide conductivity range yielded the information and data required for the calculation of the deformation force as a function of frequency and external conductivity and for the calculation of the transient development of the deformation forces during the application of a DC-field pulse. These calculations showed that (i) electric force precedes and accompanies membrane charging (up to the breakdown voltage) and (ii) that under low-conductivity conditions, the electric stretching force contributes significantly to the enlargement of ``electroleaks' in the plasma membrane generated by electric breakdown. Received: 12 December 1997/Revised: 13 March 1998     

Molecular Dynamics Simulation on the Connection Machine

An algorithm for the efficient calculation of macromolecular force fields on the Connection Machine is described. The full force field is separated into bond interactions and non-bonding interactions. Only the latter are implemented on the Connection Machine, the former, less computationally intensive tasks are performed by an existing, conventional molecular dynamics code on the front end. Parallelization of the evaluation of non-bonding interactions is achieved by the Replicated Systolic Loop algorithm introduced in this paper. The algorithm is a variant of the Systolic Loop scheme often used for the computation of 2-particle forces for the classical N-particle problem.     

A molecular dynamics study of liquid aliphatic alcohols: simulation of density and self-diffusion coefficient using a modified OPLS force field

A set of 13 aliphatic alcohols was modelled by molecular dynamics simulations at temperatures from 288 to 338 K using the optimised potential for liquid simulations (OPLS) united-atom force field, the OPLS all-atom force field and the OPLS all-atom force field with modified partial charges of the hydroxyl group. The set includes primary and secondary alcohols, and mono-, di- and trialcohols, and covers a broad range of polarities from log P = ? 0.74 (methanol) to log P = 2.9 (octanol). The density, the radial distribution function, the self-diffusion coefficient and the dielectric constant were evaluated. A long equilibration time of at least 50 ns and a large size of the molecular system of more than 75,000 atoms were used. Except for glycerol, the OPLS all-atom force field reliably reproduced the experimentally determined density with deviations of less than 4% over the whole temperature range. In contrast, the modelled self-diffusion coefficient deviated from its experimental value by up to 55%. To modify the force field, the partial charges of the hydroxyl group were varied by up to 3%. Using the modified OPLS force field, the deviation of the self-diffusion coefficients from their experimental values decreased to less than 19%, while the densities changed by less than 1%.     

Interchain hydrophobic clustering promotes rigidity in HIV-1 protease flap dynamics: new insights from Molecular Dynamics

The dynamics of HIV-1 protease (HIV-pr), a drug target for HIV infection, has been studied extensively by both computational and experimental methods. The flap dynamics of HIV-pr is considered to be more important for better ligand binding and enzymatic actions. Moreover, it has been demonstrated that the drug-induced mutations can change the flap dynamics of HIV-pr affecting the binding affinity of the ligands. Therefore, detailed understanding of flap dynamics is essential for designing better inhibitors. Previous computational investigations observed significant variation in the flap opening in nanosecond time scale indicating that the dynamics is highly sensitive to the simulation protocols. To understand the sensitivity of the flap dynamics on the force field and simulation protocol, molecular dynamics simulations of HIV-pr have been performed with two different AMBER force fields, ff99 and ff02. Two different trajectories (20?ns each) were obtained using the ff99 and ff02 force field. The results showed polarizable force field (ff02) make the flap tighter than the nonpolarizable force field (ff99). Some polar interactions and hydrogen bonds involving flap residues were found to be stronger with ff02 force field. The formation of interchain hydrophobic cluster (between flap tip of one chain and active site wall of another chain) was found to be dominant in the semi-open structures obtained from the simulations irrespective of the force field. It is proposed that an inhibitor, which will promote this interchain hydrophobic clustering, may make the flaps more rigid, and presumably the effect of mutation would be small on ligand binding.     

Peptoid conformational free energy landscapes from implicit-solvent molecular simulations in AMBER

To test the accuracy of existing AMBER force field models in predicting peptoid conformation and dynamics, we simulated a set of model peptoid molecules recently examined by Butterfoss et al. (JACS 2009, 131, 16798-16807) using QM methods as well as three peptoid sequences with experimentally determined structures. We found that AMBER force fields, when used with a Generalized Born/Surface Area (GBSA) implicit solvation model, could accurately reproduce the peptoid torsional landscape as well as the major conformers of known peptoid structures. Enhanced sampling by replica exchange molecular dynamics (REMD) using temperatures from 300 to 800 K was used to sample over cis-trans isomerization barriers. Compared to (Nrch)5 and cyclo-octasarcosyl, the free energy of N-(2-nitro-3-hydroxyl phenyl)glycine-N-(phenyl)glycine has the most "foldable" free energy landscape, due to deep trans-amide minima dictated by N-aryl sidechains. For peptoids with (S)-N (1-phenylethyl) (Nspe) side chains, we observe a discrepancy in backbone dihedral propensities between molecular simulations and QM calculations, which may be due to force field effects or the inability to capture n --> n* interactions. For these residues, an empirical phi-angle biasing potential can "rescue" the backbone propensities seen in QM. This approach can serve as a general strategy for addressing force fields without resorting to a complete reparameterization. Overall, this study demonstrates the utility of implicit-solvent REMD simulations for efficient sampling to predict peptoid conformational landscapes, providing a potential tool for first-principles design of sequences with specific folding properties.     

Analysis of force fields and BET theory for polymers of intrinsic microporosity

A detailed force field analysis for polymers of intrinsic microporosity (PIMs) was carried out in this study. The generalised amber force field (GAFF) with united atom transferable potential for phase equilibria (TraPPE-UA), and the atomistic polymer consistent force field were evaluated. Analysis carried out with PIM-1 showed that the use of GAFF for bonded interactions and TraPPE-UA for non-bonded interactions yielded a simulated sample that compared best with available experimental data (wide-angle X-ray scattering and nitrogen adsorption at 77 K). In addition, Brunauer–Emmett–Teller surface areas, calculated from simulated nitrogen isotherms as pseudo-experimental data, showed that this common method failed to measure the geometric surface area of this type of material. These findings are expected to facilitate the predictive screening of different PIM functionalities.     

Water activity coefficients in aqueous amino acid solutions by molecular dynamics simulation: 1. Force field development

New force fields for molecular dynamics (MD) simulation of aqueous zwitterionic amino acid simulations were developed. These were especially designed to calculate activity coefficient of water in amino acid solutions with high accuracy. For example, aqueous solutions of the following amino acids were considered: glycine, alanine, α-aminobutyric acid, α-aminovalerianic acid, valine and leucine. The force fields were obtained by quantum chemical calculations using B3LYP/6-31G and MP2/6-311(d,p) model theories in combination with the Merz–Kollmann–Singh scheme. To further increase the accuracy of the force field, a polarised continuum was considered in all quantum chemical calculations. Water activity coefficients obtained from MD using different all-purpose literature force fields, namely, OPLS, AMBER ff03 and GROMOS 53A6 as well as experimental data are compared with the results utilising the new force field. The new force field is shown to give better results compared with experimental data than existing force fields.     

Comparison of a QM/MM force field and molecular mechanics force fields in simulations of alanine and glycine "dipeptides" (Ace-Ala-Nme and Ace-Gly-Nme) in water in relation to the problem of modeling the unfolded peptide backbone in solution

We compare the conformational distributions of Ace-Ala-Nme and Ace-Gly-Nme sampled in long simulations with several molecular mechanics (MM) force fields and with a fast combined quantum mechanics/molecular mechanics (QM/MM) force field, in which the solute's intramolecular energy and forces are calculated with the self-consistent charge density functional tight binding method (SCCDFTB), and the solvent is represented by either one of the well-known SPC and TIP3P models. All MM force fields give two main states for Ace-Ala-Nme, beta and alpha separated by free energy barriers, but the ratio in which these are sampled varies by a factor of 30, from a high in favor of beta of 6 to a low of 1/5. The frequency of transitions between states is particularly low with the amber and charmm force fields, for which the distributions are noticeably narrower, and the energy barriers between states higher. The lower of the two barriers lies between alpha and beta at values of psi near 0 for all MM simulations except for charmm22. The results of the QM/MM simulations vary less with the choice of MM force field; the ratio beta/alpha varies between 1.5 and 2.2, the easy pass lies at psi near 0, and transitions between states are more frequent than for amber and charmm, but less frequent than for cedar. For Ace-Gly-Nme, all force fields locate a diffuse stable region around phi = pi and psi = pi, whereas the amber force field gives two additional densely sampled states near phi = +/-100 degrees and psi = 0, which are also found with the QM/MM force field. For both solutes, the distribution from the QM/MM simulation shows greater similarity with the distribution in high-resolution protein structures than is the case for any of the MM simulations.     

Properties of Polycrystalline Diamond: Multiscale Modeling Approach

Abstract

Two modeling techniques to characterize fracture behavior of polycrystalline diamond films are discussed. The first technique is a multiscale modeling method in which first-principles local density approximation calculations on selected structures are combined with an analytic mesoscale model to obtain energies and cleavage fracture energies for symmetric ?001? tilt grain boundaries (GBs) over the entire misorientation range. The second technique is large-scale atomistic simulation of the dynamics of failure in notched polycrystalline diamond samples under an applied strain. Electronic characteristics of selected ?001? symmetrical tilt GBs calculated with a semiempirical tight-binding Hamiltonian are also presented, and the possible role of graphitic defects on field emission from polycrystalline diamond is briefly discussed.     

Vibrational analysis of peptides, polypeptides, and proteins. XXXII. alpha-Poly(L-glutamic acid)

The Raman and ir spectra of α-helical poly(L -glutamic acid) have been assigned on the basis of a normal mode calculation for this structure. The force field was based on our previously refined main-chain force constants for α-poly(L -alanine) and side-chain force constants for β-calcium–poly(L -glutamate). Despite the identical backbone α-helical structures, significantly different frequencies are calculated, and observed, in the amide III and backbone stretch regions of α-poly(L -glutamic acid), as compared with α-poly(L -alanine). This clearly demonstrates the influence of side-chain structure on mainchain vibrational modes.     

Effects of different force fields on the structural character of α synuclein β-hairpin peptide (35–56) in aqueous environment

The hallmark of Parkinson’s disease (PD) is the intracellular protein aggregation forming Lewy Bodies (LB) and Lewy neuritis which comprise mostly of a protein, alpha synuclein (α-syn). Molecular dynamics (MD) simulation methods can augment experimental techniques to understand misfolding and aggregation pathways with atomistic resolution. The quality of MD simulations for proteins and peptides depends greatly on the accuracy of empirical force fields. The aim of this work is to investigate the effects of different force fields on the structural character of β hairpin fragment of α-syn (residues 35–56) peptide in aqueous solution. Six independent MD simulations are done in explicit solvent using, AMBER03, AMBER99SB, GROMOS96 43A1, GROMOS96 53A6, OPLS-AA, and CHARMM27 force fields with CMAP corrections. The performance of each force field is assessed from several structural parameters such as root mean square deviation (RMSD), root mean square fluctuation (RMSF), radius of gyration (Rg), solvent accessible surface area (SASA), formation of β-turn, the stability of folded β-hairpin structure, and the favourable conformations obtained for different force fields. In this study, CMAP correction of CHARMM27 force field is found to overestimate the helical conformation, while GROMOS96 53A6 is found to most successfully capture the conformational dynamics of α-syn β-hairpin fragment as elicited from NMR.     

部分力偏倚蒙特卡罗方法和它在计算相对自由能中的应用

提出了一种部分力偏倚蒙特卡罗模拟方法,其中力偏倚仅用于体系中每个分子质心的平动而不用于分子的转动.此方法的总体效益较之先前几种方案为佳,同时报导了该方法在自由能微扰计算中的具体应用.     

Molecular dynamics simulations of structural instability of fullerene family under tension force

Fullerene molecules are cage-like nanoscopic structures with pentagonal and hexagonal faces. In practical applications such as fullerene-reinforced nanocomposites (FRNCs), these structures may be subjected to tension force. In this research, we employ molecular dynamics (MD) simulation to compute the behaviour and deformation of different fullerene molecules, ranging from C₆₀ to C₂₀₀₀, under tension force. To model the interactions between carbon atoms in the MD simulations, the adaptive intermolecular reactive bond order (AIREBO) force field is used. The displacement–force and the displacement–strain energy curves are obtained. It is observed that a new type of structural instability occurs in the fullerene molecules when the applied tension force increases. This abnormal structural instability in the fullerenes is investigated for the first time in the literature. The critical tensile forces and the corresponding mode shapes are determined for different fullerenes. The results indicate that the critical forces and deformations strongly depend upon the number of carbon atoms.     

Force field dependence of NMR-Based,restrained molecular dynamics DNA structure calculations including an analysis of the influence of residual dipolar coupling restraints

Restrained molecular dynamics is widely used to calculate DNA structures from NMR data. Here, results of an in silico experiment show that the force field can be significant compared to the NMR restraints in driving the final structures to converge. Specifically, we observed that i) the influence of the force field leads to artificially tight convergence within final families of structures and ii) the precision and character of resulting structures depend on the choice of force field used in the calculations. A canonical B-DNA model was used as a target structure. Distances, dihedral angles, and simulated residual dipolar couplings were measured in the target structure and used as restraints. X-PLOR and Discover, which use force fields developed for CHARMM and AMBER programs, respectively, were tested and found to produce different final structures despite the use of identical distance and dihedral restraints. Incorporation of residual dipolar coupling restraints in X-PLOR improves convergence with the target structure and between families of structures indicating that the force field dependence can potentially be overcome if residual dipolar coupling restraints are employed.     

Simulating the vapour–liquid equilibria of 1,4-dioxane

1,4-dioxane, a cyclic ether, is an emerging contaminant which is difficult to remove from water with conventional water treatment methods and resistant to biodegradation. Once a reliable force field is developed for 1,4-dioxane, molecular simulation techniques can be useful to study alternative adsorbents for its removal. For this purpose, we carried out Monte Carlo simulations in a constant volume Gibbs Ensemble to generate a force field which is capable of predicting the vapour–liquid coexistence curve and critical data of 1,4-dioxane. Results are given in comparison with experimental data and results from simulations with other force fields. Liquid densities and critical temperature are predicted in excellent agreement with experimental data using the new force field. At high temperatures, predicted vapour densities are in good agreement with experimental data, however, at lower temperatures the predicted vapour densities deviate about an order of magnitude from the experimental values. The critical density is slightly underestimated with our new force field. However, overall, the results of simulations with the new parameters give much better agreement with experimental data compared to the results obtained using other force fields.     

Atom–atom force field for simulation of zirconia bulk,nanosheets and nanotubes

The set of three pair atom–atom potentials for zirconia crystals and nano objects was developed. Besides Buckingham potential and Coulomb interaction, these force fields employ the additional functional forms: Fermi-Dirac, Inverse Gaussian or Morse potential. The developed force fields are capable to reproduce the structural, mechanical and thermodynamic properties of five zirconia phases: monoclinic, tetragonal, cubic, brookite and cotunnite with acceptable accuracy. The proposed force fields give monoclinic phase as a ground state structure and predict the correct energy ordering of the phases. The Raman and IR phonon frequencies, temperature dependencies of the thermodynamic properties of monoclinic and tetragonal zirconia are predicted in good and moderate agreement with the experimentally observed data and results of first-principles calculations. To test the transferability of proposed potentials, the zirconia-based nanosheets and nanotubes were simulated. Results of structure optimisation and calculation of the nanotubes’ strain energies are well compared with the corresponding data of the first-principles calculations.     

The effect of protein conformational flexibility on the electronic properties of a chromophore

In this paper we address the question of how a protein environment can modulate the absorption spectrum of a chromophore during a molecular dynamics simulation. The effect of the protein is modeled as an external field acting on the unperturbed eigenstates of the chromophore. Using a first-principles method recently developed in our group, we calculated the perturbed electronic energies for each frame and the corresponding wavelength absorption during the simulation. We apply this method to a nanosencond timescale molecular dynamics simulation of the light-harvesting peridinin-chlorophyll-protein complex from Amphidinium carterae, where chlorophyll was selected among the chromophores of the complex for the calculation. The combination of this quantum-classical calculation with the analysis of the large amplitude motions of the protein makes it possible to point out the relationship between the conformational flexibility of the environment and the excitation wavelength of the chromophore. Results support the idea of the existence of a correlation between protein conformational flexibility and chlorophyll electronic transitions induced by light.     

Determination method of the balance of the secondary-structure-forming tendencies of force fields

We propose a new method of optimisation of backbone torsion-energy parameters in the force field for molecular simulations of protein systems. This method is based on the idea of balancing the secondary-structure-forming tendencies, namely, those of α-helix and β-sheet structures. We perform a minimisation of the backbone dihedral angle-based root-mean-square deviation of the helix and β structure regions in many protein structures. As an example, we optimised the backbone torsion-energy parameters of AMBER parm96 force field using 100 protein molecules from the Protein Data Bank. We then performed folding simulations of α-helical and β-hairpin peptides, using the optimised force field. The results imply that the new force-field parameters give structures more consistent with the experimental implications than the original AMBER parm96 force field.