Halide, Ammonium, and Alkali Metal Ion Parameters for Modeling Aqueous Solutions

OPLS3: A Force Field Providing Broad Coverage of Drug-like Small Molecules and Proteins

The parametrization and validation of the OPLS3 force field for small molecules and proteins are reported. Enhancements with respect to the previous version (OPLS2.1) include the addition of off-atom charge sites to represent halogen bonding and aryl nitrogen lone pairs as well as a complete refit of peptide dihedral parameters to better model the native structure of proteins. To adequately cover medicinal chemical space, OPLS3 employs over an order of magnitude more reference data and associated parameter types relative to other commonly used small molecule force fields (e.g., MMFF and OPLS_2005). As a consequence, OPLS3 achieves a high level of accuracy across performance benchmarks that assess small molecule conformational propensities and solvation. The newly fitted peptide dihedrals lead to significant improvements in the representation of secondary structure elements in simulated peptides and native structure stability over a number of proteins. Together, the improvements made to both the small molecule and protein force field lead to a high level of accuracy in predicting protein-ligand binding measured over a wide range of targets and ligands (less than 1 kcal/mol RMS error) representing a 30% improvement over earlier variants of the OPLS force field.

Effect of Force Field Parameters on Sodium and Potassium Ion Binding to Dipalmitoyl Phosphatidylcholine Bilayers

The behavior of electrolytes in molecular dynamics simulations of zwitterionic phospholipid bilayers is very sensitive to the force field parameters used. Here, several 200 ns molecular dynamics of simulations of dipalmitoyl phosphotidylcholine (PC) bilayers in 0.2 M sodium or potassium chloride using various common force field parameters for the cations are presented. All employed parameter sets give a larger number of Na(+) ions than K(+) ions that bind to the lipid heads, but depending on the parameter choice quite different results are seen. A wide range of coordination numbers for the Na(+) and K(+) ions is also observed. These findings have been analyzed and compared to published experimental data. Some simulations produce aggregates of potassium chloride, indicating (in accordance with published simulations) that these force fields do not reproduce the delicate balance between salt and solvated ions. The differences between the force fields can be characterized by one single parameter, the electrostatic radius of the ion, which is correlated to σMO (M represents Na(+)/K(+)), the Lennard-Jones radius. When this parameter exceeds a certain threshold, binding to the lipid heads is no longer observed. One would, however, need more accurate experimental data to judge or rank the different force fields precisely. Still, reasons for the poor performance of some of the parameter sets are clearly demonstrated, and a quality control procedure is provided.

The Dependence of Amyloid-β Dynamics on Protein Force Fields and Water Models

We studied the dynamics of Aβ40 , involved in Alzheimer's disease, by using 21 methods combined from Amber03, Amber99sb-ILDN, Charmm27, Charmm22*, OPLS-2001, OPLS-2006, OPLS-2008, Gromos96-43a1, Gromos96-53a6, Gromos96-54a7, and the water models SPC, TIP3P, TIP4P. Major differences in the structural ensembles were systematized: Amber03, Charmm27, and Gromos96-54a7 stabilize the helices; Gromos96-43a1 and Gromos53a6 favor the β-strands (with Charmm22* and Amber99sb-ILDN in between), and OPLS produces unstructured ensembles. The accuracy of the NMR chemical shifts was in the order: Charmm22*>Amber99sb-ILDN>OPLS-2008≈Gromos96-43a1>Gromos96-54a7≈OPLS-2001>OPLS-2006>Gromos96-53a6>Charmm27>Amber03. The computed (3) JHNHα -coupling constants were sensitive to experiment type and Karplus parameterization. Overall, the ensembles of Charmm22* and Amber99sb-ILDN provided the best agreement with experimental NMR and circular dichroism data, providing a model for the real Aβ monomer ensemble. Also, the polar water model TIP3P significantly favored helix and compact conformations.

