PACKMOL: a package for building initial configurations for molecular dynamics simulations

Hydration of calcium oxide surface predicted by reactive force field molecular dynamics

In this work, we present the parametrization of Ca-O/H interactions within the reactive force field ReaxFF, and its application to study the hydration of calcium oxide surface. The force field has been fitted using density functional theory calculations on gas phase calcium-water clusters, calcium oxide bulk and surface properties, calcium hydroxide, bcc and fcc Ca, and proton transfer reactions in the presence of calcium. Then, the reactive force field has been used to study the hydration of the calcium oxide {001} surface with different water contents. Calcium oxide is used as a catalyzer in many applications such as CO(2) sequestration and biodiesel production, and the degree of surface hydroxylation is a key factor in its catalytic performance. The results show that the water dissociates very fast on CaO {001} bare surfaces without any defect or vacancy. The surface structure is maintained up to a certain amount of water, after which the surface undergoes a structural rearrangement, becoming a disordered calcium hydroxyl layer. This transformation is the most probable reason for the CaO catalytic activity decrease.

Kinetics of the terminal electron transfer step in cytochrome c oxidase

Cytochrome c oxidase (cco) catalyzes the oxygen reduction reaction in most aerobically respiring organisms. Decades of research have uncovered many aspects relating to structure and function of this enzyme. However, the origin of the unusually fast terminal electron transfer step from heme a to heme a(3) in cco has been the subject of intense discussions over recent years. Yet, no satisfactory consensus has been achieved. Carrying out large-scale molecular dynamics simulation of the protein embedded in a solvated membrane, we obtain a reorganization free energy λ = 0.57 eV. Evaluation of the quantized single-mode rate equation using the experimental rate and the computed reorganization free energy gives a value of 1.5 meV for the average electronic coupling (H(ab)) between heme a and heme a(3). Thus, according to our calculations, the nanosecond electron transfer (ET) is due to a small but significant activation barrier (ΔG(‡) = 0.12 eV) in combination with effective electronic coupling between the two cofactors. The activation free energy is caused predominantly by collective reorganization of protein residues. We show that our results are consistent with the weak temperature dependence observed in experiment if one allows for very minor variations in the donor-acceptor distance as the temperature changes.

Self-assembly of DNA on a gapped carbon nanotube

We perform molecular dynamics simulations to analyze the wrapping process of a single-stranded (ss) DNA around a gapped CNT immersed in a bath of water. We observe the formation of a stable molecular junction with the ssDNA adopting a helical or circular conformation around one CNT electrode and a linear conformation around the opposite electrode. We find that DNA undergoes several conformational changes during equilibration of the self-assembled molecular junction. This process would allow a higher yield of successful CNT-DNA interconnections, which constitutes a novel structure of interest in chemical and biological sensing at the single-molecule level.

Hb S-São Paulo: a new sickling hemoglobin with stable polymers and decreased oxygen affinity

Hb S-São Paulo (SP) [HBB:c.20A>T p.Glu6Val; c.196A>G p.Lys65Glu] is a new double-mutant hemoglobin that was found in heterozygosis in an 18-month-old Brazilian male with moderate anemia. It behaves like Hb S in acid electrophoresis, isoelectric focusing and solubility testing but shows different behavior in alkaline electrophoresis, cation-exchange HPLC and RP-HPLC. The variant is slightly unstable, showed reduced oxygen affinity and also appeared to form polymers more stable than the Hb S. Molecular dynamics simulation suggests that the polymerization is favored by interfacial electrostatic interactions. This provides a plausible explanation for some of the reported experimental observations.

Separation of chemical reaction intermediates by metal-organic frameworks

HPLC columns custom-packed with metal-organic framework (MOF) materials are used for the separation of four small intermediates and byproducts found in the commercial synthesis of an important active pharmaceutical ingredient in methanol. In particular, two closely related amines can be separated in the methanol reaction medium using MOFs, but not with traditional C18 columns using an optimized aqueous mobile phase. Infrared spectroscopy, UV-vis spectroscopy, X-ray diffraction, and thermogravimetric analysis are used in combination with molecular dynamic simulations to study the separation mechanism for the best-performing MOF materials. It is found that separation with ZIF-8 is the result of an interplay between the thermodynamic driving force for solute adsorption within the framework pores and the kinetics of solute diffusion into the material pores, while the separation with Basolite F300 is achieved because of the specific interactions between the solutes and Fe(3+) sites. This work, and the exceptional ability to tailor the porous properties of MOF materials, points to prospects for using MOF materials for the continuous separation and synthesis of pharmaceutical compounds.

