Environmental remediation of per- and polyfluoroalkyl substances (PFAS) has become a significant research topic in recent years due to the fact that these materials are omnipresent, resistant to degradation and thus environmentally persistent. Unfortunately, they have also been shown to cause health concerns. PFAS are widely used in industrial applications and consumer products. Vitamin B12 (B12) has been identified as being catalytically active towards a variety of halogenated compounds such as PFAS. It has also been shown to be effective when using sulfide as a reducing agent for B12. This is promising as sulfide is readily available in the environment. However, there are many unknowns with respect to PFAS interactions with B12. These include the reaction mechanism and B12\'s specificity for PFAS with certain functionalization(s). In order to understand the specificity of B12 towards branched PFAS, we examined the atomistic interactions between B12 and eight different PFAS molecules using Density Functional Theory (B3LYP/cc-pVDZ). The PFAS test set included linear PFAS and their branched analogs, carboxylic acid and sulfonic acid headgroups, and aromatic and non-aromatic cyclic structures. Conformational analyses were carried out to determine the lowest energy configurations. This analysis showed that small chain PFAS such as perfluorobutanoic acid interact with the cobalt center of B12. Bulkier PFAS prefer to interact with the amine and carbonyl groups on the sidechains of the B12 ring system. Furthermore, computed complexation energies determined that, in general, branched PFAS (e.g. perfluoro-5-methylheptane sulfonic acid) interact more strongly than linear molecules (e.g. perfluorooctanesulfonic acid). Our results indicate that it may be possible to alter the interactions between B12 and PFAS by synthetically modifying the sidechains of the ring structure.

BACKGROUND: Molecularly imprinted polymers (MIPs) have promising applications as synthetic antibodies for protein and peptide recognition. A critical aspect of MIP design is the selection of functional monomers and their adequate proportions to achieve materials with high recognition capacity toward their targets. To contribute to this goal, we calibrated a molecular dynamics protocol to reproduce the experimental trends in peptide recognition of 13 pre-polymerization mixtures reported in the literature for the peptide toxin melittin.
METHODS: Three simulation conditions were tested for each mixture by changing the box size and the number of monomers and cross-linkers surrounding the template in a solvent-explicit environment. Fully atomistic MD simulations of 350 ns were conducted with the AMBER20 software, with ff19SB parameters for the peptide, gaff2 parameters for the monomers and cross-linkers, and the OPC water model. Template-monomer interaction energies under the LIE approach showed significant differences between high-affinity and low-affinity mixtures. Simulation systems containing 100 monomers plus cross-linkers in a cubic box of 90 Å3 successfully ranked the mixtures according to their experimental performance. Systems with higher monomer densities resulted in non-specific intermolecular contacts that could not account for the experimental trends in melittin recognition. The mixture with the best recognition capacity showed preferential binding to the 13-26-α-helix, suggesting a relevant role for this segment in melittin imprinting and recognition. Our findings provide insightful information to assist the computational design of molecularly imprinted materials with a validated protocol that can be easily extended to other templates.

This work systematically examines the interactions between a single argon atom and the edges and faces of cyclic H2O clusters containing three-five water molecules (Ar(H2O)n=3-5). Full geometry optimizations and subsequent harmonic vibrational frequency computations were performed using MP2 with a triple-ζ correlation consistent basis set augmented with diffuse functions on the heavy atoms (cc-pVTZ for H and aug-cc-pVTZ for O and Ar; denoted as haTZ). Optimized structures and harmonic vibrational frequencies were also obtained with the two-body-many-body (2b:Mb) and three-body-many-body (3b:Mb) techniques; here, high-level CCSD(T) computations capture up through the two-body or three-body contributions from the many-body expansion, respectively, while less demanding MP2 computations recover all higher-order contributions. Five unique stationary points have been identified in which Ar binds to the cyclic water trimer, along with four for (H2O)4 and three for (H2O)5. To the best of our knowledge, eleven of these twelve structures have been characterized here for the first time. Ar consistently binds more strongly to the faces than the edges of the cyclic (H2O)n clusters, by as much as a factor of two. The 3b:Mb electronic energies computed with the haTZ basis set indicate that Ar binds to the faces of the water clusters by at least 3 kJ mol-1 and by nearly 6 kJ mol-1 for one Ar(H2O)5 complex. An analysis of the interaction energies for the different binding motifs based on symmetry-adapted perturbation theory (SAPT) indicates that dispersion interactions are primarily responsible for the observed trends. The binding of a single Ar atom to a face of these cyclic water clusters can induce perturbations to the harmonic vibrational frequencies on the order of 5 cm-1 for some hydrogen-bonded OH stretching frequencies.