Adsorption Properties of Triethylene Glycol on a Hydrated {101̅4} Calcite Surface and Its Effect on Adsorbed Water

Molecular dynamics (MD) and Born-Oppenheimer MD (BOMD) simulations were employed to investigate adsorption of aqueous triethylene glycol (TEG) on a hydrated {101̅4} calcite surface at 298 K. We analyzed the orientation of TEG adsorbed on calcite, as well as the impact of TEG on the water density and adsorption free energy. The adsorption energies of TEG, free energy profiles for TEG, details of hydrogen bonding between water and adsorbed TEG, and dihedral angle distribution of adsorbed TEG were estimated. We found that while the first layer of water was mostly unaffected by the presence of adsorbed TEG, the density of the second water layer was decreased by 71% at 75% surface coverage of TEG. TEG primarily attached to the calcite surface via two adjacent adsorption sites. Hydrogen bonds between water and adsorbed TEG in the second layer almost exclusively involved the hydroxyl oxygen of TEG. The adsorption energy of TEG on calcite in a vacuum environment calculated by classical MD amounted to 217 kJ/mol, which agreed very well with estimates found by using BOMD. Adsorption on hydrated calcite yielded a drastically lower value of 33 kJ/mol, with the corresponding adsorption free energy of 55.3 kJ/mol, giving an entropy increase of 22.3 kJ/mol due to adsorption. We found that the presence of TEG resulted in a decreased magnitude of the adsorption free energy of water, thus decreasing the calcite wettability. This effect can have a profound effect on oil and gas reservoir properties and must be carefully considered when evaluating the risk of hydrate nucleation.

Small molecule solvation changes due to the presence of salt are governed by the cost of solvent cavity formation and dispersion

We are interested in the free energies of transferring nonpolar solutes into aqueous NaCl solutions with salt concentrations upwards of 2 M, the Hofmeister regime. We use the semi-explicit assembly (SEA) computational model to represent these electrolyte solutions. We find good agreement with experiments (Setschenow coefficients) on 43 nonpolar and polar solutes and with TIP3P explicit-solvent simulations. Besides being much faster than explicit solvent calculations, SEA is more accurate than the PB models we tested, successfully capturing even subtle salt effects in both the polar and nonpolar components of solvation. We find that the salt effects are mainly due to changes in the cost of forming nonpolar cavities in aqueous NaCl solutions, and not mainly due to solute-ion electrostatic interactions.

Crystal structure of the sodium-proton antiporter NhaA dimer and new mechanistic insights

A dimeric structure of the sodium–proton antiporter NhaA provides insight into the roles of Asp163 and Lys300 in the transport mechanism.

Sodium–proton antiporters rapidly exchange protons and sodium ions across the membrane to regulate intracellular pH, cell volume, and sodium concentration. How ion binding and release is coupled to the conformational changes associated with transport is not clear. Here, we report a crystal form of the prototypical sodium–proton antiporter NhaA from Escherichia coli in which the protein is seen as a dimer. In this new structure, we observe a salt bridge between an essential aspartic acid (Asp163) and a conserved lysine (Lys300). An equivalent salt bridge is present in the homologous transporter NapA, but not in the only other known crystal structure of NhaA, which provides the foundation of most existing structural models of electrogenic sodium–proton antiport. Molecular dynamics simulations show that the stability of the salt bridge is weakened by sodium ions binding to Asp164 and the neighboring Asp163. This suggests that the transport mechanism involves Asp163 switching between forming a salt bridge with Lys300 and interacting with the sodium ion. pK_a calculations suggest that Asp163 is highly unlikely to be protonated when involved in the salt bridge. As it has been previously suggested that Asp163 is one of the two residues through which proton transport occurs, these results have clear implications to the current mechanistic models of sodium–proton antiport in NhaA.

Modelling the effects of salt solutions on the hydration of calcium ions

