Well-tempered Metadynamics Simulations Research Articles

Advanced porous covalent organic framework (COF) materials for the capture of alizarin dye and its derivatives from the aquatic environment

The effective removal of dye pollutants from water and wastewater is a key environmental challenge. The present study is developed to investigate alizarin (ALI) dye and its derivations, including Alizarin blue (ABL), Alizarin purpurin (APU), Quinalizarin (AQU), Alizarin cyanin (ACY), and Alizarin Red S (ARS) removal process from water and wastewater sources, using covalent organic frameworks (COFs) nanoadsorbents. Herein, we explore the process of how dye molecules are absorbed onto COFs with precise supramolecular structures. The molecular dynamics (MD) and well-tempered metadynamics (WTMtD) simulations are used to investigate this process in aqueous solution. From the results obtained, it is clear that the intermolecular van der Waals (vdw) and π-π interactions have a significant role on accelerating the interaction between dye molecules and the COF nanostructures. This ultimately leads to the creation of a stable dye-COF complex. The dye-adsorbent average interaction energy value reaches around APU-COF1=−604.34, AQU-COF1=-515.25, ABL-COF1=−504.74, ALI-COF1=−489.48, ARS-COF1=−475.81, ACY-COF1=−273.82, AQU-COF2=−459.76, ALI-COF2=−451.46, ABL-COF2=−405.90, APU-COF2=−367.55, ACY-COF2=−287.89, ARS-COF2=−210.63 kJ/mol for dye/COF1 and dye/COF2 complexes, respectively. The primary interaction between dye and COFs is attributed to the Lennard-Jones term, resulting from the formation of a strong π-π interaction between the dye molecules and the surface of the adsorbent. Overall, our simulations confirmed that the COF1 nanostructure is more effective than the COF2 nanostructure in removing alizarin dye and its derivatives. In this study, not only the performance of two COFs in removing alizarin dye and its derivatives has been compared, but also the possibility of removing alizarin dye and its derivatives with both COFs has been examined.

Computational screening combined with well-tempered metadynamics simulations identifies potential TMPRSS2 inhibitors

Type-II transmembrane serine proteases are effective pharmacological targets for host defence against viral entry and in certain cancer cell progressions. These serine proteases cleave viral spike proteins to expose the fusion peptide for cell entry, which is essential to the life cycle of the virus. TMPRSS2 inhibitors can also fight against respiratory viruses that employ them for cell entry. Our study combining virtual screening, all-atom molecular dynamics, and well-tempered metadynamics simulation identifies vicenin-2, neohesperidin, naringin, and rhoifolin as promising TMPRSS2 antagonists. The binding energies obtained are − 16.3, − 15.4, − 13.6, and − 13.8 kcal/mol for vicenin-2, neohesperidin, naringin, and rhoifolin respectively. The RMSD, RMSF, PCA, DCCM, and binding free energy profiles also correlate with the stable binding of these ligands at the active site of TMPRSS2. The study reveals that these molecules could be promising lead molecules for combating future outbreaks of coronavirus and other respiratory viruses.

Inhibition of the Early-Stage Cross-Amyloid Aggregation of Amyloid-β and IAPP via EGCG: Insights from Molecular Dynamics Simulations.

Amyloid-β (Aβ) and islet amyloid polypeptide (IAPP) are small peptides that have the potential to not only self-assemble but also cross-assemble and form cytotoxic amyloid aggregates. Recently, we experimentally investigated the nature of Aβ-IAPP coaggregation and its inhibition by small polyphenolic molecules. Notably, we found that epigallocatechin gallate (EGCG) had the ability to reduce heteroaggregate formation. However, the precise molecular mechanism behind the reduction of heteroaggregates remains unclear. In this study, the dimerization processes of Aβ40 and IAPP peptides with and without EGCG were characterized by the enhanced sampling technique. Our results showed that these amyloid peptides exhibited a tendency to form a stable heterodimer, which represented the first step toward coaggregation. Furthermore, we also found that the EGCG regulated the dimerization process. In the presence of EGCG, well-tempered metadynamics simulation indicated a notable shift in the bound state toward a greater center of mass (COM) distance. Additionally, the presence of EGCG led to a significant increase in the free energy barrier height (∼15k B T) along the COM distance, and we observed a transition state between the bound and unbound states. Our findings also unveiled that the EGCG formed a greater number of hydrogen bonds with Aβ40, effectively obstructing the dimer formation. In addition, we carried out microseconds of all-atom conventional molecular dynamics (cMD) simulations to investigate the formation of both hetero- and homo-oligomer states by these peptides. MD simulations illustrated that EGCG played a significant role in preventing oligomer formation by reducing the content of β-sheets in the peptide. Collectively, our results offered valuable insight into the mechanism of cross-amyloid aggregation between Aβ40 and IAPP and the inhibition effect of EGCG on the heteroaggregation process.

