The ACPYPE web server for small-molecule MD topology generation.

MOTIVATION: The generation of parameter files for molecular dynamics (MD) simulations of small molecules that are suitable for force fields commonly applied to proteins and nucleic acids is often challenging. The ACPYPE software and website aid the generation of such parameter files. RESULTS: ACPYPE uses OpenBabel and ANTECHAMBER to generate MD input files in Gromacs, AMBER, CHARMM, and CNS formats. It can now take a SMILES string as input, in addition to the original PDB or mol2 coordinate files, with GAFF2 support and GLYCAM force field conversion added. It can be installed locally via Anaconda, PyPI, and Docker distributions, while the web server at https://bio2byte.be/acpype/ was updated with an API, and provides visualization of results for uploaded molecules as well as a pre-generated set of 3738 drug molecules. AVAILABILITY AND IMPLEMENTATION: The web application is freely available at https://www.bio2byte.be/acpype/ and the open-source code can be found at https://github.com/alanwilter/acpype.

First principles calculations of thermodynamics and kinetic parameters and molecular dynamics simulations of acetylcholinesterase reactivators: can mouse data provide new insights into humans?

We have applied a theoretical methodology, previously developed to evaluate the association and kinetic reactivation constants of oximes, comparing theoretical data obtained for human acetylcholinesterase (HsAChE) with in vitro results from Mus musculus AChE (MmAChE) previously reported in the literature. Our results, further checked by additional molecular dynamics simulations steps, showed a good correlation between the theoretical and experimental data, supporting the methodology as appropriate for prediction of thermodynamic and kinetic parameters and corroborated MmAChE as a suitable model for studies with HsAChE.

Molecular dynamics simulations applied to the study of subtypes of HIV-1 protease common to Brazil, Africa, and Asia.

Africa accounts for the majority of HIV-1 infections worldwide caused mainly by the A and C viral subtypes rather than B subtype, which prevails in the United States and Western Europe. In Brazil, B subtype is the major subtype, but F, C, and A also circulate. These non-B subtypes present polymorphisms, and some of them occur at sites that have been associated with drug resistance, including the HIV-1 protease (PR), one important drug target. Here, we report a Molecular Dynamics study of the B and non-B PR complexed with the inhibitor ritonavir to delineate the behavior of each subtype. We compare root mean squared deviation, binding free energy by linear interaction energy approach, hydrogen bonds, and intermolecular contact surface area between inhibitor and PR. From our results, we can provide a basis to understand the molecular mechanism of drug resistance in non-B subtypes. In this sense, we found a decrease of approx 4 kcal/mol in deltaG of binding between B and non-B subtypes. This corresponds to the loss of one hydrogen bond, which is in agreement with our H-bond analysis. Previous experimental affinity studies reported analogous results with inhibition constant values for non-B PR.