Bruno A. C. Horta

GROMOS++ Software for the Analysis of Biomolecular Simulation Trajectories

GROMOS++ is a set of C++ programs for pre- and postprocessing of molecular dynamics simulation trajectories and as such is part of the GROningen MOlecular Simulation software for (bio)molecular simulation. It contains more than 70 programs that can be used to prepare data for the production of molecular simulation trajectories and to analyze these. These programs are reviewed and the various structural, dynamic, and thermodynamic quantities that can be analyzed using time series, correlation functions, and distributions are described together with technical aspects of their implementation in GROMOS. A few examples of the use of GROMOS++ for the analysis of MD trajectories are given. A full list of all GROMOS++ programs, together with an indication of their capabilities, is given in the Appendix .

New functionalities in the GROMOS biomolecular simulation software

Since the most recent description of the functionalities of the GROMOS software for biomolecular simulation in 2005 many new functions have been implemented. In this article, the new functionalities that involve modified forces in a molecular dynamics (MD) simulation are described: the treatment of electronic polarizability, an implicit surface area and internal volume solvation term to calculate interatomic forces, functions for the GROMOS coarse‐grained supramolecular force field, a multiplicative switching function for nonbonded interactions, adiabatic decoupling of a number of degrees of freedom with temperature or force scaling to enhance sampling, and nonequilibrium MD to calculate the dielectric permittivity or viscosity. Examples that illustrate the use of these functionalities are given.

MKTOP: a program for automatic construction of molecular topologies

Molecular mechanics methods are widely-used for simulations of biomolecules. A molecular topology file containing all atom types, which depend on the chemical environment, must be constructed. MKTOP is a free-software capable of automatic atom type recognition and construction of molecular topologies for GROMACS, being a useful tool for its users.

New Interaction Parameters for Oxygen Compounds in the GROMOS Force Field: Improved Pure-Liquid and Solvation Properties for Alcohols, Ethers, Aldehydes, Ketones, Carboxylic Acids, and Esters

A new parameter set (53A6OXY) is developed for the GROMOS force field, that combines reoptimized parameters for the oxygen-containing chemical functions (alcohols, ethers, aldehydes, ketones, carboxylic acids, and esters) with the current biomolecular force field version (53A6) for all other functions. In the context of oxygen-containing functions, the 53A6OXY parameter set is obtained by optimization of simulated pure-liquid properties, namely the density ρliq and enthalpy of vaporization ΔHvap, as well as solvation properties, namely the free energies of solvation in water ΔGwat and in cyclohexane ΔGche, against experimental data for 10 selected organic compounds, and further tested for 25 other compounds. The simultaneous refinement of atomic charges and Lennard-Jones interaction parameters against the four mentioned types of properties provides a single parameter set for the simulation of both liquid and biomolecular systems. Small changes in the covalent parameters controlling the geometry of the oxygen-containing chemical functions are also undertaken. The new 53A6OXY force-field parameters reproduce the mentioned experimental data within root-mean-square deviations of 22.4 kg m(-3) (ρliq), 3.1 kJ mol(-1) (ΔHvap), 3.0 kJ mol(-1) (ΔGwat), and 1.7 kJ mol(-1) (ΔGche) for the 35 compounds considered.

Interaction of the disaccharides trehalose and gentiobiose with lipid bilayers: A comparative molecular dynamics study

The disaccharide α,α-trehalose (TRH) is known for its bioprotective action in organisms subject to stressful environmental conditions. However, the mechanisms whereby TRH stabilizes biomolecules remains matter of debate, the five main hypotheses being the water replacement (WRH), headgroup bridging (HBH), vitrification (VIH), water entrapment (WEH) and hydration forces (HFH) hypotheses. Four hypotheses (all except HFH) are in principle compatible with a preferential affinity of the sugar molecules (compared to water) for the biomolecular surface. According to the recently proposed sugar-like mechanism (Pereira and Hünenberger), preferential affinity would result from the entropy gain upon releasing many water molecules from the surface region to the bulk, at the cost of immobilizing and rigidifying fewer sugar molecules. Thus, a more flexible disaccharide such as gentiobiose (GNT) should evidence a weaker preferential affinity, limiting its bioprotective ability. In this work, molecular dynamics (MD) simulations of a dipalmitoyl-phosphatidylcholine (DPPC) bilayer patch in the presence of either pure water or aqueous solutions of GNT or TRH are performed in order to assess the validity of this suggestion. At 475 K and 1.6 m (molal), TRH indeed preserves the bilayer structure to a larger extent compared to GNT. However, the present investigation does not unambiguously indicate which of the above mechanism takes place, since the simulations reveal characteristic features of all of them. This suggests either that multiple mechanisms may be simultaneously active or that their definitions are not precise enough.