Engineered nanoparticles as Selinexor drug delivery systems across the cell membrane and related signaling pathways in cancer cells

In the present work, molecular dynamics simulation is applied to evaluate the drug carrier efficiency of graphene oxide nanoflake (GONF) for loading of Selinexor (SXR) drug as well as the drug delivery by 2D material through the membrane in aqueous solution. In addition, to investigate the adsorption and penetration of drug-nanocarrier complex into the cell membrane, well-tempered metadynamics simulations and steered molecular dynamics (SMD) simulations were performed. Based on the obtained results, it is evident that intermolecular hydrogen bonds (HBs) and π-π interactions play a significant role in expediting the interaction between drug molecules and the graphene oxide (GO) nanosheet, ultimately resulting in the formation of a stable SXR-GO complex. The Lennard-Jones (L-J) energy value for the interaction of SXR with GONF is calculated to be approximately −98.85 kJ/mol. In the SXR-GONF complex system, the dominant interaction between SXR and GONF is attributed to the L-J term, resulting from the formation of a strong π−π interaction between the drug molecules and the substrate surface. Moreover, our simulations show by decreasing the distance of GONF with respect to cell membrane, the interaction energy of GONF-membrane significantly decrease to −1500 kJ/mol resulting in fast diffusion of SXR-GONF complex toward the bilayer surface that is favored opening the way to natural drug nanocapsule.

Significantly enhanced performance for phenol compounds removal by MOF-5 nano-composite via its surface modification

The present study is focused on the use of cubic metal-organic frameworks-5 (MOF-5) and its functionalized form in the removal of phenolic pollutants by molecular dynamics (MD) and Well-tempered metadynamics (WTMD) simulation methods. It was found that the adsorption mechanism of MOF-5s/phenolic compounds is mostly due to the van der Waals and π–π interactions. However, electrostatic and hydrogen bond (HB) interactions also play a significant role in removing phenolic pollutants by MOF-5 and its functionalized form. The results show that the fluorine functional group (F-MOF-5) increases the adsorption capacity of phenol compounds on the adsorbent surface. By functionalizing the MOF-5 with a methyl functional group (CH3-MOF-5), the adsorption strength decreases. The WTMD calculation confirmed that at the most stable state, the free energy (FE) value of system II (the most stable system in functionalized systems with –F functional group) is about −289.528 kJ mol−1. This value is ~5.781 and 35.514 kJ mol−1 more negative than the FE of the I and III systems (the most stable systems in the pristine and CH3-MOF-5/pollutant systems, respectively). Altogether, the results indicate that F-MOF-5 can be considered a more suitable adsorbent than MOF-5 and CH3-MOF-5 for phenolic pollutants removal from the environment for more assessment.

Ionic liquids and Graphene: The ultimate combination for High-Performance supercapacitors

Ionic liquids (ILs), due to their tunable propensities and properties, are discussed as a new class of solvents for scientific and industrial applications. Developing novel processes and applying ionic liquids further accelerate the growth of our fundamental and engineering understanding.In this work, the structure of four ILs, 1-ethyl-3-methylimidazolium tetrafluoroborate formate (Emim+BF4−)/chloride (Emim+Cl−)/Nitrate (Emim+NO3−), and thiocyanate (Emim+SCN−) is reported. Then via employing classical molecular dynamics and well-tempered metadynamics simulations their dynamics are investigated. The dynamics of ILs are represented by studying mean squared displacements. Additionally, the increased size and weight of anions interacting with the Emim+ cation cause the dynamics to become slow. Therefore, due to the higher molecular weight, the BF4− anion has lower diffusion than any other. Indeed, radial distribution functions between cations and anions reveal the structural arrangement in ILs. On the other hand, it is unexpectedly found that a certain amount of ILs are distributed near the electrode surface. Such phenomena occur because of the relatively strong ionic interactions between cations and anions; also, the π–π stacking between cations and graphene caused their interactions to become stronger. Furthermore, the atoms in molecules (AIM) analyses identified several interatomic interactions between the ILs and graphene. The free energy values for the ILs at their global minima are about ∼ -271.79, −267.93, −251.37 kJ mol−1 and ∼ -241.98 kJ mol−1 in systems Emim+/Cl−, Emim//NO3−, Emim+/SCN−, and Emim+/BF4−, respectively. Our findings represent a crucial step towards interpreting the complex interactions and nanostructures at the electrolyte–electrode interface, which could be beneficial for designing high-performance supercapacitors in energy storage using carbon electrodes and IL electrolytes in experiments.