A refined potential for hydroxylamine clusters and the liquid phase

A detailed study including ab initio calculations and classic Monte-Carlo simulations of hydroxylamine in the gas and liquid phases is presented. A classical interaction potential for hydroxylamine, which includes polarizability, many-body effects, and intramolecular relaxation, was constructed. The results of the simulation were compared to the available experimental data in order to validate the model. We conclude that liquid hydroxylamine has a multitude of hydrogen bonds leading to a large density where the existence of cis conformers and clusters of these conformers is possible. This explains the occurrence of the classical [R. Nast and I. Z. Foppl, Z. Anorg. Allg. Chem. 263, 310 (1950)] scheme for the molecule's decomposition at room temperature and its large exothermicity and instability.

Switch of substrate specificity of hyperthermophilic acylaminoacyl peptidase by combination of protein and solvent engineering

The inherent evolvability of promiscuous enzymes endows them with great potential to be artificially evolved for novel functions. Previously, we succeeded in transforming a promiscuous acylaminoacyl peptidase (apAAP) from the hyperthermophilic archaeon Aeropyrum pernix K1 into a specific carboxylesterase by making a single mutation. In order to fulfill the urgent requirement of thermostable lipolytic enzymes, in this paper we describe how the substrate preference of apAAP can be further changed from p-nitrophenyl caprylate (pNP-C8) to p-nitrophenyl laurate (pNP-C12) by protein and solvent engineering. After one round of directed evolution and subsequent saturation mutagenesis at selected residues in the active site, three variants with enhanced activity towards pNP-C12 were identified. Additionally, a combined mutant W474V/F488G/R526V/T560W was generated, which had the highest catalytic efficiency (k (cat)/K (m)) for pNP-C12, about 71-fold higher than the wild type. Its activity was further increased by solvent engineering, resulting in an activity enhancement of 280-fold compared with the wild type in the presence of 30% DMSO. The structural basis for the improved activity was studied by substrate docking and molecular dynamics simulation. It was revealed that W474V and F488G mutations caused a significant change in the geometry of the active center, which may facilitate binding and subsequent hydrolysis of bulky substrates. In conclusion, the combination of protein and solvent engineering may be an effective approach to improve the activities of promiscuous enzymes and could be used to create naturally rare hyperthermophilic enzymes.

Coarse-grained molecular dynamics simulations of the sphere to rod transition in surfactant micelles

Surfactant molecules self-assemble in aqueous solutions to form various micellar structures such as spheres, rods, or lamellae. Although phase transitions in surfactant solutions have been studied experimentally, their molecular mechanisms are still not well understood. In this work, we show that molecular dynamics (MD) simulations using the coarse-grained (CG) MARTINI force field and explicit CG solvent, validated against atomistic MD studies, can accurately represent micellar assemblies of cetyltrimethylammonium chloride (CTAC). The effect of salt on micellar structures is studied for aromatic anionic salts, e.g., sodium salicylate, and simple inorganic salts, e.g., sodium chloride. Above a threshold concentration, sodium salicylate induces a sphere to rod transition in the micelle. CG MD simulations are shown to capture the dynamics of this shape transition and support a mechanism based on the reduction in the micelle-water interfacial tension induced by the adsorption of the amphiphilic salicylate ions. At the threshold salt concentration, the interface is nearly saturated with adsorbed salicylate ions. Predictions of the effect of salt on the micelle structure in different CG solvent models, namely, single-site standard water and three-site polarizable water, show qualitative agreement. This suggests that phase transitions in aqueous micelle solutions could be investigated by using standard CG water models which allow for 3 orders of magnitude reduction in the computational time as compared to that required for atomistic MD simulations.

Capturing spontaneous partitioning of peripheral proteins using a biphasic membrane-mimetic model

