Influence of non-protein amino-acid mimosine in peptide conformational propensities from novel Amber force field parameters

Mimosine is a non-protein amino acid derived from plants known for its ability to bind to divalent or trivalent metal cations such as Zn$^{2+}$, Ni$^{2+}$, Fe$^{2+}$ or Al$^{3+}$. This results in interesting antimicrobial and anti-cancer properties, which make mimosine a promising candidate for therapeutic applications. One possibility is to incorporate mimosine into synthetic short peptide drugs. However, our understanding of how this amino acid affects peptide structure is still limited, reducing our ability to design effective therapeutic compounds. In this work, we used computer simulations to understand this question. We first build parameters for the mimosine residue to be used in combination with two classical force fields of the Amber family. Then, we used atomistic molecular dynamics simulations with the resulting parameter sets to evaluate the influence of mimosine in the structural propensities for this amino acid. We compared the results of these simulations with identical peptides where mimosine is replaced by either phenylalanine or tyrosine. We found that the strong dipole in mimosine induces a preference for conformations where the amino acid rings are stacked over more traditional conformations. We validated our results using quantum mechanical calculations, which provide a robust foundation to the outcome of our classical simulations.

Computing the Structural Dynamics of RVFV L Protein Domain in Aqueous Glycerol Solutions

Many biological and biotechnological processes are controlled by protein–protein and protein–solvent interactions. In order to understand, predict, and optimize such processes, it is important to understand how solvents affect protein structure during protein–solvent interactions. In this study, all-atom molecular dynamics are used to investigate the structural dynamics and energetic properties of a C-terminal domain of the Rift Valley Fever Virus L protein solvated in glycerol and aqueous glycerol solutions in different concentrations by molecular weight. The Generalized Amber Force Field is modified by including restrained electrostatic potential atomic charges for the glycerol molecules. The peptide is considered in detail by monitoring properties like the root-mean-squared deviation, root-mean-squared fluctuation, radius of gyration, hydrodynamic radius, end-to-end distance, solvent-accessible surface area, intra-potential energy, and solvent–peptide interaction energies for hundreds of nanoseconds. Secondary structure analysis is also performed to examine the extent of conformational drift for the individual helices and sheets. We predict that the peptide helices and sheets are maintained only when the modeling strategy considers the solvent with lower glycerol concentration. We also find that the solvent-peptide becomes more cohesive with decreasing glycerol concentrations. The density and radial distribution function of glycerol solvent calculated when modeled with the modified atomic charges show a very good agreement with experimental results and other simulations at 298.15K.

Modeling of Self-Assembled Peptide Nanotubes and Determination of Their Chirality Sign Based on Dipole Moment Calculations

The chirality quantification is of great importance in structural biology, where the differences in proteins twisting can provide essentially different physiological effects. However, this aspect of the chirality is still poorly studied for helix-like supramolecular structures. In this work, a method for chirality quantification based on the calculation of scalar triple products of dipole moments is suggested. As a model structure, self-assembled nanotubes of diphenylalanine (FF) made of L- and D-enantiomers were considered. The dipole moments of FF molecules were calculated using semi-empirical quantum-chemical method PM3 and the Amber force field method. The obtained results do not depend on the used simulation and calculation method, and show that the D-FF nanotubes are twisted tighter than L-FF. Moreover, the type of chirality of the helix-like nanotube is opposite to that of the initial individual molecule that is in line with the chirality alternation rule general for different levels of hierarchical organization of molecular systems. The proposed method can be applied to study other helix-like supramolecular structures.

Density, Enthalpy of Vaporization and Local Structure of Neat N-Alkane Liquids