Accurate information about interactions between group I metals and nucleic acids is required to understand the roles these metals play in basic cellular functions, disease progression, and pharmaceuticals, as well as to aid the design of new energy storage materials and nucleic acid sensors that target metal contaminants, among other applications. From this perspective, this work generates a complete CCSD(T)/CBS data set of the binding energies for 64 complexes involving each group I metal (Li+, Na+, K+, Rb+, or Cs+) directly coordinated to various sites in each nucleic acid component (A, C, G, T, U, or dimethylphosphate). This data have otherwise been challenging to determine experimentally, with highly accurate information missing for many group I metal-nucleic acid combinations and no data available for the (charged) phosphate moiety. Subsequently, the performance of 61 DFT methods in combination with def2-TZVPP is tested against the newly generated CCSD(T)/CBS reference values. Detailed analysis of the results reveals that functional performance is dependent on the identity of the metal (with increased errors as group I is descended) and nucleic acid binding site (with larger errors for select purine coordination sites). Over all complexes considered, the best methods include the mPW2-PLYP double-hybrid and ωB97M-V RSH functionals (≤1.6% MPE; <1.0 kcal/mol MUE). If more computationally efficient approaches are required, the TPSS and revTPSS local meta-GGA functionals are reasonable alternatives (≤2.0% MPE; <1.0 kcal/mol MUE). Inclusion of counterpoise corrections to account for basis set superposition error only marginally improves the computed binding energies, suggesting that these corrections can be neglected with little loss in accuracy when using larger models that are necessary for describing biosystems and biomaterials. Overall, the most accurate functionals identified in this study will permit future works geared towards uncovering the impact of group I metals on the environment and human biology, designing new ways to selectively sense harmful metals, engineering modern biomaterials, and developing improved computational methods to more broadly study group I metal-nucleic acid interactions.

In this study, polyethersulfone (PES) ultrafiltration (UF) membranes were modified with GO, Ag, ZnO, Ag-GO and ZnO-GO nanoparticles to improve carbamazepine removal and fouling prevention by making membrane surfaces more hydrophilic. The fabricated membranes were characterized for surface and cross-sectional morphology, surface roughness and zeta potential, as well as hydrophilicity, functional groups, surface tension parameters and water permeability Thereafter, the membranes were evaluated for their efficiency in removing MgSO4 and carbamazepine as well as antifouling properties. To understand the role of affinity interactions in rejection and fouling, membrane-solute adhesion energies (∆Gslm) were quantified based on the Lifshitz-van der Waals/acid-base method. Unlike previous studies, which have generalized fouling prevention to be due to improvements in hydrophilicity upon adding nanoparticles, this work further explored the role of surface tension components on rejection and fouling prevention. The addition of nanoparticles improved membrane hydrophilicity (77-62°), water permeability (11.9-17.7 Lm-2 h-1 bar-1), mechanical strength (3.46-4.11 N/mm2), carbamazepine rejection (30-85%) and fouling prevention (60-23% flux decline). Rejection and antifouling properties increased as ∆Gslm became more repulsive (i.e., less negative). Membrane modification reduced irreversible fouling, and the fouled membranes were cleaned by flushing with water. Fouling related more to membrane electron donor components (γ-), while the roles of electron acceptor (γ+) and Lifshitz-van der Waals components (γLW) were less important. This work provides more insights into the role of affinity interactions in rejection and fouling and how rejection and fouling mechanisms change with nanoparticle addition.