Classical molecular dynamics simulations of several aqueous alkali halide salt solutions have been used to determine the effect of electrolytes on the structure of water and the hydration properties of calcium ions. Compared with the simulations of Ca(2+) ions in pure liquid water, the frequency of water exchange in the first hydration shell of calcium, which is a fundamental process in controlling the reactivity of calcium(ii) aqua-ions, is drastically reduced in the presence of other electrolytes in solution. The strong stabilization of the hydration shell of Ca(2+) occurs not only when the halide anions are directly coordinated to calcium, but also when the alkali and halide ions are placed at or outside the second coordination shell of Ca(2+), suggesting that the reactivity of the first solvation shell of the calcium ion can be influenced by the specific affinity of other ions in solution for the water molecules coordinated to Ca(2+). Analysis of the hydrogen-bonded structure of water in the vicinity of the calcium ion shows that the average number of hydrogen bonds per water molecules, which is 1.8 in pure liquid water, decreases as the concentration of alkali-halide salts in solution increases, and that the temporal fluctuations of hydrogen bonds are significantly larger than those obtained for Ca(2+) in pure liquid water. This effect has been explained in terms of the dynamics of reorganization of the O-H···X(-) (X = F, Cl and Br) hydrogen bond. This work shows the importance of solution composition in determining the hydrogen-bonding network and ligand-exchange dynamics around metal ions, both in solution and at the mineral-water interfaces, which in turn has implications for interactions occurring at the mineral-water interface, ultimately controlling the mobilization of ions in the environment as well as in industrial processes.

Systematic Parameterization of Monovalent Ions Employing the Nonbonded Model

Monovalent ions play fundamental roles in many biological processes in organisms. Modeling these ions in molecular simulations continues to be a challenging problem. The 12-6 Lennard-Jones (LJ) nonbonded model is widely used to model monovalent ions in classical molecular dynamics simulations. A lot of parameterization efforts have been reported for these ions with a number of experimental end points. However, some reported parameter sets do not have a good balance between the two Lennard-Jones parameters (the van der Waals (VDW) radius and potential well depth), which affects their transferability. In the present work, via the use of a noble gas curve we fitted in former work (J. Chem. Theory Comput. 2013, 9, 2733), we reoptimized the 12-6 LJ parameters for 15 monovalent ions (11 positive and 4 negative ions) for three extensively used water models (TIP3P, SPC/E, and TIP4P(EW)). Since the 12-6 LJ nonbonded model performs poorly in some instances for these ions, we have also parameterized the 12-6-4 LJ-type nonbonded model (J. Chem. Theory Comput. 2014, 10, 289) using the same three water models. The three derived parameter sets focused on reproducing the hydration free energies (the HFE set) and the ion-oxygen distance (the IOD set) using the 12-6 LJ nonbonded model and the 12-6-4 LJ-type nonbonded model (the 12-6-4 set) overall give improved results. In particular, the final parameter sets showed better agreement with quantum mechanically calculated VDW radii and improved transferability to ion-pair solutions when compared to previous parameter sets.

Toward the correction of effective electrostatic forces in explicit-solvent molecular dynamics simulations: restraints on solvent-generated electrostatic potential and solvent polarization