Engineered monomeric human histidine triad nucleotide-binding protein 1 hydrolyzes fluorogenic acyl-adenylate and lysyl-tRNA synthetase-generated lysyl-adenylate

Hint1 is a homodimeric protein and member of the ubiquitous HIT superfamily. Hint1 catalyzes the hydrolysis of purine phosphoramidates and lysyl-adenylate generated by lysyl-tRNA synthetase (LysRS). To determine the importance of homodimerization on the biological and catalytic activity of Hint1, the dimer interface of human Hint1 (hHint1) was destabilized by replacement of Val97 of hHint1 with Asp, Glu, or Arg. The mutants were shown to exist as monomers in solution by a combination of size exclusion chromatograph, static light scattering, and chemically induced dimerization studies. Circular dichroism studies revealed little difference between the stability of the V97D, V97E, and wild-type hHint1. Relative to wild-type and the V97E mutant, however, significant perturbation of the V97D mutant structure was observed. hHint1 was shown to prefer 3-indolepropionic acyl-adenylate (AIPA) over tryptamine adenosine phosphoramidate monoester (TpAd). Wild-type hHint1 was found to be 277- and 1000-fold more efficient (kcat/Km values) than the V97E and V97D mutants, respectively. Adenylation of wild-type, V97D, and V97E hHint1 by human LysRS was shown to correlate with the mutant kcat/Km values using 3-indolepropionic acyl-adenylate as a substrate, but not tryptamine adenosine phosphoramidate monoester. Significant perturbations of the active site residues were not detected by molecular dynamics simulations of the hHint1s. Taken together, these results demonstrate that for hHint1; 1) the efficiency (kcat/Km) of acylated AMP hydrolysis, but not maximal catalytic turnover (kcat), is dependent on homodimerization and 2) the hydrolysis of lysyl-AMP generated by LysRS is not dependent on homodimerization if the monomer structure is similar to the wild-type structure.

Simulating the Transition between Gel and Liquid-Crystal Phases of Lipid Bilayers: Dependence of the Transition Temperature on the Hydration Level

Explicit-solvent molecular dynamics (MD) simulations of the monoglyceride glycerol-1-monopalmitin (GMP; bilayer patch of 2 × 6 × 6 lipids) at different hydration levels (full, half, or quarter hydration) and at different temperatures (318 to 338 K) are reported. The 40 ns simulations (some extended to 200 ns) are initiated from structures appropriate for the gel (GL) or liquid-crystal (LC) phases, with the goal of investigating whether atomistic MD simulations on this time scale can be used to monitor GL ↔ LC transitions in lipid bilayers, and to evaluate the corresponding transition temperatures Tm (as well as the influence of hydration on Tm) in a reliable fashion. The main conclusions are as follows: (i) The GL → LC transition of GMP can be observed on the 40 ns time scale. (ii) The LC → GL transition is comparatively slower, requiring simulations on the 200 ns time scale. (iii) A set of simulations initiated from a structure appropriate for the GL phase and carried out at slightly different temperatures permits the determination of a reliable value for Tm. (iv) The calculated Tm values reproduce the experimentally observed increase in this transition temperature upon decreasing the bilayer hydration. (v) The Tm values calculated at the three hydration levels considered are in essentially quantitative agreement with the experimental phase diagram of GMP. To our knowledge, this study represents the first accurate determination of the Tm of a lipid via atomistic simulations of the (reversible) GL ↔ LC phase transition, as well as the first direct simulation evidence for the increase in the transition temperature upon dehydration. The possible direct determination of Tm and the characterization of environmental effects on this quantity by simulation opens up promising perspectives in the contexts of force field refinement and the investigation of dehydration-induced damages in living cells (and bioprotection by cosolutes).

