The role of disulfide bridges in the 3-D structures of the antimicrobial peptides gomesin and protegrin-1: a molecular dynamics study

Some antimicrobial peptides have a broad spectrum of action against many different kinds of microorganisms. Gomesin and protegrin-1 are examples of such antimicrobial peptides, and they were studied by molecular dynamics in this research. Both have a beta-hairpin conformation stabilized by two disulfide bridges and are active against Gram-positive and Gram-negative bacteria, as well as fungi. In this study, the role of the disulfide bridge in the maintenance of the tertiary peptide structure of protegrin-1 and gomesin is analyzed by the structural characteristics of these peptides and two of their respective variants, gomy4 and proty4, in which the four cysteines are replaced by four tyrosine residues. The absence of disulfide bridges in gomy4 and proty4 is compensated by overall reinforcement of the original hydrogen bonds and extra attractive interactions between the aromatic rings of the tyrosine residues. The net effects on the variants with respect to the corresponding natural peptides are: i) maintenance of the original beta-hairpin conformation, with great structural similarities between the mutant and the corresponding natural peptide; ii) combination of positive F and. Ramachandran angles within the hairpin head region with a qualitative change to a combination of positive (F) and negative (.) angles, and iii) significant increase in structural flexibility. Experimental facts about the antimicrobial activity of the gomesin and protegrin-1 variants have also been established here, in the hope that the detailed data provided in the present study may be useful for understanding the mechanism of action of these peptides.

Dynamic polarizability, Cauchy moments, and the optical absorption spectrum of liquid water: A sequential molecular dynamics/quantum mechanical approach

The dynamic polarizability and optical absorption spectrum of liquid water in the 6-15 eV energy range are investigated by a sequential molecular dynamics (MD)/quantum mechanical approach. The MD simulations are based on a polarizable model for liquid water. Calculation of electronic properties relies on time-dependent density functional and equation-of-motion coupled-cluster theories. Results for the dynamic polarizability, Cauchy moments, S(-2), S(-4), S(-6), and dielectric properties of liquid water are reported. The theoretical predictions for the optical absorption spectrum of liquid water are in good agreement with experimental information.

Molecular dynamics simulation of liquid trimethylphosphine

Structural and dynamical properties of liquid trimethylphosphine (TMP), (CH(3))(3)P, as a function of temperature is investigated by molecular dynamics (MD) simulations. The force field used in the MD simulations, which has been proposed from molecular mechanics and quantum chemistry calculations, is able to reproduce the experimental density of liquid TMP at room temperature. Equilibrium structure is investigated by the usual radial distribution function, g(r), and also in the reciprocal space by the static structure factor, S(k). On the basis of center of mass distances, liquid TMP behaves like a simple liquid of almost spherical particles, but orientational correlation due to dipole-dipole interactions is revealed at short-range distances. Single particle and collective dynamics are investigated by several time correlation functions. At high temperatures, diffusion and reorientation occur at the same time range as relaxation of the liquid structure. Decoupling of these dynamic properties starts below ca. 220 K, when rattling dynamics of a given TMP molecules due to the cage effect of neighbouring molecules becomes important. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3624408]

Structure optimization combining soft-core interaction functions, the diffusion equation method, and molecular dynamics

Smoothing the potential energy surface for structure optimization is a general and commonly applied strategy. We propose a combination of soft-core potential energy functions and a variation of the diffusion equation method to smooth potential energy surfaces, which is applicable to complex systems such as protein structures; The performance of the method was demonstrated by comparison with simulated annealing using the refinement of the undecapeptide Cyclosporin A as a test case. Simulations were repeated many times using different initial conditions and structures since the methods are heuristic and results are only meaningful in a statistical sense.

Molecular dynamics computer simulation of scratch resistance testing of polymers: visualization

Experimental scratch resistance testing provides two numbers: the penetration depth Rp and the healing depth Rh. In molecular dynamics computer simulations, we create a material consisting of N statistical chain segments by polymerization; a reinforcing phase can be included. Then we simulate the movement of an indenter and response of the segments during X time steps. Each segment at each time step has three Cartesian coordinates of position and three of momentum. We describe methods of visualization of results based on a record of 6NX coordinates. We obtain a continuous dependence on time t of positions of each of the segments on the path of the indenter. Scratch resistance at a given location can be connected to spatial structures of individual polymeric chains.

Polymer indentation with mesoscopic molecular dynamics

Indentation tests are used to determine the hardness of a material, e.g., Rockwell, Vickers, or Knoop. The indentation process is empirically observed in the laboratory during these tests; the mechanics of indentation is insufficiently understood. We have performed first molecular dynamics computer simulations of indentation resistance of polymers with a chain structure similar to that of high density polyethylene (HDPE). A coarse grain model of HDPE is used to simulate how the interconnected segments respond to an external force imposed by an indenter. Results include the time-dependent measurement of penetration depth, recovery depth, and recovery percentage, with respect to indenter force, indenter size, and indentation time parameters. The simulations provide results that are inaccessible experimentally, including continuous evolution of the pertinent tribological parameters during the entire indentation process.