The strategy of three-dimensional Covalent Organic Frameworks to exclude dye contaminants in aqueous solutions

The development of effective solid adsorbents is necessary to remove organic pollutants from wastewater in response to the global water crisis. The present study investigates the adsorption behavior of dye pollutants on a 1, 3, 6, 8-tetrakis (4formylphenyl) pyrene (TFPPy) for the self-assembly of an innovative three-dimensional (3D) covalent organic frameworks (COFs). Molecular dynamics (MD) and the well-tempered metadynamics simulation were used to evaluate the ability of COFs to adsorb seven types of dye pollutants (DPs). These include methyl orange (MO), gentian violet (GV), brilliant green (BG), tolonium chloride (TC), methylene blue (MB), anionic congo red (CR), and rhodamine B (Rh B). To examine this process, several descriptors are calculated from the simulation trajectories, including interaction energies, root-mean-square deviation, radial distribution function, solvent-accessible surface area, mean square displacement, and the number of hydrogen bonds to figure out the process. After analyzing the data, it was found that the MO/COFs system displayed exceptional stability. It was confirmed by a significant value of −448.54 KJmol−1 for total interaction energy. The density profile of dye molecules in the MO and Rh B-COFs systems exhibits the highest peaks at approximately 325.79, and 26.75 kg/m³, respectively. The free energy values for the DPs/COFs complexes at their global minima reached about, −512.127, −435.603, −381.244, −300.313, and −283.230 kJmol-1 for MO/GV/BG/TC/ and Rh B with the COFs, respectively. The nature of the design of TFFPy-COF’s unique, ultrafast kinetics, high adsorption capacity, and good reusability endow them with great potential for removing dyes in an aqueous environment.

Investigating the Aryl Hydrocarbon Receptor Agonist/Antagonist Conformational Switch Using Well-Tempered Metadynamics Simulations.

The aryl hydrocarbon receptor (AhR) is a ligand-dependent transcription factor that mediates biological signals to control various complicated cellular functions. It plays a crucial role in environmental sensing and xenobiotic metabolism. Dysregulation of AhR is associated with health concerns, including cancer and immune system disorders. Upon binding to AhR ligands, AhR, along with heat shock protein 90 and other partner proteins undergoes a transformation in the nucleus, heterodimerizes with the aryl hydrocarbon receptor nuclear translocator (ARNT), and mediates numerous biological functions by inducing the transcription of various AhR-responsive genes. In this manuscript, the 3-dimensional structure of the entire human AhR is obtained using an artificial intelligence tool, and molecular dynamics (MD) simulations are performed to study different structural conformations. These conformations provide insights into the protein's function and movement in response to ligand binding. Understanding the dynamic behavior of AhR will contribute to the development of targeted therapies for associated health conditions. Therefore, we employ well-tempered metadynamics (WTE-metaD) simulations to explore the conformational landscape of AhR and obtain a better understanding of its functional behavior. Our computational results are in excellent agreement with previous experimental findings, revealing the closed and open states of helix α1 in the basic helix-loop-helix (bHLH domain) in the cytoplasm at the atomic level. We also predict the inactive form of AhR and identify Arginine 42 as a key residue that regulates switching between closed and open conformations in existing AhR modulators.

In silico drug derivatives for KRAS-G12D: free-energy surfaces in aqueous solution by well-tempered metadynamics simulations

