A method to apply bond-angle constraints in molecular dynamics simulations.

An algorithm to apply bond-angle constraints in molecular dynamics simulations of macromolecules or molecular liquids is presented. It uses Cartesian coordinates and determines the Lagrange multipliers required for maintaining the constraints iteratively. It constitutes an alternative to the use of only distance constraints (DCs) between particles to maintain a particular geometry. DCs are unsuitable to maintain particular, for example, linear or flat, geometries of molecules. The proposed algorithm can easily handle bond-length, bond-angle, and dihedral-angle constraints simultaneously, as when calculating a potential of mean force along a dihedral-angle degree of freedom.

Rotational Profiler: A Fast, Automated, and Interactive Server to Derive Torsional Dihedral Potentials for Classical Molecular Simulations.

Rotational Profiler provides an analytical algorithm to compute sets of classical torsional dihedral parameters by fitting an empirical energy profile to a reference one that can be obtained experimentally or by quantum-mechanical methods. The resulting profiles are compatible with the functional forms in the most widely used biomolecular force fields (e.g., GROMOS, AMBER, OPLS, and CHARMM). The linear least-squares regression method is used to generate sets of parameters that best satisfy the fitting. Rotational Profiler is free to use, analytical, and force field/package independent. The formalism is herein described, and its usage, in an interactive and automated manner, is made available as a Web server at http://rotprof.lncc.br.

Alpha-hemolysin nanopore allows discrimination of the microcystins variants.

Microcystins (MCs) are a class of cyclic heptapeptides with more than 100 variants produced by cyanobacteria present in surface waters. MCs are potent hepatotoxic agents responsible for fatal poisoning in animals and humans. Several techniques are employed in the detection of MCs, however, there is a shortage of methods capable of discriminating variants of MCs. In this work we demonstrate that the α-hemolysin (αHL) nanopore can detect and discriminate the variants (LR, YR and RR) of MCs in aqueous solution. The discrimination process is based on the analysis of the residence times of each variant of MCs within the unitary nanopore, as well as, on the amplitudes of the blockages in the ionic current flowing through it. Simulations of molecular dynamics and calculation of the electrostatic potential revealed that the variants of MCs present different charge distribution and correlated with the three patterns on the amplitudes of the blockages in the ionic current. Additionally, molecular docking analysis indicates different patterns of interaction of the variants of MCs with two specific regions of the nanopore. We conclude that αHL nanopore can discriminate variants of microcystins by a mechanism based mainly on electrostatic interaction. Finally, we propose the use of nanopore-based technology as a promising method for analyzing microcystins in aqueous solutions.

Aromatic Rings Commonly Used in Medicinal Chemistry: Force Fields Comparison and Interactions With Water Toward the Design of New Chemical Entities.

The identification of lead compounds usually includes a step of chemical diversity generation. Its rationale may be supported by both qualitative (SAR) and quantitative (QSAR) approaches, offering models of the putative ligand-receptor interactions. In both scenarios, our understanding of which interactions functional groups can perform is mostly based on their chemical nature (such as electronegativity, volume, melting point, lipophilicity etc.) instead of their dynamics in aqueous, biological solutions (solvent accessibility, lifetime of hydrogen bonds, solvent structure etc.). As a consequence, it is challenging to predict from 2D structures which functional groups will be able to perform interactions with the target receptor, at which intensity and relative abundance in the biological environment, all of which will contribute to ligand potency and intrinsic activity. With this in mind, the aim of this work is to assess properties of aromatic rings, commonly used for drug design, in aqueous solution through molecular dynamics simulations in order to characterize their chemical features and infer their impact in complexation dynamics. For this, common aromatic and heteroaromatic rings were selected and received new atomic charge set based on the direction and module of the dipole moment from MP2/6-31G* calculations, while other topological terms were taken from GROMOS53A6 force field. Afterwards, liquid physicochemical properties were simulated for a calibration set composed by nearly 40 molecules and compared to their respective experimental data, in order to validate each topology. Based on the reliance of the employed strategy, we expanded the dataset to more than 100 aromatic rings. Properties in aqueous solution such as solvent accessible surface area, H-bonds availability, H-bonds residence time, and water structure around heteroatoms were calculated for each ring, creating a database of potential interactions, shedding light on features of drugs in biological solutions, on the structural basis for bioisosterism and on the enthalpic/entropic costs for ligand-receptor complexation dynamics.