Despite considerable advances in computing power, atomistic simulations under nonperiodic boundary conditions, with Coulombic electrostatic interactions and in systems large enough to reduce finite-size associated errors in thermodynamic quantities to within the thermal energy, are still not affordable. As a result, periodic boundary conditions, systems of microscopic size and effective electrostatic interaction functions are frequently resorted to. Ensuing artifacts in thermodynamic quantities are nowadays routinely corrected a posteriori, but the underlying configurational sampling still descends from spurious forces. The present study addresses this problem through the introduction of on-the-fly corrections to the physical forces during an atomistic molecular dynamics simulation. Two different approaches are suggested, where the force corrections are derived from special potential energy terms. In the first approach, the solvent-generated electrostatic potential sampled at a given atom site is restrained to a target value involving corrections for electrostatic artifacts. In the second approach, the long-range regime of the solvent polarization around a given atom site is restrained to the Born polarization, i.e., the solvent polarization corresponding to the ideal situation of a macroscopic system under nonperiodic boundary conditions and governed by Coulombic electrostatic interactions. The restraints are applied to the explicit-water simulation of a hydrated sodium ion, and the effect of the restraints on the structural and energetic properties of the solvent is illustrated. Furthermore, by means of the calculation of the charging free energy of a hydrated sodium ion, it is shown how the electrostatic potential restraint translates into the on-the-fly consideration of the corresponding free-energy correction terms. It is discussed how the restraints can be generalized to situations involving several solute particles. Although the present study considers a very simple system only, it is an important step toward the on-the-fly elimination of finite-size and approximate-electrostatic artifacts during atomistic molecular dynamics simulations.

Investigations of clustering of ions and diffusivity in concentrated aqueous solutions of lithium chloride by molecular dynamic simulations

The diffusion coefficient of Li⁺ ions decreases with increase in LiCl concentration which depends on the size of coordination structure of ions formed in solutions.

A short designed semi-aromatic organic nanotube--synthesis, chiroptical characterization, and host properties

The first generation of an organic nanotube based on the enantiomerically pure bicyclo[3.3.1]nonane framework is presented. The helical tube synthesised is the longest to date having its aromatic systems oriented parallel to the axis of propagation (length ∼26 Å and inner diameter ∼11 Å according to molecular dynamics simulations in chloroform). The synthesis of the tube, a heptamer, is based on a series of Friedländer condensations and the use of pyrido[3,2-d]pyrimidine units as masked 2-amino aldehydes, as a general means to propagate organic tubular structures and the introduction of a methoxy group for modification toward solubility and functionalization are described. The electronic CD spectra of the tube and molecular intermediates are correlated with theoretical spectra calculated with time-dependent density functional theory to characterize the chirality of the tube. Both experimental (NMR-titrations) and theoretical (molecular dynamics simulations) techniques are used to investigate the use of the tube as a receptor for the acetylcholine and guanidinium cations, respectively.

Molecular dynamics simulations of shock waves in hydroxyl-terminated polybutadiene melts: mechanical and structural responses

The mechanical and structural responses of hydroxyl-terminated cis-1,4-polybutadiene melts to shock waves were investigated by means of all-atom non-reactive molecular dynamics simulations. The simulations were performed using the OPLS-AA force field but with the standard 12-6 Lennard-Jones potential replaced by the Buckingham exponential-6 potential to better represent the interactions at high compression. Monodisperse systems containing 64, 128, and 256 backbone carbon atoms were studied. Supported shock waves were generated by impacting the samples onto stationary pistons at impact velocities of 1.0, 1.5, 2.0, and 2.5 km s(-1), yielding shock pressures between approximately 2.8 GPa and 12.5 GPa. Single-molecule structural properties (squared radii of gyration, asphericity parameters, and orientational order parameters) and mechanical properties (density, shock pressure, shock temperature, and shear stress) were analyzed using a geometric binning scheme to obtain spatio-temporal resolution in the reference frame centered on the shock front. Our results indicate that while shear stress behind the shock front is relieved on a ∼0.5 ps time scale, a shock-induced transition to a glass-like state occurs with a concomitant increase of structural relaxation times by several orders of magnitude.

Molecular dynamics simulations of polyamide membrane, calcium alginate gel, and their interactions in aqueous solution

We perform molecular dynamics (MD) simulations to investigate the cross-linked polyamide (PA) membrane, the aggregation of alginate molecules in the presence of Ca(2+) ions, and their molecular binding mechanism in aqueous solution. We use a steered molecular dynamics (SMD) approach to simulate the unbinding process between a PA membrane and an alginate gel complex. Simulation results show that Ca(2+) ions are strongly associated with the carboxylate groups in alginate molecules, forming a web structure. The adhesion force between alginate gel and PA surface during unbinding originates from several important molecular interactions. These include the short-range hydrogen bonding and van der Waals attraction forces, and the ionic bridge binding that extends much longer pulling distances due to the significant chain deformations of alginate gel and PA membrane.