The properties of alkanes are consequential for understanding many chemical processes in nature and industry. We use molecular dynamics simulations with the Amber force field GAFF2 to examine the structure of pure liquids at each respective normal boiling point, spanning the 15 n-alkanes from methane to pentadecane. The densities predicted from the simulations are found to agree well with reported experimental values, with an average deviation of 1.9%. The enthalpies of vaporization have an average absolute deviation from experiment of 10.4%. Radial distribution functions show that short alkanes have distinct local structures that are found to converge with each other with increasing chain length. This provides a unique perspective on trends in the n-alkane series and will be useful for interpreting similarities and differences in the n-alkane series as well as the breakdown of ideal solution behavior in mixtures of these molecules.

Comparison of different approaches to derive classical bonded force-field parameters for a transition metal cofactor: a case study for non-heme iron site of ectoine synthase

Computational investigations into the structure and function of metalloenzymes with transition metal cofactors require proper preparation of the model, which requires obtaining reliable force field parameters for the cofactor. Here, we present a test case where several methods were used to derive amber force field parameters for a bonded model of the Fe(II) cofactor of ectoine synthase. Moreover, the spin of the ground state of the cofactor was probed by DFT and post-HF methods, which consistently indicated the quintet state is lowest in energy and well separated from triplet and singlet. The performance of the obtained force field parameter sets, derived for the quintet spin state, was scrutinized and compared taking into account metrics focused on geometric features of the models as well as their energetics. The main conclusion of this study is that Hessian-based methods yield parameters which represent the geometry around the metal ion, but poorly reproduce energy variance with geometrical changes. On the other hand, the energy-based method yields parameters accurately reproducing energy-structure relationships, but with bad performance in geometry optimization. Preliminary tests show that admixing geometrical criteria to energy-based methods may allow to derive parameters with acceptable performance for both energy and geometry.

Polyacrylic Acid to Improve Flotation Tailings Management: Understanding the Chemical Interactions through Molecular Dynamics

Molecular dynamic simulations of polyacrylic acid polyelectrolyte (PAA) analyzed its interaction with the main minerals that make up characteristic tailings of the mining industry, in this case, quartz, kaolinite, and montmorillonite. The simulations were carried out with the package Gromacs 2020.3. The interaction potentials used were General AMBER Force Field (GAFF) for PAA and CLAYFF-MOH for mineral surfaces. The SPC/E model described water molecules and Lennard-Jones 12-6 parameters adjusted for SPC/E model were used for Na+ and Cl− ions. The studied systems were carried out at pH 7, obtaining stable adsorption between the PAA and the studied surfaces. Interestingly, the strongest adsorptions were for montmorillonite at both low and high salt concentrations. The effect of salinity differs according to the system, finding that it impairs the absorption of the polymer on montmorillonite surfaces. However, a saline medium favors the interaction with quartz and kaolinite. This is explained because montmorillonite has a lower surface charge density and a greater capacity to adsorb ions. This facilitated the adsorption of PAA. It was possible to identify that the main interaction by which the polymer is adsorbed is through the hydroxyl of the mineral surface and the COO−Na+ complexes. Molecular dynamics allows us to advance in the understanding of interactions that define the behavior of this promising reagent as an alternative for sustainable treatment of complex tailings in highly saline environments.

Optimized Halogen Atomic Radii for PBSA Calculations Using Off-Center Point-Charges

<div>In force field methods, the usage of off-center point-charges, also called extra-points (EPs), is a common strategy to tackle the anisotropy of the electrostatic potential of covalently-bonded halogens (X), thus allowing the description of halogen bonds (XBs) at the molecular mechanics / molecular dynamics (MM/MD) level. Diverse EP implementations exist in the literature differing on the charge sets and/or the X–EP distances. Poisson–Boltzmann and surface area (PBSA) calculations can be used to obtain solvation free energies (∆G solv ) of small molecules, often to compute binding free energies (∆G bind ) at the MM PBSA level. This method depends, among other parameters, on the empirical assignment of atomic radii (PB radii). Given the multiplicity of off-center point-charges models and the lack of specific PB radii for halogens compatible with such implementations, in this work we assessed the performance of PBSA calculations for the estimation of ∆G solv values in water (∆G hyd ), also conducting an optimization of the halogen PB radii (Cl, Br, and I) for each EP model. We not only expand the usage of EP models in the scope of the General AMBER Force Field (GAFF) but also provide the first optimized halogen PB radii in the context of the CHARMM General Force Field (CGenFF), thus contributing to improving the description of halogenated compounds in PBSA calculations.</div>