Deriving Structural Information from Experimentally Measured Data on Biomolecules.

During the past half century, the number and accuracy of experimental techniques that can deliver values of observables for biomolecular systems have been steadily increasing. The conversion of a measured value Qexp of an observable quantity Q into structural information is, however, a task beset with theoretical and practical problems: 1) insufficient or inaccurate values of Qexp , 2) inaccuracies in the function Q(r→) used to relate the quantity Q to structure r→ , 3) how to account for the averaging inherent in the measurement of Qexp , 4) how to handle the possible multiple-valuedness of the inverse r→(Q) of the function Q(r→) , to mention a few. These apply to a variety of observable quantities Q and measurement techniques such as X-ray and neutron diffraction, small-angle and wide-angle X-ray scattering, free-electron laser imaging, cryo-electron microscopy, nuclear magnetic resonance, electron paramagnetic resonance, infrared and Raman spectroscopy, circular dichroism, Förster resonance energy transfer, atomic force microscopy and ion-mobility mass spectrometry. The process of deriving structural information from measured data is reviewed with an eye to non-experts and newcomers in the field using examples from the literature of the effect of the various choices and approximations involved in the process. A list of choices to be avoided is provided.

Investigation of the structural preference and flexibility of the loop residues in amyloid fibrils of the HET-s prion.

The structural variability of a 16-residue loop (residues 246-261) which is in part disordered and connects two layers of the ß-solenoid formed by the prion-form of HET-s and its prion domain HET-s(218-289) is investigated using molecular dynamics computer simulation. A system of three HET-s(218-289) molecules in a ß-sheet structure as in the fibril is simulated in aqueous solution. The trajectory structures appear to be consistent with the Cα chemical shift data obtained. In order to delineate the influence of the ß-sheet core of the fibril upon the structural variability of the loop, the latter is also simulated without the ß-sheet core, but with its N- and C-terminal residues restrained at their positions in the fibril. The analysis of the trajectories shows that the structural variability of the loop is restricted by the ß-sheet core, least at its N-terminal end and most in the middle of the trimer.

Assessing protein conformational sampling and structural stability via de novo design and molecular dynamics simulations.

Molecular dynamics and de novo techniques, associated to quality parameter sets, have excelled at determining the structure of small proteins with high accuracy. To achieve a detailed description of protein conformations, these methods must critically assess the thermodynamic features of the molecular ensembles. Here, a comparison of the conformational ensemble generated by molecular dynamics and de novo techniques were carried out for six Top7-based proteins carrying gp41 HIV-1 epitopes. The native Top7, a highly stable computationally designed protein, was used as benchmark. Structural stability, flexibility, and secondary structure content were assessed. The consistency of the latter was compared to experimental circular dichroism spectra for all proteins. While both methods are capable to identify the stable from unstable chimeric proteins, the sampled conformational space and flexibility differ significantly in both methods. Molecular dynamics simulations seem to better describe secondary structure content and identify regions responsible for conformational instability. The de novo method, as implemented in Rosetta-a prime tool for protein design, overestimates secondary structure content. On the other hand, its empirical energy function is capable to predict the threshold for protein stability.

MDWiZ: a platform for the automated translation of molecular dynamics simulations.