Membrane binding of peripheral proteins, mediated by specialized anchoring domains, is a crucial step for their biological function. Computational studies of membrane insertion, however, have proven challenging and largely inaccessible, due to the time scales required for the complete description of the process, mainly caused by the slow diffusion of the lipid molecules composing the membrane. Furthermore, in many cases, the nature of the membrane "anchor", i.e., the part of the protein that inserts into the membrane, is also unknown. Here, we address some of these issues by developing and employing a simplified representation of the membrane by a biphasic solvent model which we demonstrate can be used efficiently to capture and describe the process of hydrophobic insertion of membrane anchoring domains in all-atom molecular dynamics simulations. Applying the model, we have studied the insertion of the anchoring domain of a coagulation protein (the GLA domain of human protein C), starting from multiple initial configurations varying with regard to the initial orientation and height of the protein with respect to the membrane. In addition to efficiently and consistently identifying the "keel" region as the hydrophobic membrane anchor, within a few nanoseconds each configuration simulated showed a convergent height (2.20 ± 1.04 Å) and angle with respect to the interface normal (23.37 ± 12.48°). We demonstrate that the model can produce the same results as those obtained from a full representation of a membrane, in terms of both the depth of penetration and the orientation of the protein in the final membrane-bound form with an order of magnitude decrease in the required computational time compared to previous models, allowing for a more exhaustive search for the correct membrane-bound configuration.

Helix 12 dynamics and thyroid hormone receptor activity: experimental and molecular dynamics studies of Ile280 mutants

Nuclear hormone receptors (NRs) form a family of transcription factors that mediate cellular responses initiated by hormone binding. It is generally recognized that the structure and dynamics of the C-terminal helix 12 (H12) of NRs' ligand binding domain (LBD) are fundamental to the recognition of coactivators and corepressors that modulate receptor function. Here we study the role of three mutations in the I280 residue of H12 of thyroid hormone receptors using site-directed mutagenesis, functional assays, and molecular dynamics simulations. Although residues at position 280 do not interact with coactivators or with the ligand, we show that its mutations can selectively block coactivator and corepressor binding, and affect hormone binding affinity differently. Molecular dynamics simulations suggest that ligand affinity is reduced by indirectly displacing the ligand in the binding pocket, facilitating water penetration and ligand destabilization. Mutations I280R and I280K link H12 to the LBD by forming salt bridges with E457 in H12, stabilizing H12 in a conformation that blocks both corepressor and coactivator recruitment. The I280M mutation, in turn, blocks corepressor binding, but appears to enhance coactivator affinity, suggesting stabilization of H12 in agonist conformation.

Mechanism of reactant and product dissociation from the anthrax edema factor: a locally enhanced sampling and steered molecular dynamics study

The anthrax edema factor is a toxin overproducing damaging levels of cyclic adenosine monophosphate (cAMP) and pyrophosphate (PPi) from ATP. Here, mechanisms of dissociation of ATP and products (cAMP, PPi) from the active site are studied using locally enhanced sampling (LES) and steered molecular dynamics simulations. Various substrate conformations and ionic binding modes found in crystallographic structures are considered. LES simulations show that PPi and cAMP dissociate through different solvent accessible channels, while ATP dissociation requires significant active site exposure to solvent. The ionic content of the active site directly affects the dissociation of ATP and products. Only one ion dissociates along with ATP in the two-Mg(2+) binding site, suggesting that the other ion binds EF prior to ATP association. Dissociation of reaction products cAMP and PPi is impaired by direct electrostatic interactions between products and Mg(2+) ions. This provides an explanation for the inhibitory effect of high Mg(2+) concentrations on EF enzymatic activity. Breaking of electrostatic interactions is dependent on a competitive binding of water molecules to the ions, and thus on the solvent accessibility of the active site. Consequently, product dissociation seems to be a two-step process. First, ligands are progressively solvated while preserving the most important electrostatic interactions, in a process that is dependent on the flexibility of the active site. Second, breakage of the electrostatic bonds follows, and ligands diffuse into solvent. In agreement with this mechanism, product protonation facilitates dissociation.

Local and bulk hydration of zwitterionic glycine and its analogues through molecular simulations

