Molecular Design and Spodumene Flotation-A Review.

Spodumene flotation stands as the most commonly used method to concentrate lithium minerals. However, it faces significant challenges related to low collector recoveries and similarity in the surface characteristics of the minerals, which make the effective separation of this valuable mineral difficult. For this reason, numerous researchers have conducted studies to address and confront this problem. In this work, an exhaustive bibliographic search was carried out using keywords and search queries, and the results were structured in three sections according to temporal, methodological, and thematic criteria. The first section covers the period from 1950 to 2004, focusing on experimental tests. The second section covers from 2004 to the present and focuses on flotation tests and measurement analysis. Simultaneously, the third section spans from 2011 to the present and is based on molecular dynamics simulations. Topics covered include spodumene surface properties, the influence of metal ions, pre-treatment techniques, and the use of collectors. Ultimately, molecular dynamics simulations are positioned as a tool that accurately represents experimental phenomena. In this context, specialized software such as Materials Studio or Gromacs prove to be reliable instruments that allow a detailed study of mineral surfaces and other elements to be carried out, which justifies their consideration for future research in this scientific field.

Visual dynamics: a WEB application for molecular dynamics simulation using GROMACS.

BACKGROUND: The molecular dynamics is an approach to obtain kinetic and thermodynamic characteristics of biomolecular structures. The molecular dynamics simulation softwares are very useful, however, most of them are used in command line form and continue with the same common implementation difficulties that plague researchers who are not computer specialists. RESULTS: Here, we have developed the VisualDynamics-a WEB tool developed to automate biological simulations performed in Gromacs using a graphical interface to make molecular dynamics simulation user-friendly task. In this new application the researcher can submit a simulation of the protein in the free form or complexed with a ligand. Can also download the graphics analysis and log files at the end of the simulation. CONCLUSIONS: VisualDynamics is a tool that will accelerate implementations and learning in the area of molecular dynamics simulation. Freely available at https://visualdynamics.fiocruz.br/login , is supported by all major web browsers. VisualDynamics was developed with Flask, which is a Python-based free and open-source framework for web development. The code is freely available for download at GitHub https://github.com/LABIOQUIM/visualdynamics .

Out of Sight, Out of Mind: The Effect of the Equilibration Protocol on the Structural Ensembles of Charged Glycolipid Bilayers.

Molecular dynamics (MD) simulations represent an essential tool in the toolbox of modern chemistry, enabling the prediction of experimental observables for a variety of chemical systems and processes and majorly impacting the study of biological membranes. However, the chemical diversity of complex lipids beyond phospholipids brings new challenges to well-established protocols used in MD simulations of soft matter and requires continuous assessment to ensure simulation reproducibility and minimize unphysical behavior. Lipopolysaccharides (LPS) are highly charged glycolipids whose aggregation in a lamellar arrangement requires the binding of numerous cations to oppositely charged groups deep inside the membrane. The delicate balance between the fully hydrated carbohydrate region and the smaller hydrophobic core makes LPS membranes very sensitive to the choice of equilibration protocol. In this work, we show that the protocol successfully used to equilibrate phospholipid bilayers when applied to complex lipopolysaccharide membranes occasionally leads to a small expansion of the simulation box very early in the equilibration phase. Although the use of a barostat algorithm controls the system dimension and particle distances according to the target pressure, fluctuation in the fleeting pressure occasionally enables a few water molecules to trickle into the hydrophobic region of the membrane, with spurious solvent buildup. We show that this effect stems from the initial steps of NPT equilibration, where initial pressure can be fairly high. This can be solved with the use of a stepwise-thermalization NVT/NPT protocol, as demonstrated for atomistic MD simulations of LPS/DPPE and lipid-A membranes in the presence of different salts using an extension of the GROMOS forcefield within the GROMACS software. This equilibration protocol should be standard procedure for the generation of consistent structural ensembles of charged glycolipids starting from atomic coordinates not previously pre-equilibrated. Although different ways to deal with this issue can be envisioned, we investigated one alternative that could be readily available in major MD engines with general users in mind.

Intramolecular energies of the cytotoxic protein CagA of Helicobacter pylori as a possible descriptor of strains' pathogenicity level.

The Helicobacter pylori cytotoxin-associated gene A (CagA) is known for causing gastroduodenal diseases, such as atrophic gastritis and peptic ulcerations. Furthermore Helicobacter pylori CagA positive strains has been reported as one of the main risk factors for gastric cancer (Parsonnet et al., 1997). Structural variations in the CagA structure can alter its affinity with the host proteins, inducing differences in the pathogenicity of H. pylori. CagA N-terminal region is characterized for be conserved among all H. pylori strains since the C-terminal region is characterized by an intrinsically disorder behavior. We generated complete structural models of CagA using different conformations of the C-terminal region for two H. pylori strains. These models contain the same EPIYA (ABC1C2) motifs but different level of pathogenicity: gastric cancer and duodenal ulcer. Using these structural models we evaluated the pathogenicity level of the H. pylori strain, based on the affinity of the interaction with SHP-2 and Grb2 receptors and on the number of interactions with the EPIYA motif. We found that the main differences in the interaction was due to the contributions of certain types of energies from each strain and not from the total energy of the molecule. Specifically, the electrostatic energy, helix dipole energy, Wander Waals clashes, torsional clash, backbone clash and cis bond energy allowed a separation between severe and mild pathology for the interaction of only CagA with SHP2.