A variety of popular molecular dynamics (MD) simulation packages were independently developed in the last decades to reach diverse scientific goals. However, such non-coordinated development of software, force fields, and analysis tools for molecular simulations gave rise to an array of software formats and arbitrary conventions for routine preparation and analysis of simulation input and output data. Different formats and/or parameter definitions are used at each stage of the modeling process despite largely contain redundant information between alternative software tools. Such Babel of languages that cannot be easily and univocally translated one into another poses one of the major technical obstacles to the preparation, translation, and comparison of molecular simulation data that users face on a daily basis. Here, we present the MDWiZ platform, a freely accessed online portal designed to aid the fast and reliable preparation and conversion of file formats that allows researchers to reproduce or generate data from MD simulations using different setups, including force fields and models with different underlying potential forms. The general structure of MDWiZ is presented, the features of version 1.0 are detailed, and an extensive validation based on GROMACS to LAMMPS conversion is presented. We believe that MDWiZ will be largely useful to the molecular dynamics community. Such fast format and force field exchange for a given system allows tailoring the chosen system to a given computer platform and/or taking advantage of a specific capabilities offered by different software engines.

PITOMBA: Parameter Interface for Oligosaccharide Molecules Based on Atoms.

A novel four-bead coarse-grained (CG) model for carbohydrates denoted PITOMBA was devised using a bottom-up approach based on the atomistic GROMOS 53A6GLYC force field and on experimental thermodynamical data. The model was developed to be used in conjunction with the SPC CG water model (J. Chem. Phys. 2011, 134, 084110) and the GROMOS force field functional form. Explicit electrostatic interactions are considered by assigning point charges to each CG bead. Validation of the model is presented to a variety of structural and thermodynamic properties for mono- and oligosaccharides in solution. In addition, the model development philosophy allows for prompt extensions to include hexopyranose chains with diverse glycosidic linkages and branches.

The Effect of Temperature, Cations, and Number of Acyl Chains on the Lamellar to Non-Lamellar Transition in Lipid-A Membranes: A Microscopic View.

Lipopolysaccharides (LPS) are the main constituent of the outer bacterial membrane of Gram-negative bacteria. Lipid-A is the structural region of LPS that interacts with the innate immune system and induces inflammatory responses. It is formed by a phosphorylated ß-d-glucosaminyl-(1→6)-α-N-glucosamine disaccharide backbone containing ester-linked and amide-linked long-chain fatty acids, which may vary in length and number depending on the bacterial strains and the environment. Phenotypical variation (i.e., number of acyl chains), cation type, and temperature influence the phase transition, aggregate structure, and endotoxic activity of Lipid-A. We have applied an extension of the GROMOS force field 45a4 carbohydrate parameter set to investigate the behavior of hexa- and pentaacylated Lipid-A of Pseudomonas aeruginosa at two temperatures (300 and 328 K) and in the presence of mono- and divalent cations (represented by Ca(2+) and Na(+), respectively) through molecular dynamics simulations. The distinct phase of Lipid-A aggregates was characterized by structural properties, deuterium order parameters, the molecular shape of the lipid units (conical versus cylindrical), and molecular packing. Our results show that Na(+) ions induce a transition from the lamellar to nonlamellar phase. In contrast, the bilayer integrity is maintained in the presence of Ca(2+) ions. Through these findings, we present microscopic insights on the influence of different cations on the molecular behavior of Lipid-A associated with the lamellar to nonlamellar transition.

GROMOS 53A6GLYC, an Improved GROMOS Force Field for Hexopyranose-Based Carbohydrates.

An improved parameter set for explicit-solvent simulations of carbohydrates (referred to as GROMOS 53A6GLYC) is presented, allowing proper description of the most stable conformation of all 16 possible aldohexopyranose-based monosaccharides. This set includes refinement of torsional potential parameters associated with the determination of hexopyranose rings conformation by fitting to their corresponding quantum-mechanical profiles. Other parameters, as the rules for third and excluded neighbors, are taken directly from the GROMOS 53A6 force field. Comparisons of the herein presented parameter set to our previous version (GROMOS 45A4), the GLYCAM06 force field, and available NMR data are presented in terms of ring puckering free energies, conformational distribution of the hydroxymethyl group, and glycosidic linkage geometries for 16 selected monosaccharides and eight disaccharides. The proposed parameter modifications have shown a significant improvement for the above-mentioned quantities over the two tested force fields, while retaining full compatibility with the GROMOS 53A6 and 54A7 parameter sets for other classes of biomolecules.