Histidine hydrogen bonding in MHC at pH 5 and pH 7 modeled by molecular docking and molecular dynamics simulations

Hydrogen bonds play important roles in maintaining the structure of proteins and in the formation of most biomolecular protein-ligand complexes. All amino acids can act as hydrogen bond donors and acceptors. Among amino acids, Histidine is unique, as it can exist in neutral or positively charged forms within the physiological pH range of 5.0 to 7.0. Histidine can thus interact with other aromatic residues as well as forming hydrogen bonds with polar and charged residues. The ability of His to exchange a proton lies at the heart of many important functional biomolecular interactions, including immunological ones. By using molecular docking and molecular dynamics simulation, we examine the influence of His protonation/deprotonation on peptide binding affinity to MHC class II proteins from locus HLA-DP. Peptide-MHC interaction underlies the adaptive cellular immune response, upon which the next generation of commercially-important vaccines will depend. Consistent with experiment, we find that peptides containing protonated His residues bind better to HLA-DP proteins than those with unprotonated His. Enhanced binding at pH 5.0 is due, in part, to additional hydrogen bonds formed between peptide His+ and DP proteins. In acidic endosomes, protein His79β is predominantly protonated. As a result, the peptide binding cleft narrows in the vicinity of His79β, which stabilizes the peptide - HLA-DP protein complex. © 2014 Bentham Science Publishers.

Distribution of neutral prilocaine in a phospholipid bilayer: Insights from molecular dynamics simulations

In this work, we report a 20-ns constant pressure molecular dynamics simulation of prilocaine (PLC), in amine-amide local anesthetic, in a hydrated liquid crystal bilayer of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine. The partition of PLC induces the lateral expansion of the bilayer and a concomitant contraction in its thickness. PLC molecules are preferentially found in the hydrophobic acyl chains region, with a maximum probability at similar to 12 angstrom from the center of the bilayer (between the C(4) and C(5) methylene groups). A decrease in the acyl chain segmental order parameter, vertical bar S-CD vertical bar, compared to neat bilayers, is found, in good agreement with experimental H-2-NMR studies. The decrease in vertical bar S-CD vertical bar induced by PLC is attributed to a larger accessible volume per lipid in the acyl chain region. (C) 2008 Wiley Periodicals, Inc.

Preferential location of prilocaine and etidocaine in phospholipid bilayers: A molecular dynamics study

In this work, we report a 20-ns constant pressure molecular dynamics simulation of the uncharged form of two amino-amide local anesthetics (LA). etidocaine and prilocaine, present at 1:3 LA:lipid, molar ratio inside the membrane, in the hydrated liquid crystal bilayer phase of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (POPC). Both LAs induced lateral expansion and a concomitant contraction in the bilayer thickness. A decrease in the acyl chain segment order parameter, -S(CD), compared to neat bilayers, was also observed. Besides, both LA molecules got preferentially located in the hydrophobic acyl chains region, with a maximum probability at similar to 12 and similar to 10 angstrom from the center of the bilayer for prilocaine and etidocaine, respectively. (C) 2009 Elsevier B.V. All rights reserved.

Molecular dynamics modeling and simulation of void growth in two dimensions

The mechanisms of growth of a circular void by plastic deformation were studied by means of molecular dynamics in two dimensions (2D). While previous molecular dynamics (MD) simulations in three dimensions (3D) have been limited to small voids (up to ≈10 nm in radius), this strategy allows us to study the behavior of voids of up to 100 nm in radius. MD simulations showed that plastic deformation was triggered by the nucleation of dislocations at the atomic steps of the void surface in the whole range of void sizes studied. The yield stress, defined as stress necessary to nucleate stable dislocations, decreased with temperature, but the void growth rate was not very sensitive to this parameter. Simulations under uniaxial tension, uniaxial deformation and biaxial deformation showed that the void growth rate increased very rapidly with multiaxiality but it did not depend on the initial void radius. These results were compared with previous 3D MD and 2D dislocation dynamics simulations to establish a map of mechanisms and size effects for plastic void growth in crystalline solids.

