RESUMO

In the context of classical molecular simulations, the accuracy of a force field is highly influenced by the values of the relevant simulation parameters. In this work, a parameter-space mapping (PSM) workflow is proposed to aid in the calibration of force-field parameters, based mainly on the following features: (i) regular-grid discretization of the search space; (ii) partial sampling of the search-space grid; (iii) training of surrogate models to predict the estimates of the target properties for nonsampled parameter sets; (iv) post hoc interpretation of the results in terms of multiobjective optimization concepts; (v) attenuation of statistical errors achieved via empiric extension of the duration of the simulations; (vi) iterative search-space translation according to a user-defined scalar objective function that measures the accuracy of the force field (e.g., the weighted root-mean-square deviation of the target properties relative to the reference data). This combination of features results in a hybrid of a single- and a multiobjective optimization strategy, allowing for the approximate determination of both a local minimum of the chosen objective function and its neighboring Pareto efficient points. The PSM workflow is implemented in the extensible Python program gmak, which is made available in the Git repository at http://github.com/mssm-labmmol/gmak. Using this implementation, the PSM workflow was tested in a proof-of-concept fashion in the recalibration of the Lennard-Jones parameters of the 3-point Optimal Point Charge (OPC3) water model for compatibility with the GROMOS treatment of nonbonded interactions. The recalibrated model reproduces typical pure-liquid properties with an accuracy similar to the original OPC3 model and represents a significant improvement relative to the Simple Point Charge (SPC) model, which is the official recommendation for simulations using GROMOS force fields.

RESUMO

The calibration of torsional interaction terms by fitting relative gas-phase conformational energies against their quantum-mechanical values is a common procedure in force-field development. However, much less attention has been paid to the optimization of third-neighbor nonbonded interaction parameters, despite their strong coupling with the torsions. This article introduces an algorithm termed LLS-SC, aimed at simultaneously parametrizing torsional and third-neighbor interaction terms based on relative conformational energies. It relies on a self-consistent (SC) procedure where each iteration involves a linear least-squares (LLS) regression followed by a geometry optimization of the reference structures. As a proof-of-principle, this method is applied to obtain torsional and third-neighbor interaction parameters for aliphatic chains in the context of the GROMOS 53A6 united-atom force field. The optimized parameter set is compared to the original one, which has been fitted manually against thermodynamic properties for small linear alkanes. The LLS-SC implementation is freely available under http://github.com/mssm-labmmol/profiler.

RESUMO

Amine transaminases (ATAs) are pyridoxal-5'-phosphate (PLP)-dependent enzymes that catalyze the transfer of an amino group from an amino donor to an aldehyde and/or ketone. In the past decade, the enzymatic reductive amination of prochiral ketones catalyzed by ATAs has attracted the attention of researchers, and more traditional chemical routes were replaced by enzymatic ones in industrial manufacturing. In the present work, the influence of the presence of an α,ß-unsaturated system in a methylketone model substrate was investigated, using a set of five wild-type ATAs, the (R)-selective from Aspergillus terreus (Atr-TA) and Mycobacterium vanbaalenii (Mva-TA), the (S)-selective from Chromobacterium violaceum (Cvi-TA), Ruegeria pomeroyi (Rpo-TA), V. fluvialis (Vfl-TA) and an engineered variant of V. fluvialis (ATA-256 from Codexis). The high conversion rate (80 to 99%) and optical purity (78 to 99% ee) of both (R)- and (S)-ATAs for the substrate 1-phenyl-3-butanone, using isopropylamine (IPA) as an amino donor, were observed. However, the double bond in the α,ß-position of 4-phenylbut-3-en-2-one dramatically reduced wild-type ATA reactivity, leading to conversions of <10% (without affecting the enantioselectivity). In contrast, the commercially engineered V. fluvialis variant, ATA-256, still enabled an 87% conversion, yielding a corresponding amine with >99% ee. Computational docking simulations showed the differences in orientation and intermolecular interactions in the active sites, providing insights to rationalize the observed experimental results.

Assuntos

RESUMO

The construction of a molecular topology file is a prerequisite for any classical molecular dynamics simulation. However, the generation of such a file may be very challenging at times, especially for large supramolecules. While many tools are available to provide topologies for large proteins and other biomolecules, the scientific community researching nonbiological systems is not equally well equipped. Here, we present a practical tool to generate topologies for arbitrary supramolecules: The pyPolyBuilder. In addition to linear polymer chains, it also provides the possibility to generate topologies of arbitrary, large, branched molecules, such as, e.g., dendrimers. Furthermore, it also generates reasonable starting structures for simulations of these molecules. pyPolyBuilder is a standalone command-line tool implemented in python. Therefore, it may be easily incorporated in persisting simulation pipelines on any operating systems and with different simulation engines. pyPolyBuilder is freely available on github: https://github.com/mssm-labmmol/pypolybuilder.