The amino group of 2-amino-5-(4-halophenyl)-1,3,4-chalcogenadiazole has been replaced with bromo/iodo substituents to obtain a library of four compositionally related compounds. These are 2-iodo-5-(4-iodophenyl)-1,3,4-thiadiazole, C8H4I2N2S, 2-bromo-5-(4-bromophenyl)-1,3,4-selenadiazole, C8H4Br2N2Se, 2-bromo-5-(4-iodophenyl)-1,3,4-selenadiazole, C8H4BrIN2Se, and 2-bromo-5-(4-iodophenyl)-1,3,4-thiadiazole, C8H4BrIN2S. All were isostructural and contained bifurcated Ch...N (Ch is chalcogen) and X...X (X is halogen) interactions forming a zigzag packing motif. The noncovalent Ch...N interaction between the chalcogen-bond donor and the best-acceptor N atom appeared preferentially instead of a possible halogen bond to the same N atom. Hirshfeld surface analysis and energy framework calculations showed that, collectively, a bifurcated chalcogen bond was stronger than halogen bonding and this is more structurally influential in this system.

Possible polymorphic forms of the chemotherapy drug, temozolomide were predicted from the ab initio and DFT methods. The lattice minimization via distributed multipole analysis was carried out for the hypothetical generated structures. A crystal with unit cell parameters close to the real one and of same space group was retrieved, with partly similar packing and interactions. The analysis of inter molecular interaction (through Hirshfeld surface) and electrostatic potential reveals the complementary sites in the molecule. The 26 predicted structures were analyzed with respect to two computed lattice energies and hydrogen-bond propensity. The lattice energy of the real crystal [EXP] packing ranked number 6 compared on the basis of DMACRYS software and number 3 on the basis of the total lattice energy issued from the Crystalexplorer17 software at the B3LYP/6-31G∗∗ level of theory. The molecule has two strong hydrogen bond donors and five strong acceptors. The predicted packings are stabilized by one or two strong N-H…O/N-H…N as well as weak C-H…O/C-H…N and H…π hydrogen bonds. While the real structure with Z\' = 1, EXP, forms only one strong H-bond (N-H…O=C), several of the predicted packings form two strong H-bonds. Two predicted crystal packings have unit cell parameters close to the real structure, one of them shares several common intermolecular interactions.

The D2 subunit dopamine receptor represents a key factor in modulating dopamine release. Moreover, the investigated radiopharmaceutical ligands used in positron emission tomography imaging techniques are known to bind D2 receptors, allowing for dopaminergic pathways quantification in the living human brain. Thus, the biophysical characterization of these radioligands is expected to provide additional insights into the interaction mechanisms between the vehicle molecules and their targets. Using molecular dynamics simulations and QM calculations, the present study aimed to investigate the potential positions in which the D2 dopamine receptor would most likely interact with the three distinctive synthetic 11C-labeled compounds (raclopride (3,5-dichloro-N-[[(2S)-1-ethylpyrrolidin-2-yl]methyl]-2-hydroxy-6-methoxybenzamide)-RACL, FLB457 (5-bromo-N-[[(2S)-1-ethylpyrrolidin-2-yl]methyl]-2,3-dimethoxybenzamide)-FLB457 and SCH23390 (R(+)-7-Chloro-8-hydroxy-3-methyl-1-phenyl-2,3,4,5-tetrahydro-1H-3-benzazepine)-SCH)), as well as to estimate the binding affinities of the ligand-receptor complexes. A docking study was performed prior to multiple 50 ns molecular dynamics productions for the ligands situated at the top and bottom interacting pockets of the receptor. The most prominent motions for the RACL ligand were described by the high fluctuations of the peripheral aliphatic -CH3 groups and by its C-Cl aromatic ring groups. In good agreement with the experimental data, the D2 dopamine receptor-RACL complex showed the highest interacting patterns for ligands docked at the receptor\'s top position.

