Hybrid quantum mechanics/molecular mechanics (QM/MM) simulations have become an essential tool in computational chemistry, particularly for analyzing complex biological and condensed phase systems. Building on this foundation, our work presents a novel implementation of the Gaussian Electrostatic Model (GEM), a polarizable density-based force field, within the QM/MM framework. This advancement provides seamless integration, enabling efficient and optimized QM/GEM calculations in a single step using the LICHEM Code.
View Article and Find Full Text PDFHorseradish peroxidase (HRP) is an enzyme that oxidizes pollutants from wastewater. A previous report indicated that peroxidases can have an enhancement in initial enzymatic activity in an aqueous solution of 0.26 M 1-ethyl-3-methylimidazolium ethyl sulfate ([EMIm][EtSO]) at neutral pH.
View Article and Find Full Text PDFQM/MM methods have been used to study electronic structure properties and chemical reactivity in complex molecular systems where direct electronic structure calculations are not feasible. In our previous work, we showed that non-polarizable force fields, by design, describe intermolecular interactions through pairwise interactions, overlooking many-body interactions involving three or more particles. In contrast, polarizable force fields account partially for many-body effects through polarization, but still handle van der Waals and permanent electrostatic interactions pairwise.
View Article and Find Full Text PDFUnderstanding cooperativity and frustration is crucial for studying biological processes such as molecular recognition and protein aggregation. Force fields have been extensively utilized to explore cooperativity in the formation of protein secondary structures and self-assembled systems. Multiple studies have demonstrated that polarizable force fields provide more accurate descriptions of this phenomenon compared to fixed-charge pairwise nonpolarizable force fields, thanks to the incorporation of polarization effects.
View Article and Find Full Text PDFThe main protease (Mpro) of SARS-CoV-2 is an essential enzyme for the replication of the virus causing the COVID-19 pandemic. Because there is no known homologue in humans, it has been proposed as a primary target for antiviral drug development. Here, we explore the potential of five acrylamide-based molecules as possible covalent inhibitors, leading to target MPro by docking, followed by polarizable molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) calculations.
View Article and Find Full Text PDFVarious photosynthetic organisms have evolved to absorb light in different regions of the visible light spectrum, thus adapting to the various lighting conditions available on Earth. While most of these autotrophic organisms absorb wavelengths around the 700-800 nm region, some are capable of red-shifted absorptions above this range, but none as remarkably as whose main absorption is observed at 1015 nm, approximately 220 nm (0.34 eV) lower in energy than their main constituent pigments, BChl-, whose main absorption is observed at 795 nm.
View Article and Find Full Text PDFPrions have been linked to neurodegenerative diseases that affect various species of mammals including humans. The prion protein, located mainly in neurons, is believed to play the role of metal ion transporter. High levels of copper ions have been related to structural changes.
View Article and Find Full Text PDFWiley Interdiscip Rev Comput Mol Sci
January 2021
Quantum mechanics/molecular mechanics (QM/MM) simulations are a popular approach to study various features of large systems. A common application of QM/MM calculations is in the investigation of reaction mechanisms in condensed-phase and biological systems. The combination of QM and MM methods to represent a system gives rise to several challenges that need to be addressed.
View Article and Find Full Text PDFThe origin of non-additivity in hydrogen bonds (H-bonds), usually termed as H-bond cooperativity, is investigated in H-bonded linear chains. It is shown that H-bond cooperativity originates solely from classical electrostatics. The latter is corroborated by comparing the H-bond cooperativity in infinitely-long H-bonded hydrogen cyanide, 4-pyridone and formamide chains, assessed using density functional theory (DFT), against the strengthening of the dipole-dipole interaction upon the formation of an infinite chain of effective point-dipoles.
View Article and Find Full Text PDFSome health disturbances like neurodegenerative diseases are associated to the presence of amyloids. GNNQQNY and NNQQNY peptides are considered as prototypical examples for studying the formation of amyloids. These exhibit quite different aggregation behaviors despite they solely differ in size by one residue.
View Article and Find Full Text PDF