Atomic-defect engineering in thin membranes provides opportunities for ionic and molecular filtration and analysis. While molecular-dynamics (MD) calculations have been used to model conductance through atomic vacancies, corresponding experiments are lacking. We create sub-nanometer vacancies in suspended single-layer molybdenum disulfide (MoS) via Ga ion irradiation, producing membranes containing ∼300 to 1200 pores with average and maximum diameters of ∼0.5 and ∼1 nm, respectively. Vacancies exhibit missing Mo and S atoms, as shown by aberration-corrected scanning transmission electron microscopy (AC-STEM). The longitudinal acoustic band and defect-related photoluminescence were observed in Raman and photoluminescence spectroscopy, respectively. As the irradiation dose is increased, the median vacancy area remains roughly constant, while the number of vacancies (pores) increases. Ionic current versus voltage is nonlinear and conductance is comparable to that of ∼1 nm diameter single MoS pores, proving that the smaller pores in the distribution display negligible conductance. Consistently, MD simulations show that pores with diameters <0.6 nm are almost impermeable to ionic flow. Atomic pore structure and geometry, studied by AC-STEM, are critical in the sub-nanometer regime in which the pores are not circular and the diameter is not well-defined. This study lays the foundation for future experiments to probe transport in large distributions of angstrom-size pores.

Solid-state nanopores are single-molecule sensors that detect changes in ionic conductance (ΔG) when individual molecules pass through them. Producing high signal-to-noise ratio for the measurement of molecular structure in applications such as DNA sequencing requires low noise and large ΔG. The latter is achieved by reducing the nanopore diameter and membrane thickness. While the minimum diameter is limited by the molecule size, the membrane thickness is constrained by material properties. We use molecular dynamics simulations to determine the theoretical thickness limit of amorphous Si membranes to be ∼1 nm, and we designed an electron-irradiation-based thinning method to reach that limit and drill nanopores in the thinned regions. Double-stranded DNA translocations through these nanopores (down to 1.4 nm in thickness and 2.5 nm in diameter) provide the intrinsic ionic conductance detection limit in Si-based nanopores. In this regime, where the access resistance is comparable to the nanopore resistance, we observe the appearance of two conductance levels during molecule translocation. Considering the overall performance of Si-based nanopores, our work highlights their potential as a leading material for sequencing applications.

The performance of graphene oxide framework (GOF) membranes for water desalination is assessed using classical molecular dynamics (MD) simulations. The coupling between water permeability and salt rejection of GOF membranes is studied as a function of linker concentration n, thickness h and applied pressure ΔP. The simulations reveal that water permeability in GOF-(n,h) membranes can be tuned from ∼5 (n = 32 and h = 6.5 nm) to 400 L cm(-2) day(-1) MPa(-1) (n = 64 and h = 2.5 nm) and follows a Cnh(-αn) law. For a given pore size (n = 16 or 32), water permeability of GOF membranes increases when the pore spacing decreases, whereas for a given pore spacing (n = 32 or 64), water permeability increases by up to two orders of magnitude when the pore size increases. Furthermore, for linker concentrations n ≤ 32, the high water permeability corresponds to a 100% salt rejection, elevating this type of GOF membrane as an ideal candidate for water desalination. Compared to experimental performance of reverse osmosis membranes, our calculations suggest that under the same conditions of applied pressure and characteristics of membranes (ΔP ∼ 10 MPa and h ∼ 100 nm), one can expect a perfect salt rejection coupled to a water permeability two orders of magnitude higher than existing technologies, i.e., from a few cL cm(-2) day(-1) MPa(-1) to a few L cm(-2) day(-1) MPa(-1).

We report on the surface-catalyzed formal [2+2] and [2+2+2] cycloadditions of ortho-activated tetracene species on a Ag(111) substrate under ultrahigh vacuum conditions. Three different products are obtained: tetracene dimers, trimers, and tetramers. The former results from the formation of a four-membered ring while the other two arise from cyclization into six-membered rings. These on-surface reactions have been monitored by scanning tunneling microscopy and rationalized by density functional theory calculations. Our approach, based on the reaction of ortho-dihalo precursor monomers via formal cycloadditions, establishes an additional method for the highly active field of on-surface synthesis and enables the development of novel 1D and 2D covalent carbon nanostructures.