UNASSIGNED: The rapid spread of new SARS-CoV-2 variants has highlighted the crucial role played in the infection by mutations occurring at the SARS-CoV-2 spike receptor binding domain (RBD) in the interactions with the human ACE2 receptor. In this context, it urgently needs to develop new rapid tools for quickly predicting the affinity of ACE2 for the SARS-CoV-2 spike RBD protein variants to be used with the ongoing SARS-CoV-2 genomic sequencing activities in the clinics, aiming to gain clues about the transmissibility and virulence of new variants, to prevent new outbreaks and to quickly estimate the severity of the disease in the context of the 3PM.
UNASSIGNED: In our study, we used a computational pipeline for calculating the interaction energies at the SARS-CoV-2 spike RBD/ACE2 protein-protein interface for a selected group of characterized infectious variants of concern/interest (VoC/VoI). By using our pipeline, we built 3D comparative models of the SARS-CoV-2 spike RBD/ACE2 protein complexes for the VoC B.1.1.7-United Kingdom (carrying the mutations of concern/interest N501Y, S494P, E484K at the RBD), P.1-Japan/Brazil (RBD mutations: K417T, E484K, N501Y), B.1.351-South Africa (RBD mutations: K417N, E484K, N501Y), B.1.427/B.1.429-California (RBD mutations: L452R), the B.1.141 (RBD mutations: N439K), and the recent B.1.617.1-India (RBD mutations: L452R; E484Q) and the B.1.620 (RBD mutations: S477N; E484K). Then, we used the obtained 3D comparative models of the SARS-CoV-2 spike RBD/ACE2 protein complexes for predicting the interaction energies at the protein-protein interface.
UNASSIGNED: Along SARS-CoV-2 mutation database screening and mutation localization analysis, it was ascertained that the most dangerous mutations at VoC/VoI spike proteins are located mainly at three regions of the SARS-CoV-2 spike \"boat-shaped\" receptor binding motif, on the RBD domain. Notably, the P.1 Japan/Brazil variant present three mutations, K417T, E484K, N501Y, located along the entire receptor binding motif, which apparently determines the highest interaction energy at the SARS-CoV-2 spike RBD/ACE2 protein-protein interface, among those calculated. Conversely, it was also observed that the replacement of a single acidic/hydrophilic residue with a basic residue (E484K or N439K) at the \"stern\" or \"bow\" regions, of the boat-shaped receptor binding motif on the RBD, appears to determine an interaction energy with ACE2 receptor higher than that observed with single mutations occurring at the \"hull\" region or with other multiple mutants. In addition, our pipeline allowed searching for ACE2 structurally related proteins, i.e., THOP1 and NLN, which deserve to be investigated for their possible involvement in interactions with the SARS-CoV-2 spike protein, in those tissues showing a low expression of ACE2, or as a novel receptor for future spike variants. A freely available web-tool for the in silico calculation of the interaction energy at the SARS-CoV-2 spike RBD/ACE2 protein-protein interface, starting from the sequences of the investigated spike and/or ACE2 variants, was made available for the scientific community at: https://www.mitoairm.it/covid19affinities.
UNASSIGNED: In the context of the PPPM/3PM, the employment of the described pipeline through the provided webservice, together with the ongoing SARS-CoV-2 genomic sequencing, would help to predict the transmissibility of new variants sequenced from future patients, depending on SARS-CoV-2 genomic sequencing activities and on the specific amino acid replacement and/or on its location on the SARS-CoV-2 spike RBD, to put in play all the possible counteractions for preventing the most deleterious scenarios of new outbreaks, taking into consideration that a greater transmissibility has not to be necessarily related to a more severe manifestation of the disease.
UNASSIGNED: The online version contains supplementary material available at 10.1007/s13167-021-00267-w.

In this article we use MP2 and CCSD(T) calculations for the A24 and S66 data sets to explore how midbond functions can be used to generate cost effective counterpoise corrected supramolecular interaction energies of noncovalent complexes. We use the A24 data set to show that the primary role of midbond functions is not to approach the complete basis set limit, but rather to ensure a balanced description of the molecules and the interaction region (unrelated to the basis set superposition error). The need for balance is a consequence of using atom centered basis sets. In the complete basis set limit, the error will disappear, but reaching the complete basis set limit is not feasible beyond small systems. For S66 we investigate the need for increasing the number of midbond centers. Results show that adding a second midbond center increases the accuracy, but the effect is secondary to changing the atom centered basis set. Further, by comparing calculations using the 3s3p2d1f1g midbond set with using aug-cc-pVDZ and aug-cc-pVTZ as midbond sets, we see that the requirements for the midbond set to be effective, is not just that it contains diffuse functions, but also that high angular momentum functions are included. By comparing two approaches for placing midbond centers we show that results are not particularly sensitive to placement as long as the placement is reasonable.