Force field independent metal parameters using a nonbonded dummy model

The cationic dummy atom approach provides a powerful nonbonded description for a range of alkaline-earth and transition-metal centers, capturing both structural and electrostatic effects. In this work we refine existing literature parameters for octahedrally coordinated Mn²⁺, Zn²⁺, Mg²⁺, and Ca²⁺, as well as providing new parameters for Ni²⁺, Co²⁺, and Fe²⁺. In all the cases, we are able to reproduce both M²⁺–O distances and experimental solvation free energies, which has not been achieved to date for transition metals using any other model. The parameters have also been tested using two different water models and show consistent performance. Therefore, our parameters are easily transferable to any force field that describes nonbonded interactions using Coulomb and Lennard-Jones potentials. Finally, we demonstrate the stability of our parameters in both the human and Escherichia coli variants of the enzyme glyoxalase I as showcase systems, as both enzymes are active with a range of transition metals. The parameters presented in this work provide a valuable resource for the molecular simulation community, as they extend the range of metal ions that can be studied using classical approaches, while also providing a starting point for subsequent parametrization of new metal centers.

Combining the Fluctuating Charge Method, Non-Periodic Boundary Conditions and Meta-Dynamics: Aqua Ions as case studies

We present the current status of development of our code for performing Molecular Dynamics (MD) simulations exploiting a polarizable force field based on the Fluctuating Charge (FQ) method and non-Periodic Boundary Conditions (NPBC). Continuing on the path set in a previous work, we increased the capabilities of the code by implementing a number of new features, including: a non-iterative algorithm for rigid trigonal molecule simulations; two additional temperature coupling schemes; a meta-dynamics based approach for effective free energy evaluations. Although these are well known algorithms, each present in one or more widely used MD packages, they have now been tested, for the first time, in the context of the FQ model coupled with NPBC. As case studies, we considered three aqueous ions of increasing charge, namely Na⁺, Ca²⁺ and La³⁺, at infinite dilution. In particular, by exploiting a computational approach recently proposed by our group and based on the metadynamics technique, we focused on the important role played by solvent polarization on ionic hydration structures, also investigating the free energy landscapes of ion coordination and the water exchange rates. Such an approach, previously tested with standard non-polarizable models, was applied here to evaluate the effects of explicit polarization on water exchange barriers between different solvent coordination structures. Moreover, we have analyzed and discussed in some detail non-linear electrostatic effects arising from solvent polarization while going from a mono- to a di- and trivalent ion.

Anomalous water diffusion in salt solutions

The dynamics of water exhibits anomalous behavior in the presence of different electrolytes. Recent experiments [Kim JS, Wu Z, Morrow AR, Yethiraj A, Yethiraj A (2012) J Phys Chem B 116(39):12007-12013] have found that the self-diffusion of water (Dw) can either be enhanced or suppressed around CsI and NaCl, respectively, relative to that of neat water. Here we show that unlike classical empirical potentials, ab initio molecular dynamics simulations successfully reproduce the qualitative trends observed experimentally. These types of phenomena have often been rationalized in terms of the "structure-making" or "structure-breaking" effects of different ions on the solvent, although the microscopic origins of these features have remained elusive. Rather than disrupting the network in a significant manner, the electrolytes studied here cause rather subtle changes in both structural and dynamical properties of water. In particular, we show that water in the ab initio molecular dynamics simulations is characterized by dynamic heterogeneity, which turns out to be critical in reproducing the experimental trends.

Amide I IR probing of core and shell hydrogen-bond structures in reverse micelles