T-cell epitope prediction and immune complex simulation using molecular dynamics:state of the art and persisting challenges

Atomistic Molecular Dynamics provides powerful and flexible tools for the prediction and analysis of molecular and macromolecular systems. Specifically, it provides a means by which we can measure theoretically that which cannot be measured experimentally: the dynamic time-evolution of complex systems comprising atoms and molecules. It is particularly suitable for the simulation and analysis of the otherwise inaccessible details of MHC-peptide interaction and, on a larger scale, the simulation of the immune synapse. Progress has been relatively tentative yet the emergence of truly high-performance computing and the development of coarse-grained simulation now offers us the hope of accurately predicting thermodynamic parameters and of simulating not merely a handful of proteins but larger, longer simulations comprising thousands of protein molecules and the cellular scale structures they form. We exemplify this within the context of immunoinformatics.

PARAMETRIC STUDY OF REAXFF SIMULATION PARAMETERS FOR MOLECULAR DYNAMICS MODELING OF REACTIVE CARBON GASES

Abstract The development of innovative carbon-based materials can be greatly facilitated by molecular modeling techniques. Although the Reax Force Field (ReaxFF) can be used to simulate the chemical behavior of carbon-based systems, the simulation settings required for accurate predictions have not been fully explored. Using the ReaxFF, molecular dynamics (MD) simulations are used to simulate the chemical behavior of pure carbon and hydrocarbon reactive gases that are involved in the formation of carbon structures such as graphite, buckyballs, amorphous carbon, and carbon nanotubes. It is determined that the maximum simulation time step that can be used in MD simulations with the ReaxFF is dependent on the simulated temperature and selected parameter set, as are the predicted reaction rates. It is also determined that different carbon-based reactive gases react at different rates, and that the predicted equilibrium structures are generally the same for the different ReaxFF parameter sets, except in the case of the predicted formation of large graphitic structures with the Chenoweth parameter set under specific conditions.

Polarization effects in molecular dynamics simulations of glass-formers Ca(NO(3))(2)center dot nH(2)O, n=4, 6, and 8

Thermodynamics, equilibrium structure, and dynamics of glass-forming liquids Ca(NO(3))(2)center dot nH(2)O, n=4, 6, and 8, have been investigated by molecular dynamics (MD) simulations. A polarizable model was considered for H(2)O and NO(3)- on the basis of previous fluctuating charge models for pure water and the molten salt 2Ca(NO(3))(2)center dot 3KNO(3). Similar thermodynamic properties have been obtained with nonpolarizable and polarizable models. The glass transition temperature, T(g), estimated from MD simulations was dependent on polarization, in particular the dependence of T(g) with electrolyte concentration. Significant polarization effects on equilibrium structure were observed in cation-cation, cation-anion, and water-water structures. Polarization increases the diffusion coefficient of H(2)O, but does not change significantly the diffusion coefficients of ions. Viscosity decreases upon inclusion of polarization, but the conductivity calculated with the polarizable model is smaller than the nonpolarizable model because polarization enhances anion-cation interactions.

Thermal diffusion in molecular dynamics simulations of thin film diamond deposition

Molecular dynamics simulations of carbon atom depositions are used to investigate energy diffusion from the impact zone. A modified Stillinger-Weber potential models the carbon interactions for both sp2 and sp3 bonding. Simulations were performed on 50 eV carbon atom depositions onto the (111) surface of a 3.8 x 3.4 x 1.0 nm diamond slab containing 2816 atoms in 11 layers of 256 atoms each. The bottom layer was thermostated to 300 K. At every 100th simulation time step (27 fs), the average local kinetic energy, and hence local temperature, is calculated. To do this the substrate is divided into a set of 15 concentric hemispherical zones, each of thickness one atomic diameter (0.14 nm) and centered on the impact point. A 50-eV incident atom heats the local impact zone above 10 000 K. After the initial large transient (200 fs) the impact zone has cooled below 3000 K, then near 1000 K by 1 ps. Thereafter the temperature profile decays approximately as described by diffusion theory, perturbed by atomic scale fluctuations. A continuum model of classical energy transfer is provided by the traditional thermal diffusion equation. The results show that continuum diffusion theory describes well energy diffusion in low energy atomic deposition processes, at distance and time scales larger than 1.5 nm and 1-2 ps, beyond which the energy decays essentially exponentially. (C) 1998 Published by Elsevier Science S.A. All rights reserved.

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.