Force Field Parameters for Fe2+4S2−4 Clusters of Dihydropyrimidine Dehydrogenase, the 5-Fluorouracil Cancer Drug Deactivation Protein: A Step towards In Silico Pharmacogenomics Studies

The dimeric dihydropyrimidine dehydrogenase (DPD), metalloenzyme, an adjunct anti-cancer drug target, contains highly specialized 4 × Fe2+4S2−4 clusters per chain. These clusters facilitate the catalysis of the rate-limiting step in the pyrimidine degradation pathway through a harmonized electron transfer cascade that triggers a redox catabolic reaction. In the process, the bulk of the administered 5-fluorouracil (5-FU) cancer drug is inactivated, while a small proportion is activated to nucleic acid antimetabolites. The occurrence of missense mutations in DPD protein within the general population, including those of African descent, has adverse toxicity effects due to altered 5-FU metabolism. Thus, deciphering mutation effects on protein structure and function is vital, especially for precision medicine purposes. We previously proposed combining molecular dynamics (MD) and dynamic residue network (DRN) analysis to decipher the molecular mechanisms of missense mutations in other proteins. However, the presence of Fe2+4S2−4 clusters in DPD poses a challenge for such in silico studies. The existing AMBER force field parameters cannot accurately describe the Fe2+ center coordination exhibited by this enzyme. Therefore, this study aimed to derive AMBER force field parameters for DPD enzyme Fe2+ centers, using the original Seminario method and the collation features Visual Force Field Derivation Toolkit as a supportive approach. All-atom MD simulations were performed to validate the results. Both approaches generated similar force field parameters, which accurately described the human DPD protein Fe2+4S2−4 cluster architecture. This information is crucial and opens new avenues for in silico cancer pharmacogenomics and drug discovery related research on 5-FU drug efficacy and toxicity issues.

Data for Molecular Dynamics Simulations of Escherichia Coli Cytochrome Bd Oxidase with the Amber Force Field

<div> <div> <div> <div> <p>Cytochrome <i>bd</i>-type quinol oxidase is an important metalloenzyme that allows many bacteria to survive in low oxygen conditions. Since bd oxidase is found in many prokaryotes but not in eukaryotes, it has emerged as a promising bacterial drug target. Examples of organisms containing bd oxidases include the <i>Mycobacterium tuberculosis (Mtb)</i> bacterium that causes tuberculosis (TB) in humans, the <i>Vibrio cholerae</i> bacterium that causes cholera, the <i>Pseudomonas aeruginosa</i> bacterium that contributes to antibiotic resistance and sepsis, and the <i>Campylobacter jejuni</i> bacterium that causes food poisoning. <i>Escherichia coli (E. coli)</i> is another organism exhibiting the cytochrome <i>bd</i> oxidase. Since it has the highest sequence identity to <i>Mtb</i> (36 %) and we are ultimately interested in finding drug targets for TB, we have built parameters for the <i>E. coli bd </i>oxidase (Protein Data Bank ID number: 6RKO) that are compatible with the all-atom Amber ff14SB force field for molecular dynamics (MD) simulations. Specifically, we built parameters for the three heme cofactors present in all species of bacterial cytochrome <i>bd</i>-type oxidases (heme b₅₅₈, heme b₅₉₅, and heme d) along with their axial ligands. This data report includes the parameter files that can be used with Amber's LEaP program to generate input files for MD simulations using the Amber software package. We also provide the PDB data files of the initial model both by itself and solvated with TIP3P water molecules and counterions. </p> </div> </div> </div> </div>