The properties of N-methylacetamide (NMA) molecules encapsulated in the reverse micelles (RMs) formed by anionic surfactant aerosol OT (AOT), are studied with vibrational spectroscopy and computation. Vibrational spectra of the amide I′ mode of the fully deuterated NMA-d7 show gradual increase of peak frequencies and line broadening as the size of RMs decreases. Analyses of the spectral features reveal the presence of three states of NMA-d7 that correspond to NMA located in the core of water phase (absorption frequency of 1606 cm–1) and two types of interfacial NMA near the surfactant layer (1620 and 1644 cm–1). In larger RMs with water content w0 = [D2O]/[AOT] ≥ 10, only the first two states are observed, whereas in smaller RMs, the population of the third state grows up to 25 % at w0 = 2. These results indicate the general validity of the two-state core/shell model for the confined aqueous solution of NMA, with small modifications due to the system-dependent solute-interface interaction. However, simulations of small RM systems with w0 ≤ 15 show continuous variations of the population, frequency shifts, and the solute-solvent interaction strengths at solute-interface distance less than 4 Å. Thus, the distinction of solute core/shell states tends to be blurred in small RMs but is still effective in interpreting the average spectroscopic observables.

Net charge changes in the calculation of relative ligand-binding free energies via classical atomistic molecular dynamics simulation

The calculation of binding free energies of charged species to a target molecule is a frequently encountered problem in molecular dynamics studies of (bio-)chemical thermodynamics. Many important endogenous receptor-binding molecules, enzyme substrates, or drug molecules have a nonzero net charge. Absolute binding free energies, as well as binding free energies relative to another molecule with a different net charge will be affected by artifacts due to the used effective electrostatic interaction function and associated parameters (e.g., size of the computational box). In the present study, charging contributions to binding free energies of small oligoatomic ions to a series of model host cavities functionalized with different chemical groups are calculated with classical atomistic molecular dynamics simulation. Electrostatic interactions are treated using a lattice-summation scheme or a cutoff-truncation scheme with Barker-Watts reaction-field correction, and the simulations are conducted in boxes of different edge lengths. It is illustrated that the charging free energies of the guest molecules in water and in the host strongly depend on the applied methodology and that neglect of correction terms for the artifacts introduced by the finite size of the simulated system and the use of an effective electrostatic interaction function considerably impairs the thermodynamic interpretation of guest-host interactions. Application of correction terms for the various artifacts yields consistent results for the charging contribution to binding free energies and is thus a prerequisite for the valid interpretation or prediction of experimental data via molecular dynamics simulation. Analysis and correction of electrostatic artifacts according to the scheme proposed in the present study should therefore be considered an integral part of careful free-energy calculation studies if changes in the net charge are involved.

State of hydrophobic and hydrophilic ionic liquids in aqueous solutions: are the ions fully dissociated?

Molecular dynamics simulations were performed for aqueous solutions of five ionic liquids (ILs): 1-ethyl-3-methylimidazolium ([C2mim]) bis(trifluoromethanesulfonyl) imide ([NTf2]), 1-n-butyl-3-methylimidazolium ([C4mim]) [NTf2], 1-n-hexyl-3-methylimidazolium ([C6mim]) [NTf2], [C2mim] ethylsulfate ([C2H5SO4]), and [C2mim] chloride (Cl) to determine whether the ions of these ILs are associated at relatively high dilutions and whether the association is governed by hydrophobicity/hydrophilicity of the ILs. The adaptive biasing force technique was applied to calculate the potential of mean force (PMF) for each IL ion pair. For all of the ILs, the PMF is characterized by two distinct contact minima in which the ions have different relative conformations. The hydrophobic ILs bearing the anion [NTf2](-) exist predominantly in the associative state; the strength of the association of these ILs increases with increase in the alkyl chain length. The most hydrophilic IL [C2mim] Cl was determined to be almost fully dissociated at the concentration examined in the study. [C2mim] [C2H5SO4] showed hydration behavior that was intermediate between that exhibited by the ILs in which the anion is substituted with either Cl(-) or [NTf2](-) paired with [C2mim](+). Association constants for these ILs were also computed. Radial distribution functions calculated by constraining the ions at the contact minima showed that hydration of the anion plays the dominant role in determining the microscopic behavior of these ILs in aqueous solutions.