KRAS oncogenes are the largest family of mutated RAS isoforms, participating in about 30% of all cancers. Due to their paramount medical importance, enormous effort is being devoted to the development of inhibitors using clinical tests, wet-lab experiments and drug design, being this a preliminary step in the process of creating new drugs, prior to synthesis and clinical testing. One central aspect in the development of new drugs is the characterisation of all species that can be used for treatment. In this aim, we propose a computational framework based on combined all-atom molecular dynamics and metadynamics simulations in order to accurately access the most stable conformational variants for several derivatives of a recently proposed small-molecule, called DBD15-21-22. Free energy calculations are essential to unveil mechanisms at the atomic scale like binding affinities or dynamics of stable states. Considering specific atom-atom distances and torsional angles as reliable reaction coordinates we have obtained free-energy landscapes by well-tempered metadynamics simulations, revealing local and global minima of the free-energy hypersurface. We have observed that a variety of stable states together with transition states are clearly detected depending on the particular species, leading to predictions on the behaviour of such compounds in ionic aqueous solution.

The Structural Role of N170 in Substrate-Assisted Deacylation in KPC-2 β-Lactamase.

The amino acid substitutions in Klebsiella pneumoniae carbapenemase 2 (KPC-2) that have arisen in the clinic are observed to lead to the development of resistance to ceftazidime-avibactam, a preferred treatment for KPC bearing Gram-negative bacteria. Specific substitutions in the omega loop (R164-D179) result in changes in the structure and function of the enzyme, leading to alterations in substrate specificity, decreased stability, and more recently observed, increased resistance to ceftazidime/avibactam. Using accelerated rare-event sampling well-tempered metadynamics simulations, we explored in detail the structural role of R164 and D179 variants that are described to confer ceftazidime/avibactam resistance. The buried conformation of D179 substitutions produce a pronounced structural disorder in the omega loop - more than R164 mutants, where the crystallographic omega loop structure remains mostly intact. Our findings also reveal that the conformation of N170 plays an underappreciated role impacting drug binding and restricting deacylation. The results further support the hypothesis that KPC-2 D179 variants employ substrate-assisted catalysis for ceftazidime hydrolysis, involving the ring amine of the aminothiazole group to promote deacylation and catalytic turnover. Moreover, the shift in the WT conformation of N170 contributes to reduced deacylation and an altered spectrum of enzymatic activity.

The potential of some functional group compounds substituted 8-Manzamine A as RSK1 inhibitors: molecular docking and molecular dynamics simulations

Cancer, an incurable global disease, demands urgent anti-cancer drug development. Marine alkaloids like Manzamine and its derivatives show promise as RSK inhibitors against cancer cell invasion. Replacing the hydrogen at the 8-position of Manzamine A with a hydroxyl group has been shown to significantly enhance its biological activity. In this article, we designed various functional group compounds (A1–A21) substituted 8-Manzamine A by docking, MM-GBSA, molecular dynamics (MD) simulation, and well-tempered metadynamics (WT-MetaD) simulations to evaluate their potential as RSK1 inhibitors. Ligands A1–A21 were docked in the RSK1 N-terminal kinase domain (PDB ID: 2Z7Q) using the Glide module. The calculation of binding energy was performed using Prime MM-GB/SA, while MD simulations were conducted with the Desmond module of Schrodinger suite 2023. Compound A5 exhibits the highest G-score (−7.01) compared to 8-Hydroxymanzamine A (−6.08). Additionally, compounds A6, A10, A12, A17, A11, A4, and A13 demonstrate increased activity against RSK1 when compared to both 8-Hydroxymanzamine A and Manzamine A. Residues LEU68, VAL76, LEU141, PHE143, LEU144, PHE150, ASP148, GLU191, and LEU194 of RSK1 protein play a key role in binding with ligands. An MD simulation of Compound A5 was carried out to explore the dynamic interactions within the protein-ligand complex. Furthermore, WT-MetaD simulations validated the docking study results and identified the most energetically favored conformations for the A5/RSK1 complex. Ligands A5, A6, A10, A12, A17, A11, A4, and A13, featuring diverse functional groups and good Glide scores, may have the potential for significant RSK1 activity and merit further development. Communicated by Ramaswamy H. Sarma

Stable and accurate atomistic simulations of flexible molecules using conformationally generalisable machine learned potentials.