Association of the anti-tuberculosis drug rifampicin with a PAMAM dendrimer

The association of the anti-tuberculosis drug rifampicin (RIF) with a 4th-generation poly(amidoamine) (G4-PAMAM) dendrimer was investigated by means of molecular dynamics simulations. The RIF load capacity was estimated to be around 20 RIF per G4-PAMAM at neutral pH. The complex formed by 20 RIF molecules and the dendrimer (RIF20-PAMAM) was subjected to 100 ns molecular dynamics (MD) simulations at two different pH conditions (neutral and acidic). The complex was found to be significantly more stable in the simulation at neutral pH compared to the simulation at low pH in which the RIF molecules were rapidly and almost simultaneously expelled to the solvent bulk. The high stability of the RIF-PAMAM complex under physiological pH and the rapid release of RIF molecules under acidic medium provide an interesting switch for drug targeting since the Mycobacterium resides within acidic domains of the macrophage. Altogether, these results suggest that, at least in terms of stability and pH-dependent release, PAMAM-like dendrimers may be considered suitable drug delivery systems for RIF and derivatives.

Origin of the spectral shifts among the early intermediates of the rhodopsin photocycle.

A combined strategy based on the computation of absorption energies, using the ZINDO/S semiempirical method, for a statistically relevant number of thermally sampled configurations extracted from QM/MM trajectories is used to establish a one-to-one correspondence between the structures of the different early intermediates (dark, batho, BSI, lumi) involved in the initial steps of the rhodopsin photoactivation mechanism and their optical spectra. A systematic analysis of the results based on a correlation-based feature selection algorithm shows that the origin of the color shifts among these intermediates can be mainly ascribed to alterations in intrinsic properties of the chromophore structure, which are tuned by several residues located in the protein binding pocket. In addition to the expected electrostatic and dipolar effects caused by the charged residues (Glu113, Glu181) and to strong hydrogen bonding with Glu113, other interactions such as π-stacking with Ala117 and Thr118 backbone atoms, van der Waals contacts with Gly114 and Ala292, and CH/π weak interactions with Tyr268, Ala117, Thr118, and Ser186 side chains are found to make non-negligible contributions to the modulation of the color tuning among the different rhodopsin photointermediates.

A GROMOS-Compatible Force Field for Small Organic Molecules in the Condensed Phase: The 2016H66 Parameter Set.

This article reports on the calibration and validation of a new GROMOS-compatible parameter set 2016H66 for small organic molecules in the condensed phase. The calibration is based on 62 organic molecules spanning the chemical functions alcohol, ether, aldehyde, ketone, carboxylic acid, ester, amine, amide, thiol, sulfide, and disulfide, as well as aromatic compounds and nucleic-acid bases. For 57 organic compounds, the calibration targets are the experimental pure-liquid density ρliq and the vaporization enthalpy ΔHvap, as well as the hydration free energy ΔGwat and the solvation free energy ΔGche in cyclohexane, at atmospheric pressure and at (or close to) room temperature. The final root-mean-square deviations (RMSD) for these four quantities over the set of compounds are 32.4 kg m(-3), 3.5 kJ mol(-1), 4.1 kJ mol(-1), and 2.1 kJ mol(-1), respectively, and the corresponding average deviations (AVED) are 1.0 kg m(-3), 0.2 kJ mol(-1), 2.6 kJ mol(-1), and 1.0 kJ mol(-1), respectively. For the five nucleic-acid bases, the parametrization is performed by transferring the final 2016H66 parameters from analogous organic compounds followed by a slight readjustment of the charges to reproduce the experimental water-to-chloroform transfer free energies ΔGtrn. The final RMSD for this quantity over the five bases is 1.7 kJ mol(-1), and the corresponding AVED is 0.8 kJ mol(-1). As an initial validation of the 2016H66 set, seven additional thermodynamic, transport, and dielectric properties are calculated for the 57 organic compounds in the liquid phase. The agreement with experiment in terms of these additional properties is found to be reasonable, with significant deviations typically affecting either a specific chemical function or a specific molecule. This suggests that in most cases, a classical force-field description along with a careful parametrization against ρliq, ΔHvap, ΔGwat, and ΔGche results in a model that appropriately describes the liquid in terms of a wide spectrum of its physical properties.