Simulating scratch behavior of polymers with mesoscopic molecular dynamics

Part replacement and repair is needed in structures with moving parts because of scratchability and wear. In spite of some accumulation of experimental evidence, scratch resistance is still not well understood. We have applied molecular dynamics to study scratch resistance of amorphous polymeric materials through computer simulations. As a first approach, a coarse grain model was created for high density polyethylene at the mesoscale. We have also extended the traditional approach and used real units rather than reduced units (to our knowledge, for the first time), which enable an improved quantification of simulation results. The obtained results include analysis of penetration depth, residual depth and recovery percentage related to indenter force and size. Our results show there is a clear effect from these parameters on the tribological properties. We also discuss a "crooked smile" effect on the scratched surface and the reasons for its appearance.

Ligand-Induced structural changes in TEM‑1 probed by molecular dynamics and relative binding free energy calculations

The TEM family of enzymes has had a crucial impact on the pharmaceutical industry due to their important role in antibiotic resistance. Even with the latest technologies in structural biology and genomics, no 3D structure of a TEM- 1/antibiotic complex is known previous to acylation. Therefore, the comprehension of their capability in acylate antibiotics is based on the protein macromolecular structure uncomplexed. In this work, molecular docking, molecular dynamic simulations, and relative free energy calculations were applied in order to get a comprehensive and thorough analysis of TEM-1/ampicillin and TEM-1/amoxicillin complexes. We described the complexes and analyzed the effect of ligand binding on the overall structure. We clearly demonstrate that the key residues involved in the stability of the ligand (hot-spots) vary with the nature of the ligand. Structural effects such as (i) the distances between interfacial residues (Ser70−Oγ and Lys73−Nζ, Lys73−Nζ and Ser130−Oγ, and Ser70−Oγ−Ser130−Oγ), (ii) side chain rotamer variation (Tyr105 and Glu240), and (iii) the presence of conserved waters can be also influenced by ligand binding. This study supports the hypothesis that TEM-1 suffers structural modifications upon ligand binding.

The effects of solvent composition on the affinity of a peptide towards hair keratin: experimental and molecular dynamics data

The study of the interaction between hair filaments and formulations or peptides is of utmost importance in fields like cosmetic research. Keratin intermediate filaments structure is not fully described, limiting the molecular dynamics (MD) studies in this field although its high potential to improve the area. We developed a computational model of a truncated protofibril, simulated its behavior in alcoholic based formulations and with one peptide. The simulations showed a strong interaction between the benzyl alcohol molecules of the formulations and the model, leading to the disorganization of the keratin chains, which regress with the removal of the alcohol molecules. This behavior can explain the increase of peptide uptake in hair shafts evidenced in fluorescence microscopy pictures. The model developed is valid to computationally reproduce the interaction between hair and alcoholic formulations and provide a robust base for new MD studies about hair properties. It is shown that the MD simulations can improve hair cosmetic research, improving the uptake of a compound of interest.

Protein-protein interaction investigated by steered molecular dynamics: the TCR-pMHC complex.

We present a novel steered molecular dynamics scheme to induce the dissociation of large protein-protein complexes. We apply this scheme to study the interaction of a T cell receptor (TCR) with a major histocompatibility complex (MHC) presenting a peptide (p). Two TCR-pMHC complexes are considered, which only differ by the mutation of a single amino acid on the peptide; one is a strong agonist that produces T cell activation in vivo, while the other is an antagonist. We investigate the interaction mechanism from a large number of unbinding trajectories by analyzing van der Waals and electrostatic interactions and by computing energy changes in proteins and solvent. In addition, dissociation potentials of mean force are calculated with the Jarzynski identity, using an averaging method developed for our steering scheme. We analyze the convergence of the Jarzynski exponential average, which is hampered by the large amount of dissipative work involved and the complexity of the system. The resulting dissociation free energies largely underestimate experimental values, but the simulations are able to clearly differentiate between wild-type and mutated TCR-pMHC and give insights into the dissociation mechanism.

Influence of ionization state on the activation of temocapril by hCES1: a molecular-dynamics study.

Temocapril is a prodrug whose hydrolysis by carboxylesterase 1 (CES1) yields the active ACE inhibitor temocaprilat. This molecular-dynamics (MD) study uses a resolved structure of the human CES1 (hCES1) to investigate some mechanistic details of temocapril hydrolysis. The ionization constants of temocapril (pK1 and pK3) and temocaprilat (pK1, pK2, and pK3) were determined experimentally and computationally using commercial algorithms. The constants so obtained were in good agreement and revealed that temocapril exists mainly in three ionic forms (a cation, a zwitterion, and an anion), whereas temocaprilat exists in four major ionic forms (a cation, a zwitterion, an anion, and a dianion). All these ionic forms were used as ligands in 5-ns MS simulations. While the cationic and zwitterionic forms of temocapril were involved in an ion-pair bond with Glu255 suggestive of an inhibitor behavior, the anionic form remained in a productive interaction with the catalytic center. As for temocaprilat, its cation appeared trapped by Glu255, while its zwitterion and anion made a slow departure from the catalytic site and a partial egress from the protein. Only its dianion was effectively removed from the catalytic site and attracted to the protein surface by Lys residues. A detailed mechanism of product egress emerges from the simulations.