Computational simulation methods based on machine learned potentials (MLPs) promise to revolutionise shape prediction of flexible molecules in solution, but their widespread adoption has been limited by the way in which training data is generated. Here, we present an approach which allows the key conformational degrees of freedom to be properly represented in reference molecular datasets. MLPs trained on these datasets using a global descriptor scheme are generalisable in conformational space, providing quantum chemical accuracy for all conformers. These MLPs are capable of propagating long, stable molecular dynamics trajectories, an attribute that has remained a challenge. We deploy the MLPs in obtaining converged conformational free energy surfaces for flexible molecules via well-tempered metadynamics simulations; this approach provides a hitherto inaccessible route to accurately computing the structural, dynamical and thermodynamical properties of a wide variety of flexible molecular systems. It is further demonstrated that MLPs must be trained on reference datasets with complete coverage of conformational space, including in barrier regions, to achieve stable molecular dynamics trajectories.

Can local heating and molecular crowders disintegrate amyloid aggregates?

The present study employs a blend of molecular dynamics simulations and a theoretical model to explore the potential disintegration mechanism of a matured Aβ octamer, aiming to offer a strategy to combat Alzheimer's disease. We investigate local heating and crowding effects on Aβ disintegration by selectively heating key Aβ segments and varying the concentration of sodium dodecyl sulphate (SDS), respectively. Despite initiation of disruption, Aβ aggregates resist complete disintegration during local heating due to rapid thermal energy distribution to the surrounding water. Conversely, although SDS molecules effectively inhibit Aβ aggregation at higher concentration through micelle formation, they fail to completely disintegrate the aggregate due to the exceedingly high energy barrier. To address the sampling challenge posed by the formidable energy barrier, we have performed well-tempered metadynamics simulations. Simulations reveal a multi-step disintegration mechanism for the Aβ octamer, suggesting a probable sequence: octamer → pentamer/hexamer ⇌ tetramer → monomer, with a rate-determining step constituting 45 kJ mol-1 barrier during the octamer to pentamer/hexamer transition. Additionally, we have proposed a novel two-state mean-field model based on Ising spins that offers an insight into the kinetics of the Aβ growth process and external perturbation effects on disintegration. Thus, the current simulation study, coupled with the newly introduced mean-field model, offers an insight into the detailed mechanisms underlying the Aβ aggregation process, guiding potential strategies for effective disintegration of Aβ aggregates.

Assessing the role of residue Phe108 of cytochrome P450 3A4 in allosteric effects of midazolam metabolism.

Cytochrome P450 3A4 (CYP3A4) is involved in the metabolism of more drugs in clinical use than any other xenobiotic-metabolizing enzyme. CYP3A4-mediated drug metabolism is usually allosterically modulated by substrate concentration (homotropic allostery) and other drugs (heterotropic allostery), exhibiting unusual kinetic profiles and regiospecific metabolism. Recent studies suggest that residue Phe108 (F108) of CYP3A4 may have an important role in drug metabolism. In this work, residue mutations were coupled with well-tempered metadynamics simulations to assess the importance of F108 in the allosteric effects of midazolam metabolism. Comparing the simulation results of the wild-type and mutation systems, we identify that the π-π interaction and steric effect between the F108 side chain and midazolam is favorable for the stable binding of substrate in the active site. F108 also plays an important role in the transition of substrate binding mode, which mainly induces the transition of substrate binding mode by forming π-π interactions with multiple aromatic rings of the substrate. Moreover, the side chain of F108 is closely related to the radius and depth of the 2a and 2f channels, and F108 may further regulate drug metabolism by affecting the pathway, orientation, or time of substrate entry into the CYP3A4 active site or product egress from the active site. Altogether, we suggest that F108 affects drug metabolism and regulatory mechanisms by affecting substrate binding stability, binding mode transition, and channel characteristics of CYP3A4. Our findings could promote the understanding of complicated allosteric mechanisms in CYP3A4-mediated drug metabolism, and the knowledge could be used for drug development and disease treatment.

Attempting Well-Tempered Funnel Metadynamics Simulations for the Evaluation of the Binding Kinetics of Methionine Aminopeptidase-II Inhibitors.