LiGRO: a graphical user interface for protein-ligand molecular dynamics.

To speed up the drug-discovery process, molecular dynamics (MD) calculations performed in GROMACS can be coupled to docking simulations for the post-screening analyses of large compound libraries. This requires generating the topology of the ligands in different software, some basic knowledge of Linux command lines, and a certain familiarity in handling the output files. LiGRO-the python-based graphical interface introduced here-was designed to overcome these protein-ligand parameterization challenges by allowing the graphical (non command line-based) control of GROMACS (MD and analysis), ACPYPE (ligand topology builder) and PLIP (protein-binder interactions monitor)-programs that can be used together to fully perform and analyze the outputs of complex MD simulations (including energy minimization and NVT/NPT equilibration). By allowing the calculation of linear interaction energies in a simple and quick fashion, LiGRO can be used in the drug-discovery pipeline to select compounds with a better protein-binding interaction profile. The design of LiGRO allows researchers to freely download and modify the software, with the source code being available under the terms of a GPLv3 license from http://www.ufrgs.br/lasomfarmacia/ligro/ .

A molecular dynamics study of the BACE1 conformational change from Apo to closed form induced by hydroxyethylamine derived compounds.

BACE1 is an aspartyl protease which is a therapeutic target for Alzheimer's disease (AD) because of its participation in the rate-limiting step in the production of Aß-peptide, the accumulation of which produces senile plaques and, in turn, the neurodegenerative effects associated with AD. The active site of this protease is composed in part by two aspartic residues (Asp93 and Asp289). Additionally, the catalytic site has been found to be covered by an antiparallel hairpin loop called the flap. The dynamics of this flap are fundamental to the catalytic function of the enzyme. When BACE1 is inactive (Apo), the flap adopts an open conformation, allowing a substrate or inhibitor to access the active site. Subsequent interaction with the ligand induces flap closure and the stabilization of the macromolecular complex. Further, the protonation state of the aspartic dyad is affected by the chemical nature of the species entering the active site, so that appropriate selection of protonation states for the ligand and the catalytic residues will permit the elucidation of the inhibitory pathway for BACE1. In the present study, comparative analysis of different combinations of protonation states for the BACE1-hydroxyethylamine (HEA) system is reported. HEAs are potent inhibitors of BACE1 with favorable pharmacological and kinetic properties, as well as oral bioavailability. The results of Molecular Dynamics (MD) simulations and population density calculations using 8 different parameters demonstrate that the LnAsp289 configuration (HEA with a neutral amine and the Asp289 residue protonated) is the only one which permits the expected conformational change in BACE1, from apo to closed form, after flap closure. Additionally, differences in their capacities to establish and maintain interactions with residues such as Asp93, Gly95, Thr133, Asp289, Gly291, and Asn294 during this step allow differentiation among the inhibitory activities of the HEAs. The results and methodology here reported will serve to elucidate the inhibitory pathway of other families of compounds that act as BACE1 inhibitors, as well as the design of better leader compounds for the treatment of AD.

Determination of the protonation state for the catalytic dyad in ß-secretase when bound to hydroxyethylamine transition state analogue inhibitors: A molecular dynamics simulation study.

BACE1 is an aspartyl protease of pharmacological interest for its direct participation in Alzheimer's disease (AD) through ß-amyloid peptide production. Two aspartic acid residues are present in the BACE1 catalytic region which can adopt multiple protonation states depending on the chemical nature of its inhibitors, i.e., monoprotonated, diprotonated and di-deprotonated states. In the present study a series of protein-ligand molecular dynamics (MD) simulations was carried out to identify the most feasible protonation state adopted by the catalytic dyad in the presence of hydroxyethylamine transition state analogue inhibitors. The MD trajectories revealed that the di-deprotonated state is most prefered in the presence of hydroxyethilamine (HEA) family inhibitors. This appears as a result after evaluating, for all 9 protonation state configurations during the simulation time, the deviations of a set of distances and dihedral angles measured on the ligand, protein and protein-ligand complex with reference to an X-ray experimental BACE1/HEA crystallographic structure. These results will help to clarify the phenomena related to the HEAs inhibitory pathway, and improve HEAs databases' virtual screening and ligand design processes targeting ß-secretase protein.