NO rebinding to myoglobin: a reactive molecular dynamics study.

The rebinding of NO to myoglobin after photolysis is studied using the 'reactive molecular dynamics' method. In this approach the energy of the system is evaluated on two potential energy surfaces that include the heme-ligand interactions which change between liganded and unliganded myoglobin. This makes it possible to take into account in a simple way, the high dimensionality of the transition seam connecting the reactant and product states. The dynamics of the dissociated NO molecules are examined, and the geometrical and energetic properties of the transition seam are studied. Analysis of the frequency of recrossing shows that the height of the effective rebinding barrier is dependent on the time after photodissociation. This effect is due mainly to protein relaxation and may contribute to the experimentally observed non-exponential rebinding rate of NO, as has been suggested previously.

Exploring unfrequent events with accelerated molecular dynamics simulations: from photochemistry to drug design

La meva incorporació al grup de recerca del Prof. McCammon (University of California San Diego) en qualitat d’investigador post doctoral amb una beca Beatriu de Pinós, va tenir lloc el passat 1 de desembre de 2010; on vaig dur a terme les meves tasques de recerca fins al darrer 1 d’abril de 2012. El Prof. McCammon és un referent mundial en l’aplicació de simulacions de dinàmica molecular (MD) en sistemes biològics d’interès humà. La contribució més important del Prof. McCammon en la simulació de sistemes biològics és el desenvolupament del mètode de dinàmiques moleculars accelerades (AMD). Les simulacions MD convencionals, les quals estan limitades a l’escala de temps del nanosegon (~10-9s), no son adients per l’estudi de sistemes biològics rellevants a escales de temps mes llargues (μs, ms...). AMD permet explorar fenòmens moleculars poc freqüents però que son clau per l’enteniment de molts sistemes biològics; fenòmens que no podrien ser observats d’un altre manera. Durant la meva estada a la “University of California San Diego”, vaig treballar en diferent aplicacions de les simulacions AMD, incloent fotoquímica i disseny de fàrmacs per ordinador. Concretament, primer vaig desenvolupar amb èxit una combinació dels mètodes AMD i simulacions Car-Parrinello per millorar l’exploració de camins de desactivació (interseccions còniques) en reaccions químiques fotoactivades. En segon lloc, vaig aplicar tècniques estadístiques (Replica Exchange) amb AMD en la descripció d’interaccions proteïna-lligand. Finalment, vaig dur a terme un estudi de disseny de fàrmacs per ordinador en la proteïna-G Rho (involucrada en el desenvolupament de càncer humà) combinant anàlisis estructurals i simulacions AMD. Els projectes en els quals he participat han estat publicats (o estan encara en procés de revisió) en diferents revistes científiques, i han estat presentats en diferents congressos internacionals. La memòria inclosa a continuació conté més detalls de cada projecte esmentat.

How T cell receptors interact with peptide-MHCs: A multiple steered molecular dynamics study.

The T-cell receptor (TCR) interaction with antigenic peptides (p) presented by the major histocompatibility complex (MHC) molecule is a key determinant of immune response. In addition, TCR-pMHC interactions offer examples of features more generally pertaining to protein-protein recognition: subtle specificity and cross-reactivity. Despite their importance, molecular details determining the TCR-pMHC binding remain unsolved. However, molecular simulation provides the opportunity to investigate some of these aspects. In this study, we perform extensive equilibrium and steered molecular dynamics simulations to study the unbinding of three TCR-pMHC complexes. As a function of the dissociation reaction coordinate, we are able to obtain converged H-bond counts and energy decompositions at different levels of detail, ranging from the full proteins, to separate residues and water molecules, down to single atoms at the interface. Many observed features do not support a previously proposed two-step model for TCR recognition. Our results also provide keys to interpret experimental point-mutation results. We highlight the role of water both in terms of interface resolvation and of water molecules trapped in the bound complex. Importantly, we illustrate how two TCRs with similar reactivity and structures can have essentially different binding strategies. Proteins 2011; © 2011 Wiley-Liss, Inc.