Linear viscoelasticity from molecular dynamics simulation of entangled polymers

The linear viscoelastic (LVE) spectrum is one of the primary fingerprints of polymer solutions and melts, carrying information about most relaxation processes in the system. Many single chain theories and models start with predicting the LVE spectrum to validate their assumptions. However, until now, no reliable linear stress relaxation data were available from simulations of multichain systems. In this work, we propose a new efficient way to calculate a wide variety of correlation functions and mean-square displacements during simulations without significant additional CPU cost. Using this method, we calculate stress−stress autocorrelation functions for a simple bead−spring model of polymer melt for a wide range of chain lengths, densities, temperatures, and chain stiffnesses. The obtained stress−stress autocorrelation functions were compared with the single chain slip−spring model in order to obtain entanglement related parameters, such as the plateau modulus or the molecular weight between entanglements. Then, the dependence of the plateau modulus on the packing length is discussed. We have also identified three different contributions to the stress relaxation:  bond length relaxation, colloidal and polymeric. Their dependence on the density and the temperature is demonstrated for short unentangled systems without inertia.

Grid warehousing of molecular dynamics protein unfolding data

With the increasing awareness of protein folding disorders, the explosion of genomic information, and the need for efficient ways to predict protein structure, protein folding and unfolding has become a central issue in molecular sciences research. Molecular dynamics computer simulations are increasingly employed to understand the folding and unfolding of proteins. Running protein unfolding simulations is computationally expensive and finding ways to enhance performance is a grid issue on its own. However, more and more groups run such simulations and generate a myriad of data, which raises new challenges in managing and analyzing these data. Because the vast range of proteins researchers want to study and simulate, the computational effort needed to generate data, the large data volumes involved, and the different types of analyses scientists need to perform, it is desirable to provide a public repository allowing researchers to pool and share protein unfolding data. This paper describes efforts to provide a grid-enabled data warehouse for protein unfolding data. We outline the challenge and present first results in the design and implementation of the data warehouse.

Scanning tunneling microscopy and molecular dynamics study of the Li2TiO3(001) surface

We have investigated the (001) surface structure of lithium titanate (Li2TiO3) using auger electron spectroscopy (AES), low-energy electron diffraction (LEED), and scanning tunneling microscopy (STM). Li2TiO3 is a potential fusion reactor blanket material. After annealing at 1200 K, LEED demonstrated that the Li2TiO3(001) surface was well ordered and not reconstructed. STM imaging showed that terraces are separated in height by about 0.3 nm suggesting a single termination layer. Moreover, hexagonal patterns with a periodicity of ∼0.4 nm are observed. On the basis of molecular dynamics (MD) simulations, these are interpreted as a dynamic arrangement of Li atoms.

Combined Crystallographic and Solution Molecular Dynamics Study of Allosteric Effects in Ester and Ketone p-tert-Butylcalix[4]arene Derivatives and Their Complexes with Acetonitrile, Cd(II), and Pb(II)

We describe here a procedure to bridge the gap in the field of calixarene physicochemistry between solid-state atomic-resolution structural information and the liquid-state low-resolution thermodynamics and spectroscopic data. We use MD simulations to study the kinetics and energetics involved in the complexation of lower rim calix[4]arene derivatives (L), containing bidentate ester (1) and ketone (2) pendant groups, with acetonitrile molecule (MeCN) and Cd2+ and Pb2+ ions (M2+) in acetonitrile solution. On one hand, we found that the prior inclusion of MeCN into the calix to form a L(MeCN) adduct has only a weak effect in preorganizing the hydrophilic cavity toward metal ion binding. On the other hand, the strong ion-hydrophilic cavity interaction produces a wide open calix which enhances the binding of one MeCN molecule (allosteric effect) to stabilize the whole (M2+)1(MeCN) bifunctional complex. We reach two major conclusions: (i) the MD results for the (M2+)1(MeCN) binding are in close agreement with the ""endo"", fully encapsulated, metal complex found by X-ray diffraction and in vacuo MD calculations, and (ii) the MD structure for the more flexible 2 ligand, however, differs from the also endo solid-state molecule. In fact, it shows strong solvation effects at the calixarene lower bore by competing MeCN molecules that share the metal coordination sphere with the four C=O oxygens of an ""exo"" (M2+)2(MeCN) complex.