RESUMO
Molecularly imprinted polymers (MIPs) have emerged as bespoke materials with versatile molecular applications. In this study, we propose a proof of concept for a methodology employing molecular dynamics (MD) simulations to guide the selection of functional monomers for curcuminoid binding in MIPs. Curcumin, demethoxycurcumin, and bisdemethoxycurcumin are phenolic compounds widely employed as spices, pigments, additives, and therapeutic agents, representing the three main curcuminoids of interest. Through MD simulations, we investigated prepolymerization mixtures composed of various functional monomers, including acrylamide (ACA), acrylic acid (AA), methacrylic acid (MAA), and N-vinylpyrrolidone (NVP), with ethylene glycol dimethacrylate (EGDMA) as the cross-linker and acetonitrile as the solvent. Curcumin was selected as the template molecule due to its structural similarity to the other curcuminoids. Notably, the prepolymerization mixture containing NVP as the functional monomer demonstrated superior molecular recognition capabilities toward curcumin. This observation was supported by higher functional monomer molecules surrounding the template, a lower total nonbonded energy between the template and monomer, and a greater number of hydrogen bonds in the aggregate. These findings suggest a stronger affinity between the functional monomer NVP and the template. We synthesized, characterized, and conducted binding tests on the MIPs to validate the MD simulation results. The experimental binding tests confirmed that the MIP-NVP exhibited higher binding capacity. Consequently, based on MD simulations, our computational methodology effectively guided the selection of the functional monomer, leading to MIPs with binding capacity for curcuminoids. The outcomes of this study provide a valuable reference for the rational design of MIPs through MD simulations, facilitating the selection of components for MIPs. This computational approach holds the potential for extension to other templates, establishing a robust methodology for the rational design of MIPs.
Assuntos
Curcumina , Simulação de Dinâmica Molecular , Polímeros Molecularmente Impressos , Curcumina/química , Curcumina/análogos & derivados , Curcumina/metabolismo , Polímeros Molecularmente Impressos/química , Desenho de Fármacos , Impressão Molecular , Metacrilatos/química , Diarileptanoides/química , Conformação MolecularRESUMO
This work reports on new software for automatic conformer energy benchmarking calculations for flexible molecules. The software workflow consists of four parts: conformational search, preoptimization, optimization, and frequency calculations at a higher level and last calculations using several theoretical levels. The software was written to be user-friendly and versatile to be used by nonexperts in computational chemistry. Any theoretical levels available in either Gaussian 16 or ORCA 5 may be applied in the benchmarking study. The workflow will automatically run conformational search calculations and deal with conformers that converge to the same minimum and those that show a negative frequency. At the end of the workflow, the user will have the mean absolute deviations and the most accurate method/DFT functional and basis set in comparison to the benchmark to be applied for the molecular system of interest. Case examples are given at the end of the paper that may help users to get insight into the software's main features.
Assuntos
Benchmarking , Software , Conformação Molecular , Automação , Elétrons , Modelos MolecularesRESUMO
Staurosporine and its analogs (STA-analogs) are indolocarbazoles (ICZs) compounds able to inhibit kinase proteins in a non-specific way, while present antimicrobial and cytostatic properties. The knowledge of molecular features associated to the complexation, including the ligand shape in solution and thermodynamics of complexation, is substantial to the development of new bioactive ICZs with improved therapeutic properties. In this context, the empirical approach of GROMOS force field is able to accurately reproduce condensed phase physicochemical properties of molecular systems after parameterization. Hence, through parameterization under GROMOS force field and molecular simulations, we assessed STA-analogs dynamics in aqueous solution, as well as its interaction with water to probe conformational and structural features involved in complexation to therapeutic targets. The coexistence of multiple conformers observed in simulations, and confirmed by metadynamics calculations, expanding the conformational space knowledge of these ligands with potential implications in understanding the ligand conformational selection during complexation. Also, changes in availability to H-bonding concerning the different substituents and water can reflect on effects at complexation free energy due to variation at the desolvation energetic costs. Based on these results, we expect the obtained structural data provide systemic framework for rational chemical modification of STA-analogs.
Assuntos
Modelos Teóricos , Água , Estaurosporina/farmacologia , Ligantes , Água/química , Conformação Molecular , Termodinâmica , Simulação de Dinâmica MolecularRESUMO
One of the main concerns of Anfinsen was to reveal the connection between the amino-acid sequence and their biologically active conformation. This search gave rise to two crucial questions in structural biology, namely, why the proteins fold and how a sequence encodes its folding. As to the why, he proposes a plausible answer, namely, the thermodynamic hypothesis. As to the how, this remains an unsolved challenge. Consequently, the protein folding problem is examined here from a new perspective, namely, as an 'analytic whole'. Conceiving the protein folding in this way enabled us to (i) examine in detail why the force-field-based approaches have failed, among other purposes, in their ability to predict the three-dimensional structure of a protein accurately; (ii) propose how to redefine them to prevent these shortcomings, and (iii) conjecture on the origin of the state-of-the-art numerical-methods success to predict the tridimensional structure of proteins accurately.
Assuntos
Dobramento de Proteína , Proteínas , Proteínas/química , Conformação Molecular , Sequência de Aminoácidos , Termodinâmica , Conformação ProteicaRESUMO
For biomedical applications, gelatin is usually modified with methacryloyl groups to obtain gelatin methacryloyl (GelMA), which can be crosslinked by a radical reaction induced by low wavelength light to form mechanically stable hydrogels. The potential of GelMA hydrogels for tissue engineering has been well established, however, one of the main disadvantages of mammalian-origin gelatins is that their sol-gel transitions are close to room temperature, resulting in significant variations in viscosity that can be a problem for biofabrication applications. For these applications, cold-water fish-derived gelatins, such as salmon gelatin, are a good alternative due to their lower viscosity, viscoelastic and mechanical properties, as well as lower sol-gel transition temperatures, when compared with mammalian gelatins. However, information regarding GelMA (with special focus on salmon GelMA as a model for cold-water species) molecular conformation and the effect of pH prior to crosslinking, which is key for fabrication purposes since it will determine final hydrogel's structure, remains scarce. The aim of this work is to characterize salmon gelatin (SGel) and salmon methacryloyl gelatin (SGelMA) molecular configuration at two different acidic pHs (3.6 and 4.8) and to compare them to commercial porcine gelatin (PGel) and methacryloyl porcine gelatin (PGelMA), usually used for biomedical applications. Specifically, we evaluated gelatin and GelMA samples' molecular weight, isoelectric point (IEP), their molecular configuration by circular dichroism (CD), and determined their rheological and thermophysical properties. Results showed that functionalization affected gelatin molecular weight and IEP. Additionally, functionalization and pH affected gelatin molecular structure and rheological and thermal properties. Interestingly, the SGel and SGelMA molecular structure was more sensitive to pH changes, showing differences in gelation temperatures and triple helix formation than PGelMA. This work suggests that SGelMA presents high tunability as a biomaterial for biofabrication, highlighting the importance of a proper GelMA molecular configuration characterization prior to hydrogel fabrication.
Assuntos
Gelatina , Engenharia Tecidual , Animais , Gelatina/química , Temperatura de Transição , Viscosidade , Suspensões , Engenharia Tecidual/métodos , Metacrilatos/química , Salmão , Hidrogéis/química , Conformação Molecular , Água , MamíferosRESUMO
The use of quantum-based NMR methods to complement and guide the connectivity and stereochemical assignment of natural and unnatural products has grown enormously. One of the unsolved problems is related to the improper calculation of the conformational landscape of flexible molecules that have functional groups capable of generating a complex network of intramolecular H-bonding (IHB) interactions. Here the authors present MESSI (Multi-Ensemble Strategy for Structural Identification), a method inspired by the wisdom of the crowd theory that breaks with the traditional mono-ensemble approach. By including independent mappings of selected artificially manipulated ensembles, MESSI greatly improves the sense of the assignment by neutralizing potential energy biases.
Assuntos
Hidrogênio , Teoria Quântica , Modelos Moleculares , Conformação Molecular , Espectroscopia de Ressonância Magnética , Ligação de HidrogênioRESUMO
Previous studies propose that genetic mutations and post-translational modifications in protein crystallins promote protein aggregation and are considered significant risk factors for cataract formation. The ßB2-crystallin (HßB2C) forms a high proportion of proteins in the human eye lens. Different congenital mutations and post-translational deamidations in ßB2-crystallin have been reported and linked to cataract formation. In this work, we employed extensive all-atom molecular dynamics simulations to evaluate the conformational stability of deamidated and mutated HßB2C. Our results show critical changes in the protein surface and its native contacts due to a modification in the conformational equilibrium of these proteins. The double deamidated (Q70E/Q162E) and single deamidated (Q70E) impact the well compact conformation of the HßB2C. These post-translational modifications allow the exposure of the protein hydrophobic interface, which lead to the exposure of electronegative residues. On the other hand, our mutational studies showed that the S143F mutation modifies the hydrogen-bond network of an antiparallel ß-sheet, unfolding the C-terminal domain. Interestingly, the chain termination mutation (Q155X) does not unfold the N-terminal domain. However, the resultant conformation is more compact and avoids the exposure of the hydrophobic interface. Our results provide valuable information about the first steps of HßB2C unfolding in the presence of deamidated amino acids that have been reported to appear during aging. The findings reported in this work are essential for the general knowledge of the initial steps in the cataract formation mechanism, which may be helpful for the further development of molecules with pharmacological potential against cataract disease.
Assuntos
Catarata , Cristalino , beta-Cristalinas , Humanos , Cristalino/química , Conformação Molecular , Mutação , beta-Cristalinas/metabolismoRESUMO
In the search for new glycosidase inhibitors, a set of benzyl ß-D-Gal-S-(1â4)-3-deoxy-4-thio-α-D-hexopyranosides was synthesized. Diverse configurations were installed at C-2 and C-4 of the glucose residue. The benzyl glycosidic group was kept intact or substituted by an electron-donating or electron-withdrawing group that could also participate in hydrogen bonding. All thiodisaccharides were found to be inhibitors of E. coli ß-galactosidase. In general, benzyl thiodisaccharides were better inhibitors than those substituted (NO2 or NH2) on the benzyl ring. Thiodisaccharides containing a hexopyranoside, instead of a pentopyranoside, showed a weaker inhibitory activity, except for those having the α-D-xylo configuration, which exhibited inhibition constants of the same order of magnitude. These and previous results indicated that the inhibition process by thiodisaccharides is strongly dependent on the configuration of the 3-deoxy-4-thiopyranoside, as well as its substitution pattern (such as the presence of a benzyl glycoside). The enzyme-inhibitor interaction during the hydrolysis process involves a conformational selection resulting from rotation around the thioglycosidic bond and the flexibility of the terminal six-membered ring. Thus, the mentioned structural features of the inhibitor could give rise to favorable ground state conformations for the interaction with the enzyme, similar to those found for selected thiodisaccharides in the bound state. These studies demonstrated that the performance of thiodisaccharides as enzyme inhibitors could be increased by selecting the appropriate configuration and substitution of the hexopyranoside replacing the glucose moiety of 4-thiolactose.
Assuntos
Escherichia coli , Glicosídeos , Escherichia coli/metabolismo , Conformação Molecular , Glicosídeos/farmacologia , beta-Galactosidase/metabolismo , Inibidores Enzimáticos/química , GlucoseRESUMO
In this work the diastereoisomers (2S) and (2R)-naringenin-6-C-ß-D-glucopyroside, isolated for the first time from Clitoria guianensis, were studied using the density functional theory. The frontier molecular orbitals and structural properties showed that the diastereoisomers exhibit the same energy gap 166.61 kcal mol-1 and structural properties different, where in the S diastereoisomer, the bond length between the chiral carbon and the phenolic group is greater (difference of 0.0126 Å). The HPLC data showed that the retention time of the S-diastereoisomer (16.7 min) is shorter than that of R, suggesting that the S compound is more polar than R. The HPLC results corroborates with the molecular electrostatic potential which showed that in the S configuration, the electronegative density was more intense overall, particularly in the glucose molecule. The reactivity indices showed that the diastereoisomers are good electrophiles and reactive species. Finally, the absolute configuration of the diastereoisomers were determined using electronic circular dichroism (ECD) spectroscopy and the theoretical spectra were similar to the experimental. METHODS : All calculations of Density Functional Theory (DFT) and Time-Dependent Density Functional Theory (TDDFT) were performed using the program Gaussian 09 and the structures of the diastereoisomers were generated and analyzed using the GaussView program. The optimization and vibrational frequency calculations were performed using the functional CAM-B3LYP and 6-311 + + G(2d,2p) basis set. Conformational searches were performed for R configuration, by molecular mechanics using the MM + , MMFF, and OPLS05 force fields; the entire molecular mechanics simulation was performed using the Maestro/MacroModel software. The calculations for the simulations of the ECD spectra were performed for the eight lowest energy conformers obtained in the geometric optimization step, and the TDDFT at the CAM-B3LYP/6-311 + + G(2d,2p) theory level used. The effects of methanol and chloroform were calculated using the SMD implicit solvent model.
Assuntos
Clitoria , Estrutura Molecular , Dicroísmo Circular , Conformação MolecularRESUMO
Experimental and theoretical data have been revisited to shed light onto the aspects of hydration and chain expansion of pectic acid (galacturonan) upon charging. The prediction of the variation of the number of solvation water molecules between the two limit ionization states from theoretical calculations was confirmed to a very high accuracy by the corresponding number evaluated form dilatometric measurements. The relevance of hydration to the mechanism of bonding of calcium ions by sodium pectate is discussed. Characterization of polymer expansion has been obtained by calculating the values of the characteristic ratio and/or the persistence length on the respective populations and comparing the theoretical predictions with experimental data. The results show that a charged chain in typical conditions of ionic strength is more expanded than its neutral counterpart, whereas the ideal limit (31 and 21) helical conformations in the uncharged and totally charged conditions, respectively, share the same value of the linear advance of the helical repeat, when the ionic strength tends to infinite. Total divergence between theoretical predictions and experimental evidence rules out the possibility that carboxylate charge reduction by protonation and by methyl esterification are equivalent in determining the solution behavior of galacturonan.