Brownian dynamics simulation of DNA condensation.

DNA condensation observed in vitro with the addition of polyvalent counterions is due to intermolecular attractive forces. We introduce a quantitative model of these forces in a Brownian dynamics simulation in addition to a standard mean-field Poisson-Boltzmann repulsion. The comparison of a theoretical value of the effective diameter calculated from the second virial coefficient in cylindrical geometry with some experimental results allows a quantitative evaluation of the one-parameter attractive potential. We show afterward that with a sufficient concentration of divalent salt (typically approximately 20 mM MgCl(2)), supercoiled DNA adopts a collapsed form where opposing segments of interwound regions present zones of lateral contact. However, under the same conditions the same plasmid without torsional stress does not collapse. The condensed molecules present coexisting open and collapsed plectonemic regions. Furthermore, simulations show that circular DNA in 50% methanol solutions with 20 mM MgCl(2) aggregates without the requirement of torsional energy. This confirms known experimental results. Finally, a simulated DNA molecule confined in a box of variable size also presents some local collapsed zones in 20 mM MgCl(2) above a critical concentration of the DNA. Conformational entropy reduction obtained either by supercoiling or by confinement seems thus to play a crucial role in all forms of condensation of DNA.

Structural modelling of liquid ZrO2 and glassy B2O3 by using AB INITIO Molecular Dynamics method

Atomic structure of ZrO2 and B2O3 was investigated in this work. New data under extreme conditions (T = 3100 K) was obtained for the liquid ZrO2 structure. A fractional number of boron was investigated for glassy structure of B2O3. It was shown that it is possible to obtain an agreement for the fractional number between NMR and DFT techniques using a suitable initial configuration.

Selectivity of bis(calix[4]diquinone) ionophores towards metal ions in solvent dimethylsulfoxide: a molecular mechanics and molecular dynamics study

Molecular modelling studies have been carried out on two bis(calix[4]diqu(inone) ionophores, each created from two (calix[4]diquinone)arenes bridged at their bottom rims via alkyl chains (CH2)(n), 1: n = 3, 2; n = 4, in order to understand the reported selectivity of these ligands towards different sized metal ions such as Na+, K+, Rb+, and Cs+ in dmso solution. Conformational. analyses have been carried out which show that in the lowest energy conformations of the two macrocycles, the individual calix[4]diquinones exhibit a combination of partial cone, 1,3-alternate and cone conformations. The interactions of these alkali metals with the macrocycles have been studied in the gas phase and in a periodic box of solvent dmso by molecular mechanics and molecular dynamics calculations. Molecular mechanics calculations have been carried out on the mode of entry of the ions into the macrocycles and suggest that this is likely to occur from the side of the central cavity, rather than through the main axis of the calix[4]diquinones. There are energy barriers of ca. 19 kcal mol(-1) for this entry path in the gas phase, but in solution no energy barrier is found. Molecular dynamics simulations show that in both 1 and 2, though particularly in the latter macrocycle, one or two solvent molecules are bonded to the metal throughout the course of the simulation, often to the exclusion, of one or more of the ether oxygen atoms. By contrast the carbonyl oxygen atoms remain bonded to the metal atoms throughout with bond lengths that remain significantly less than those to the ether oxygen atoms. Free energy perturbation studies have been carried out in dmso and indicate that for 1, the selectivity follows the order Rb+ approximate to K+ > Cs+ >> Na+, which is partially in agreement with the experimental results. The energy differences are small and indeed the ratio between stability constants found for Cs+ and K+ complexes is only 0.60, showing that 1 has only a slight preference for K+. For the larger receptor 2, which is better suited to metal complexation, the binding affinity follows the pattern Cs+ >> Rb+ >> K+ >> Na+, with energy differences of 5.75, 2.61, 2.78 kcal mol(-1) which is perfectly consistent with experimental results.