Understanding the unbinding kinetics of protein-ligand complexes is considered a significant approach for the design of ligands with desired specificity and safety. In recent years, enhanced sampling methods have emerged as effective tools for studying the unbinding kinetics of protein-ligand complexes at the atomistic level. MetAP-II is a target for the treatment of cancer for which not a single effective drug is available yet. The identification of the dissociation rate of ligands from the complexes often serves as a better predictor for in vivo efficacy than the ligands' binding affinity. Here, funnel-based restraint well-tempered metadynamics simulations were applied to predict the residence time of two ligands bound to MetAP-II, along with the ligand association and dissociation mechanism involving the identification of the binding hotspot during ligand egress. The ligand-egressing route revealed by metadynamics simulations also correlated with the identified pathways from the CAVER analysis and by the enhanced sampling simulation using PLUMED. Ligand 1 formed a strong H-bond interaction with GLU364 estimating a higher residence time of 28.22 ± 5.29 ns in contrast to ligand 2 with a residence time of 19.05 ± 3.58 ns, which easily dissociated from the binding pocket of MetAP-II. The results obtained from the simulations were consistent to reveal ligand 1 being superior to ligand 2; however, the experimental data related to residence time were close for both ligands, and no kinetic data were available for ligand 2. The current study could be considered the first attempt to apply an enhanced sampling method for the evaluation of the binding kinetics and thermodynamics of two different classes of ligands to a binuclear metalloprotein.

Molecular mechanism of calcitriol enhances membrane water permeability

Helicobacter pylori (H. pylori) exhibits a unique membrane lipid composition, including dimyristoyl phosphatidylethanolamine (DMPE) and cholesterol, unlike other Gram-negative bacteria. Calcitriol has antimicrobial activity against H. pylori, but cholesterol enhances antibiotics resistance in H. pylori. This study explored the changes in membrane structure and the molecular mechanisms of cholesterol/calcitriol translocation using well-tempered metadynamics (WT-MetaD) simulations and microsecond conventional molecular dynamics (CMD) simulations. Calcitriol facilitated water transport across the membrane, while cholesterol had the opposite effect. The differing effects might result from the tail 25-hydroxyl group and a wider range of orientations of calcitriol in the DMPE/dimyristoyl phosphatidylglycerol (DMPG) (3:1) membrane. Calcitriol moves across the bilayer center without changing its orientation along the membrane Z-axis, becomes parallel to the membrane surface at the membrane-water interface, and then rotates approximately 90° in this interface. The translocation mechanism of calcitriol is quite different from the flip-flop of cholesterol. Moreover, calcitriol crossed from one layer to another more easily than cholesterol, causing successive perturbations to the hydrophobic core and increasing water permeation. These results improve our understanding of the relationship between cholesterol/calcitriol concentrations and the lipid bilayer structure and the role of lipid composition in water permeation.

Natural protein engineering in the Ω-loop: the role of Y221 in ceftazidime and ceftolozane resistance in Pseudomonas-derived cephalosporinase.

A wide variety of clinically observed single amino acid substitutions in the Ω-loop region have been associated with increased minimum inhibitory concentrations and resistance to ceftazidime (CAZ) and ceftolozane (TOL) in Pseudomonas-derived cephalosporinase and other class C β-lactamases. Herein, we demonstrate the naturally occurring tyrosine to histidine substitution of amino acid 221 (Y221H) in Pseudomonas-derived cephalosporinase (PDC) enables CAZ and TOL hydrolysis, leading to similar kinetic profiles (k cat = 2.3 ± 0.2 µM and 2.6 ± 0.1 µM, respectively). Mass spectrometry of PDC-3 establishes the formation of stable adducts consistent with the formation of an acyl enzyme complex, while spectra of E219K (a well-characterized, CAZ- and TOL-resistant comparator) and Y221H are consistent with more rapid turnover. Thermal denaturation experiments reveal decreased stability of the variants. Importantly, PDC-3, E219K, and Y221H are all inhibited by avibactam and the boronic acid transition state inhibitors (BATSIs) LP06 and S02030 with nanomolar IC50 values and the BATSIs stabilize all three enzymes. Crystal structures of PDC-3 and Y221H as apo enzymes and complexed with LP06 and S02030 (1.35-2.10 Å resolution) demonstrate ligand-induced conformational changes, including a significant shift in the position of the sidechain of residue 221 in Y221H (as predicted by enhanced sampling well-tempered metadynamics simulations) and extensive hydrogen bonding between the enzymes and BATSIs. The shift of residue 221 leads to the expansion of the active site pocket, and molecular docking suggests substrates orientate differently and make different intermolecular interactions in the enlarged active site compared to the wild-type enzyme.

Multifunctional Roles of Clathrate Hydrate Nanoreactors for CO2 Reduction.