Hydrated polyamide membrane and its interaction with alginate: a molecular dynamics study

The properties of the hydrated amorphous polyamide (PA) membrane and its binding with alginate are investigated through molecular dynamics simulations. The density of the hydrated membrane, surface morphology, and water diffusion near and inside the membrane are compared to other studies. Particular focus is given to the steered molecular dynamics (SMD) simulation of the binding between the PA membrane and an alginate model. The PA surface composition is determined on the basis of experimental measurements of the oxygen/nitrogen (O/N) ratio. The surface model is built using a configurational-bias Monte Carlo technique. The consistent valence force field (CVFF) is used to describe the atomic interactions in the membrane-foulant system. Simulation results show that the carboxylate groups in both the PA surface and alginate exhibit strong binding with metal ions. This binding mechanism plays a major role in the PA-alginate fouling through the formation of an ionic binding bridge. Specifically, Ca(2+) ions have stronger binding with the carboxylate group than Na(+) ions, while the binding breakdown time is shorter for Ca(2+) than Na(+) because of the comparably higher hydration free energy of Ca(2+) ions with water molecules.

Stability Mechanisms of Laccase Isoforms using a Modified FoldX Protocol Applicable to Widely Different Proteins

A recent computational protocol that accurately predicts and rationalizes protein multisite mutant stabilities has been extended to handle widely different isoforms of laccases. We apply the protocol to four isoenzymes of Trametes versicolor laccase (TvL) with variable lengths (498-503 residues) and thermostability (Topt ∼ 45-80 °C) and with 67-77% sequence identity. The extended protocol uses (i) statistical averaging, (ii) a molecular-dynamics-validated "compromise" homology model to minimize bias that causes proteins close in sequence to a structural template to be too stable due to having the benefits of the better sampled template (typically from a crystal structure), (iii) correction for hysteresis that favors the input template to overdestabilize, and (iv) a preparative protocol to provide robust input sequences of equal length. The computed ΔΔG values are in good agreement with the major trends in experimental stabilities; that is, the approach may be applicable for fast estimates of the relative stabilities of proteins with as little as 70% identity, something that is currently extremely challenging. The computed stability changes associated with variations are Gaussian-distributed, in good agreement with experimental distributions of stability effects from mutation. The residues causing the differential stability of the four isoforms are consistent with a range of compiled laccase wild type data, suggesting that we may have identified general drivers of laccase stability. Several sites near Cu, notably 79, 241, and 245, or near substrate, mainly 265, are identified that contribute to stability-function trade-offs, of relevance to the search for new proficient and stable variants of these important industrial enzymes.

Stability mechanisms of a thermophilic laccase probed by molecular dynamics

Laccases are highly stable, industrially important enzymes capable of oxidizing a large range of substrates. Causes for their stability are, as for other proteins, poorly understood. In this work, multiple-seed molecular dynamics (MD) was applied to a Trametes versicolor laccase in response to variable ionic strengths, temperatures, and glycosylation status. Near-physiological conditions provided excellent agreement with the crystal structure (average RMSD ∼0.92 Å) and residual agreement with experimental B-factors. The persistence of backbone hydrogen bonds was identified as a key descriptor of structural response to environment, whereas solvent-accessibility, radius of gyration, and fluctuations were only locally relevant. Backbone hydrogen bonds decreased systematically with temperature in all simulations (∼9 per 50 K), probing structural changes associated with enthalpy-entropy compensation. Approaching T_opt (∼350 K) from 300 K, this change correlated with a beginning “unzipping” of critical β-sheets. 0 M ionic strength triggered partial denucleation of the C-terminal (known experimentally to be sensitive) at 400 K, suggesting a general salt stabilization effect. In contrast, F⁻ (but not Cl⁻) specifically impaired secondary structure by formation of strong hydrogen bonds with backbone NH, providing a mechanism for experimentally observed small anion destabilization, potentially remedied by site-directed mutagenesis at critical intrusion sites. N-glycosylation was found to support structural integrity by increasing persistent backbone hydrogen bonds by ∼4 across simulations, mainly via prevention of F⁻ intrusion. Hydrogen-bond loss in distinct loop regions and ends of critical β-sheets suggest potential strategies for laboratory optimization of these industrially important enzymes.