We report on the bottom-up fabrication of BN-substituted heteroaromatic networks achieved by surface-assisted polymerization and subsequent cyclodehydrogenation of specifically designed BN-substituted precursor monomers based on a borazine core structural element. To get insight into the cyclodehydrogenation pathway and the influence of molecular flexibility on network quality, two closely related precursor monomers with different degrees of internal cyclodehydrogenation have been employed. Scanning tunneling microscopy shows that, for both monomers, surface-assisted cyclodehydrogenation allows for complete monomer cyclization and the formation of covalently interlinked BN-substituted polyaromatic hydrocarbon networks on the Ag(111) surface. In agreement with experimental observations, density functional theory calculations reveal a significantly lower energy barrier for the cyclodehydrogenation of the conformationally more rigid precursor monomer, which is also reflected in a higher degree of long-range order of the obtained heteroaromatic network. Our proof-of-concept study will allow for the fabrication of atomically precise substitution patterns within BNC heterostructures.

Glutathione transferases (GSTs) are ubiquitous key enzymes that catalyse the conjugation of glutathione to xenobiotic compounds in the detoxification process. GSTs have been proposed to play a dual role in the signal termination of insect chemodetection by modifying odorant and tasting molecules and by protecting the chemosensory system. Among the 40 GSTs identified in Drosophila melanogaster, the Delta and Epsilon groups are insect-specific. GSTs Delta and Epsilon may have evolved to serve in detoxification, and have been associated with insecticide resistance. Here, we report the heterologous expression and purification of the D. melanogaster GST Delta 2 (GSTD2). We investigated the capacity of GSTD2 to bind tasting molecules. Among them, we found that isothiocyanates (ITC), insecticidal compounds naturally present in cruciferous plant and perceived as bitter, are good substrates for GSTD2. The X-ray structure of GSTD2 was solved, showing the absence of the classical Ser catalytic residue, conserved in the Delta and Epsilon GSTs. Using molecular dynamics, the interaction of ITC with the GSTD2 three-dimensional structure is analysed and discussed. These findings allow us to consider a biological role for GSTD2 in chemoperception, considering GSTD2 expression in the chemosensory organs and the potential consequences of insect exposure to ITC.

Computational studies of organic systems are frequently limited to static pictures that closely align with textbook style presentations of reaction mechanisms and isomerization processes. Of course, in reality chemical systems are dynamic entities where a multitude of molecular conformations exists on incredibly complex potential energy surfaces (PES). Here, we borrow a computational technique originally conceived to be used in the context of biological simulations, together with empirical force fields, and apply it to organic chemical problems. Replica‐exchange molecular dynamics (REMD) permits thorough exploration of the PES. We combined REMD with density functional tight binding (DFTB), thereby establishing the level of accuracy necessary to analyze small molecular systems. Through the study of four prototypical problems: isomer identification, reaction mechanisms, temperature‐dependent rotational processes, and catalysis, we reveal new insights and chemistry that likely would be missed using static electronic structure computations. The REMD‐DFTB methodology at the heart of this study is powered by i‐PI, which efficiently handles the interface between the DFTB and REMD codes. © 2015 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

ATP regulates the function of many proteins in the cell by transducing its binding and hydrolysis energies into protein conformational changes by mechanisms which are challenging to identify at the atomic scale. Based on molecular dynamics (MD) simulations, a method is proposed to analyze the structural changes induced by ATP binding to a protein by computing the effective free-energy landscape (FEL) of a subset of its coordinates along its amino-acid sequence. The method is applied to characterize the mechanism by which the binding of ATP to the nucleotide-binding domain (NBD) of Hsp70 propagates a signal to its substrate-binding domain (SBD). Unbiased MD simulations were performed for Hsp70-DnaK chaperone in nucleotide-free, ADP-bound and ATP-bound states. The simulations revealed that the SBD does not interact with the NBD for DnaK in its nucleotide-free and ADP-bound states whereas the docking of the SBD was found in the ATP-bound state. The docked state induced by ATP binding found in MD is an intermediate state between the initial nucleotide-free and final ATP-bound states of Hsp70. The analysis of the FEL projected along the amino-acid sequence permitted to identify a subset of 27 protein internal coordinates corresponding to a network of 91 key residues involved in the conformational change induced by ATP binding. Among the 91 residues, 26 are identified for the first time, whereas the others were shown relevant for the allosteric communication of Hsp70 s in several experiments and bioinformatics analysis. The FEL analysis revealed also the origin of the ATP-induced structural modifications of the SBD recently measured by Electron Paramagnetic Resonance. The pathway between the nucleotide-free and the intermediate state of DnaK was extracted by applying principal component analysis to the subset of internal coordinates describing the transition. The methodology proposed is general and could be applied to analyze allosteric communication in other proteins.