In this study, we explore a possible platform for the CO2 reduction (CO2 R) in one of water's solid phases, namely clathrate hydrates (CHs), by ab initio molecular dynamics and well-tempered metadynamics simulations with periodic boundary conditions. We found that the stacked H2 O nanocages in CHs help to initialize CO2 R by increasing the electron-binding ability of CO2 . The substantial CO2 R processes are further influenced by the hydrogen bond networks in CHs. The first intermediate CO2 - in this process can be stabilized through cage structure reorganization into the H-bonded [CO2 - ⋅⋅⋅H-OHcage ] complex. Further cooperative structural dynamics enables the complex to convert into a vital transient [CO2 2- ⋅⋅⋅H-OHcage ] intermediate in a low-barrier disproportionation-like process. Such a highly reactive intermediate spontaneously triggers subsequent double proton transfer along its tethering H-bonds, finally converting it into HCOOH. These hydrogen-bonded nanoreactors feature multiple functions in facilitating CO2 R such as confining, tethering, H-bond catalyzing and proton pumping. Our findings have a general interest and extend the knowledge of CO2 R into porous aqueous systems.

Free-Energy Landscape and Rate Estimation of the Aromatic Ring Flips in Basic Pancreatic Trypsin Inhibitors Using Metadynamics.

Aromatic side chains (phenylalanine and tyrosine) of a protein flip by 180° around the Cβ-Cγ axis (χ2 dihedral of the side chain), producing two symmetry-equivalent states. The study of ring flip dynamics with nuclear magnetic resonance (NMR) experiments helps to understand local conformational fluctuations. Ring flips are categorized as slow (milliseconds and onward) or fast (nanoseconds to near milliseconds) based on timescales accessible to NMR experiments. In this study, we investigated the ability of the infrequent metadynamics approach to estimate the flip rate and discriminate between slow and fast ring flips for eight individual aromatic side chains (F4, Y10, Y21, F22, Y23, F33, Y35, and F45) of the basic pancreatic trypsin inhibitor. Well-tempered metadynamics simulations were performed to estimate the ring-flipping free-energy surfaces for all eight aromatic residues. The results indicate that χ2 as a standalone collective variable (CV) is not sufficient to obtain computationally consistent results. Inclusion of a complementary CV, such as χ1(Cα-Cβ), solved the problem for most residues and enabled us to classify fast and slow ring flips. This indicates the importance of librational motions in ring flips. Multiple pathways and mechanisms were observed for residues F4, Y10, and F22. Recrossing events were observed for residues F22 and F33, indicating a possible role of friction effects in ring flipping. The results demonstrate the successful application of infrequent metadynamics to estimate ring flip rates and identify certain limitations of the approach.

Effects of C-Terminal Lys-Arg Residue of AapA1 Protein on Toxicity and Structural Mechanism.

Previous experimental investigations have established the indispensability of the C-terminal Lys-Arg residues in the toxic activity of the AapA1 toxin protein. AapA1 is classified as a type I toxin-antitoxin (TA) bacterial toxin, and the precise impact of the C-terminal Lys-Arg residues on its structure and mechanism of action remains elusive. To address this knowledge gap, the present study employed molecular dynamics (MD) and enhanced sampling Well-tempered Two-dimensional Metadynamics (2D-MetaD) simulations to examine the behavior of the C-terminal Lys-Arg residues of truncated AapA1 toxin (AapA1-28) within the inner membrane of Escherichia coli. Specifically, the study focused on the elucidation of possible conformation states of AapA1-28 protein in POPE/POPG (3:1) bilayers and their interactions between the protein and POPE/POPG (3:1) bilayers. The findings of our investigation indicate that the AapA1-28 protein does not adopt a vertical orientation upon membrane insertion; rather, it assumes an angled conformation, with the side chain of Lys-23 directed toward the upper layer of the membrane. This non-transmembrane conformation of AapA1-28 protein impedes its ability to form pores within the membrane, resulting in reduced toxicity towards Escherichia coli. These results suggest that C-Terminal positively charged residues are essential for electrostatic binding to the negatively charged head group of bottom bilayer membrane, which stabilize the transmembrane conformation. These outcomes contribute to our comprehension of the impact of C-terminal charged residues on the structure and functionality of membrane-associated proteins, and provide an improved understanding of how protein sequence influences the antimicrobial effect.