Molecular dynamics simulations were used to characterize the differences in hydration between glycine and two of its zwitterionic analogues: N,N-dimethylglycine and N,N,N-trimethylglycine (glycine betaine). The hydration of dodecane and oligo(ethylene glycol) was studied for reference. Both structuring and dynamics of bulk and bound water were examined using a variety of properties and at multiple concentrations. Metrics, such as radial distribution functions and residence times, were used to characterize hydration. Also, we used more specialized metrics that can discriminate between subtle differences in hydration, such as condensed phase order parameters, Voronoi tessellations, and multidimensional pair-pair correlation functions. Trimethylglycine was found to have a unique hydration shell that extends across the entire molecule and has no specific interactions between solute molecules. Also, we found that dimethylglycine has a similar hydration structure to that of trimethylglycine despite its hydrogen-bond donor. Glycine was found to aggregate and have a more disjoint hydration shell. All three zwitterions were found to structurally affect water within 1.5-2.0 coordination shells. Lastly, trimethylglycine is disperse in solution even at very high concentrations, and water rapidly moves between trimethylglycine amine groups. This work has meaningful implications for protein stability where trimethylglycine is known to prevent protein aggregation and nonfouling interfaces where trimethylglycine prevents nonspecific protein adsorption.

Analysis of agonist and antagonist effects on thyroid hormone receptor conformation by hydrogen/deuterium exchange

Thyroid hormone receptors (TRs) are ligand-gated transcription factors with critical roles in development and metabolism. Although x-ray structures of TR ligand-binding domains (LBDs) with agonists are available, comparable structures without ligand (apo-TR) or with antagonists are not. It remains important to understand apo-LBD conformation and the way that it rearranges with ligands to develop better TR pharmaceuticals. In this study, we conducted hydrogen/deuterium exchange on TR LBDs with or without agonist (T(3)) or antagonist (NH3). Both ligands reduce deuterium incorporation into LBD amide hydrogens, implying tighter overall folding of the domain. As predicted, mass spectroscopic analysis of individual proteolytic peptides after hydrogen/deuterium exchange reveals that ligand increases the degree of solvent protection of regions close to the buried ligand-binding pocket. However, there is also extensive ligand protection of other regions, including the dimer surface at H10-H11, providing evidence for allosteric communication between the ligand-binding pocket and distant interaction surfaces. Surprisingly, C-terminal activation helix H12, which is known to alter position with ligand, remains relatively protected from solvent in all conditions suggesting that it is packed against the LBD irrespective of the presence or type of ligand. T(3), but not NH3, increases accessibility of the upper part of H3-H5 to solvent, and we propose that TR H12 interacts with this region in apo-TR and that this interaction is blocked by T(3) but not NH3. We present data from site-directed mutagenesis experiments and molecular dynamics simulations that lend support to this structural model of apo-TR and its ligand-dependent conformational changes.

Small-angle scattering and the structure of ambient liquid water

Structural polyamorphism has been promoted as a means for understanding the anomalous thermodynamics and dynamics of water in the experimentally inaccessible supercooled region. In the metastable liquid region, theory has hypothesized the existence of a liquid-liquid critical point from which a dividing line separates two water species of high and low density. A recent small-angle X-ray scattering study has claimed that the two structural species postulated in the supercooled state are seen to exist in bulk water at ambient conditions. We analyze new small-angle X-ray scattering data on ambient liquid water taken at third generation synchrotron sources, and large 32,000 water molecule simulations using the TIP4P-Ew model of water, to show that the small-angle region measures standard number density fluctuations consistent with water's isothermal compressibility temperature trends. Our study shows that there is no support or need for heterogeneities in water structure at room temperature to explain the small-angle scattering data, as it is consistent with a unimodal density of the tetrahedral liquid at ambient conditions.

ATP conformations and ion binding modes in the active site of anthrax edema factor: a computational analysis

The Edema Factor (EF), one of the virulence factors of anthrax, is an adenylyl cyclase that promotes the overproduction of cyclic-AMP (cAMP) from ATP, and therefore perturbs cell signaling. Crystallographic structures of EF bound to ATP analogs and reaction products, cyclic-AMP, and Pyrophosphate (PPi), revealed different substrate conformations and catalytic-cation binding modes, one or two cations being observed in the active site. To shed light into the biological significance of these crystallographic structures, the energetics, geometry, and dynamics of the active site are analyzed using molecular dynamics simulations. The ATP conformation observed in the one-metal-ion structure allows stronger interactions with the catalytic ion, and ATP is more restrained than in the structure containing two Mg(2+) ions. Therefore, we propose that the conformation observed in the one-ion crystal structure is a more probable starting point for the reaction. The simulations also suggest that a C3'-endo sugar pucker facilitates nucleophilic attack. Additionally, the two-cation binding mode restrains the mobility of the reaction products, and thus their tendency to dissociate.