Assuntos

RESUMO

The complexation of quercetin molecules with poly(amidoamine) (PAMAM) dendrimers of generation 0-3 was studied by molecular dynamics simulations. Three main points were addressed: (i) the effect of starting from different initial structures; (ii) the performance of the 2016H66 force field (recently validated in the context of dendrimer simulations) in predicting the experimental drug(quercetin)-loading capacity of PAMAM dendrimers; and (iii) the stability of quercetin-PAMAM complexes and their interactions. Initial structures generated by different restraint protocols led to faster convergence compared to initial structures generated by randomly placing the drug molecules in the simulation box. The simulations yielded meta-stable complexes where the loading numbers have converged to average values and were compared to experimentally obtained values. Once the first meta-stable state was reached, the drug-dendrimer complexes did not deviate significantly throughout the simulation. They were characterized in terms of structural properties, such as the radius of gyration and radial distribution functions. The results suggest that quercetin molecules interact mostly with the internal dendrimer monomers rather than to their surface.

Assuntos

RESUMO

Alkanes are a fundamental part in empirical force fields (FF) not only due to their technological relevance, but also due to the prevalence of alkane moieties in organic molecules, e.g., compounds containing a saturated carbon chain. Therefore, a good description of alkane interactions is crucial for determining the quality of a FF. In this study, the performance of 12 empirical force fields (FF) was evaluated in the context of reproducing liquid properties of alkanes. More specifically, n-octane was chosen as a reference compound since it is a liquid in a broad temperature range and it has numerous experimental data for thermodynamic, transport, and structural properties, as well as for their temperature dependencies. A normalized root-mean-square deviation (NRMSD) analysis was used to rank the force fields in their ability to reproduce the experimental data. Five out of the six best force fields considered were united-atom models. The GROMOS force field showed the smallest deviation in terms of NRMSD, followed by TRAPPE-EH, NERD, CHARMM-UA, TRAPPE-UA, and OPLS-UA. This overall better performance of the united-atom force fields indicates that complexity does not always bring quality.

RESUMO

A novel methodology for the 1,1-dichlorocyclopropanation of dicarbonyl conjugated olefins was described. The developed protocol is simple and uses readily accessible starting materials, allowing the isolation of the desired adducts in moderate to excellent yields (up to 99 %). Furthermore, the reaction tolerated scale up to the gram scale; thus highlighting the synthetic potential of this transformation. Control experiments and DFT studies revealed that the reaction proceeded through a Michael-initiated ring-closure process, in which reaction temperature played a crucial role. Finally, these gem-dichlorocyclopropanes were also employed in the preparation of a trisubstituted naphthyl derivative and a diastereoselective reduction was also demonstrated.

RESUMO

A systematic evaluation of the accuracy of the GROMOS-compatible 2016H66 force field in the simulation of dendrimers is performed. More specifically, the poly(amido amine) (PAMAM) and the poly(propyleneimine) (PPI) are considered because of the availability of experimental data and simulation results in the literature. A total of 36 molecular systems are simulated and the radius of gyration, asphericity, density profiles, and the self-diffusion coefficient are monitored in terms of the generation number and pH (low, medium, and high) condition. Overall, the results support the recommendation of the 2016H66 force field for the simulation of dendrimer systems. The natural building-block based strategy adopted in the definition of 2016H66, together with a careful parametrization of the chemical functional groups to reproduce thermodynamic properties in environments of different polarity, and also the ability to accurately reproduce the expected structural and dynamic features of dendrimers, as shown in the present work, make this force field an attractive option for the simulation of such systems.

Assuntos

RESUMO

The effect of different treatments of the nonbonded interactions in simulations employing the recently introduced GROMOS-compatible 2016H66 force field is evaluated based on calculations carried out with the GROMACS software. This is done considering four thermodynamic and transport properties (pure liquid density, vaporization enthalpy, surface-tension coefficient, and self-diffusion constant) of 58 organic liquids representative of the chemical groups alcohol, ether, aldehyde, ketone, carboxylic acid, ester, amine, amide, thiol, sulfide, disulfide, and aromatic compounds, also including water (SPC model). A dipalmitoylphosphatidylcholine bilayer system is considered as well. The simulated properties are found to be very sensitive to the treatment of the long-range dispersion interactions, notably for the least polar systems. In general, the treatment of the long-range electrostatic or Lennard-Jones interactions using homogeneous correction terms or lattice-sum approaches yield similar results, with punctual discrepancies. The combination of a lattice-sum approach for the electrostatic interactions with a straight-cutoff truncation of the Lennard-Jones interactions at a distance of at least 1.2 nm is found to represent a good compromise setup within GROMACS for achieving compatibility with the reference results obtained using GROMOS as well as a comparable level of agreement with the experimental data. This study also reveals two potential issues with the GROMACS software, related to an incorrect calculation of the pressure when using LINCS in version 4.0.7 and an inadequate implementation of the twin-range scheme in version 5.1.2.