Ab initio modeling and experimental assessment of Janus Kinase 2 (JAK2) kinase-pseudokinase complex structure

The Janus Kinase 2 (JAK2) plays essential roles in transmitting signals from multiple cytokine receptors, and constitutive activation of JAK2 results in hematopoietic disorders and oncogenesis. JAK2 kinase activity is negatively regulated by its pseudokinase domain (JH2), where the gain-of-function mutation V617F that causes myeloproliferative neoplasms resides. In the absence of a crystal structure of full-length JAK2, how JH2 inhibits the kinase domain (JH1), and how V617F hyperactivates JAK2 remain elusive. We modeled the JAK2 JH1-JH2 complex structure using a novel informatics-guided protein-protein docking strategy. A detailed JAK2 JH2-mediated auto-inhibition mechanism is proposed, where JH2 traps the activation loop of JH1 in an inactive conformation and blocks the movement of kinase αC helix through critical hydrophobic contacts and extensive electrostatic interactions. These stabilizing interactions are less favorable in JAK2-V617F. Notably, several predicted binding interfacial residues in JH2 were confirmed to hyperactivate JAK2 kinase activity in site-directed mutagenesis and BaF3/EpoR cell transformation studies. Although there may exist other JH2-mediated mechanisms to control JH1, our JH1-JH2 structural model represents a verifiable working hypothesis for further experimental studies to elucidate the role of JH2 in regulating JAK2 in both normal and pathological settings.

Atomistic Simulations Predict a Surprising Variety of Morphologies in Precise Ionomers

The nature of ionic aggregates in ionomers remains an important open question, particularly considering its significance to their unique electrical and mechanical properties. We have carried out fully atomistic molecular dynamics simulations of melts of lithium-neutralized precise ionomers that reveal the structural features of ionic aggregates in unprecedented detail. In particular, we observe a rich variety of aggregate morphologies depending on neutralization level and ionic content, including string-like and percolated aggregates. The traditional assumption of spherical ionic aggregates with liquid-like ordering that is typically used to interpret experimental scattering data is too simplistic; a more rich and complex set of structures exist that also fit the scattering data.

Exploring Adsorption of Water and Ions on Carbon Surfaces using a Polarizable Force Field

Graphene, carbon nanotubes, and fullerenes are of great interest due to their unique properties and diverse applications in biology, molecular electronics, and materials science. Therefore, there is demand for methods that can accurately model the interface between carbon surfaces and their environment. In this letter we compare results for complexes of water, potassium ion, and chloride ion with graphene, carbon nanotube, and fullerene surfaces using a standard non-polarizable force field (OPLS-AA), a polarizable force field (OPLS-AAP), DFT, and ab initio theory. For interactions with water, OPLS-AA with the TIP3P or TIP4P water models describes the interactions with benzene (C(6)H(6)) and coronene (C(24)H(12)) well; however, for acenes larger than circumcoronene (C(54)H(18)) and especially for C(60), the interaction energies are somewhat too weak and polarization is needed. For ions interacting with carbon surfaces, inclusion of polarization is essential, and OPLS-AAP is found to perform well in comparison to the highest-level quantum mechanical methods. Overall, OPLS-AAP provides an accurate and computationally efficient force field for modeling condensed-phase systems featuring carbon surfaces.