RESUMO

Diseases caused by flaviviruses, such as dengue and zika, are globally recognized as major threats. During infection, a critical point in their replicative cycle is the maturation step, which occurs throughout the cellular exocytic pathway. This step is a pH-dependent process that involves the modification of the viral envelope by converting prM (pre-membrane) into M (membrane) proteins with the release of a "pr peptide". After this reaction, the pr peptides remain bound to the viral envelope while the virions cross the acidic trans-Golgi network, and are released only at neutral pH after secretion of the virus particles. Despite this current knowledge, the molecular basis of the flavivirus maturation step is largely unknown. Here, based on the crystal structure of the dengue pr-E complex ("pr peptide" bound to virus envelope protein) and using molecular dynamics simulations, we found that the pH shift from acidic to neutral yields considerable structural changes in the system. Dynamic cross correlation maps and root mean square deviation analyses revealed that the pr-E junction is clearly unstable under neutral pH. Secondary structure analysis also revealed that the fusion loop region, present in the E protein, is sensitive to pH and tends to unstructure at a neutral environment. Moreover, we found that five residues present in the E protein, Gly102, His244, Thr70, Thr68 and Asn67 are critical to confer stability to the pr-E complex while inside the Golgi apparatus. This work brings details about the dynamical behavior of the pr-E system, helps to better understand the flavivirus biology and may also be of use in the development of novel antiviral strategies.

Assuntos

RESUMO

Chemokine ligand 2 (CCL2), also known as monocyte chemoattractant protein 1 (MCP-1), is a chemokine that recruits immune cells to inflammatory sites by interacting with G protein-coupled receptor CCR2. The CCL2/CCR2 axis is also involved in pathological processes such as tumor growth and metastasis and hence is currently considered as an important drug target. CCL2 exists in a dynamic monomer-dimer equilibrium that is modulated by CCR2 binding. We used solution nuclear magnetic resonance (NMR) spectroscopy and molecular dynamics simulations to study the interactions between CCL2 and a sulfopeptide corresponding to the N-terminal sequence of CCR2 (CCR218-31). Peptide binding induced the dissociation of CCL2 into monomers, forming stable CCL2/CCR218-31 complexes. NMR relaxation measurements indicated that residues around the CCR218-31 binding site, which are located at the dimer interface, undergo a complex regime of motions. NMR data were used to construct a three-dimensional structural model of the CCL2/CCR218-31 complex, revealing that CCR218-31 occupies a binding site juxtaposed to the dimer interface, partially replacing monomer-monomer contacts, explaining why CCR218-31 binding weakens the dimer interface and induces dissociation. We found that the main interactions governing receptor binding are highly stable salt bridges with conserved chemokine residues as well as hydrophobic interactions. These data provide new insights into the structure-function relationship of the CCL2-CCR2 interaction and may be helpful for the design of novel antichemotactic agents.

Assuntos

RESUMO

The flavivirus non-structural protein 1 (NS1) is a conserved glycoprotein with as yet undefined biological function. This protein dimerizes when inside infected cells or associated to cell membranes but also forms lipid-associated hexamers when secreted to the extracellular space. A single amino acid substitution (P250L) is capable of preventing the dimerization of NS1 resulting in lower virulence and slower virus replication. In this work, based on molecular dynamics simulations of the dengue-2 virus NS1 [Formula: see text]-ladder monomer as a core model, we found that this mutation can induce several conformational changes that importantly affect critical monomer-monomer interactions. Based on additional simulations, we suggest a mechanism by which a highly orchestrated sequence of events propagate the local perturbations around the mutation site towards the dimer interface. The elucidation of such a mechanism could potentially support new strategies for rational production of live-attenuated vaccines and highlights a step forward in the development of novel anti-flavivirus measures.

Assuntos

RESUMO

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.

Assuntos

RESUMO

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.

Assuntos

RESUMO

The vascular endothelial growth factor (VEGF) seems to be the most important regulator of physiological and pathological angiogenesis, being, for this reason, a favorite target for therapies against angiogenesis-related diseases. VEGF is a homodimer in which the monomers are formed by beta-strands interconnected on the poles by three loops. A recent work showed that an intimate relationship between loops-1 and -3 is required for high affinity binding to the receptors (Kiba et al., J Biol Chem 2003;278:13453-13461). In this work, we report the results of a 10-ns molecular dynamics simulation of VEGF. We analyzed the dynamical behavior of the protein (using a dynamical cross-correlation map) and found that it is governed by a high degree of correlation between the motions of the loops. We also performed a principal component analysis and found an overall motion in which the opposite poles are projected against each other, just like the movement of the wings of a butterfly. From the biological point of view, it is likely that this motion would facilitate receptor binding since VEGF must enter a restricted cavity formed by the two subunits of the receptor.

Assuntos