Natures of benzene-water and pyrrole-water interactions in the forms of σ and π types: theoretical studies from clusters to liquid mixture

A Theoretical Study of Organotin Binding in Aromatase

The widely used organotin compounds are notorious for their acute toxicity. Experiments revealed that organotin might cause reproductive toxicity by reversibly inhibiting animal aromatase functioning. However, the inhibition mechanism is obscure, especially at the molecular level. Compared to experimental methods, theoretical approaches via computational simulations can help to gain a microscopic view of the mechanism. Here, in an initial attempt to uncover the mechanism, we combined molecular docking and classical molecular dynamics to investigate the binding between organotins and aromatase. The energetics analysis indicated that the van der Waals interaction is the primary driving force of binding the organic tail of organotin and the aromatase center. The hydrogen bond linkage trajectory analysis revealed that water plays a significant role in linking the ligand-water-protein triangle network. As an initial step in studying the mechanism of organotin inhibiting aromatase, this work provides an in-depth understanding of the binding mechanism of organotin. Further, our study will help to develop effective and environmentally friendly methods to treat animals that have already been contaminated by organotin, as well as sustainable solutions for organotin degradation.

Two-Sided Impact of Water on the Relaxation of Inner-Valence Vacancies of Biologically Relevant Molecules

After ionization of an inner-valence electron of molecules, the resulting cation-radicals store substantial internal energy which, if sufficient, can trigger ejection of an additional electron in an Auger decay usually followed by molecule fragmentation. In the environment, intermolecular Coulombic decay (ICD) and electron-transfer mediated decay (ETMD) are also operative, resulting in one or two electrons being ejected from a neighbor, thus preventing the fragmentation of the initially ionized molecule. These relaxation processes are investigated theoretically for prototypical heterocycle-water complexes of imidazole, pyrrole, and pyridine. It is found that the hydrogen-bonding site of the water molecule critically influences the nature and energetics of the electronic states involved, opening or closing certain relaxation processes of the inner-valence ionized system. Our results indicate that the relaxation mechanisms of biologically relevant systems with inner-valence vacancies on their carbon atoms can strongly depend on the presence of the electron-density donating or accepting neighbor, either water or another biomolecule.

Enhancement of Atmospheric Nucleation by Highly Oxygenated Organic Molecules: A Density Functional Theory Study

New particle formation (NPF) by gas-particle conversion is the main source of atmospheric aerosols. Highly oxygenated organic molecules (HOMs) and sulfuric acid (SA) are important NPF participants. 2-Methylglyceric acid (MGA), a kind of HOMs, is a tracer of isoprene-derived secondary organic aerosols. The nucleation mechanisms of MGA with SA were studied using density functional theory and atmospheric cluster dynamics simulation in this study, along with that of MGA with methanesulfonic acid (MSA) as a comparison. Our theoretical works indicate that the (MGA)(SA) and (MGA)(MSA) clusters are the most stable ones in the (MGA) _i(SA) _j ( i = 1-2, j = 1-2) and (MGA) _i(MSA) _j ( i = 1-2, j = 1-2) clusters, respectively. Both the formation rates of (MGA)(SA) and (MGA)(MSA) clusters are quite large and could have significant contributions to NPF. The results imply that the homomolecular nucleation of MGA is unlikely to occur in the atmosphere, and MGA and SA can effectively contribute to heteromolecular nucleation mainly in the form of heterodimers. MSA exhibits properties similar to SA in its ability to form clusters with MGA but is slightly weaker than SA.

Conformations of diethyl ether and its interaction with pyrrole at low temperatures

Conformations of diethyl ether (DEE) were studied at low temperatures in N₂ and Ar matrixes. Computations performed at B3LYP/aug-cc-pVDZ level of theory yielded three minima corresponding to tt, tg^± and g^±g^± conformers of DEE. Of the three, the tt and tg^± conformers of DEE were experimentally identified in N₂ and Ar matrixes. Furthermore, hydrogen bonded complexes of pyrrole (py) with DEE have been investigated using Density Functional Theory (DFT) and matrix isolation infrared spectroscopy. Computations performed at B3LYP level of theory using aug-cc-pVDZ basis set on pyrrole with tt and tg^± conformers of DEE gave py-DEE-tt and py-DEE-tg^± complexes, both characterized by NH⋯O interaction. Experimental evidence for the formation of py-DEE-tt and py-DEE-tg^± complexes was affirmed from the shifts in the NH stretching, NH bending regions of pyrrole and COC and CH stretching regions of DEE. NBO analysis was carried out to understand the charge-transfer delocalization interactions in the conformers of DEE and its hydrogen bonded complexes.

Dynamics of water in conical solvation shells around a benzene solute under different thermodynamic conditions

Water molecules in different parts of the anisotropic hydration shell of an aromatic molecule experience different interactions. In the present study, we investigate the anisotropic dynamics of water molecules in different non-overlapping conical shells around a benzene solute at sub- and supercritical conditions by means of molecular dynamics simulations using both non-polarizable and polarizable models. In addition to the dynamical properties, the effects of polarizability on the hydration structure of benzene at varying thermodynamic conditions are also investigated in the current study. The presence of πH-bonding in the solvation shell is found to be an important factor that influences the anisotropic dynamics of the benzene hydration shell. The water molecules located axial to the benzene plane are found to be maximally influenced by the πH-bonding. The extent of πH-bonding is found to be somewhat reduced on inclusion of polarizability. The πH-bonded water molecules are found to reorient through large-amplitude angular jumps where the jump-angle amplitude increases at higher temperatures and lower densities. For both non-polarizable and polarizable models, it is found that the water molecules in the axial conical shells possess faster orientational and hydrogen bond dynamics compared to those in the equatorial plane. Water molecules in the axial conical shells are also found to diffuse at a faster rate than bulk molecules due to the relatively weaker benzene-water πH-bonding interactions in the axial region of the hydration shell. The residence dynamics of water molecules in different conical solvation shells around the solute is also investigated in the current study.

Hydration of 3-hydroxy-4,4-dimethylglutaric acid with dimethylamine complex and its atmospheric implications

Changes in temperature affects the distribution of isomers, which facilitates the understanding of new particle formation in the atmosphere.

MD simulation study of the diffusion and local structure of n-alkanes in liquid and supercritical methanol at infinite dilution

The diffusion coefficients of 14 n-alkanes (ranging from methane to n-tetradecane) in liquid and supercritical methanol at infinite dilution (at a pressure of 10.5 MPa and at temperatures of 299 K and 515 K) were deduced via molecular dynamics simulations. Values for the radial distribution function, coordination number, and number of hydrogen bonds were then calculated to explore the local structure of each fluid. The flexibility of the n-alkane (as characterized by the computed dihedral distribution, end-to-end distance, and radius of gyration) was found to be a major influence and hydrogen bonding to be a minor influence on the local structure. Hydrogen bonding reduces the flexibility of the n-alkane, whereas increasing the temperature enhances its flexibility, with temperature having a greater effect than hydrogen bonding on flexibility. Graphical abstract The flexibility of the alkane is a major influence and the hydrogen bonding is a minor influence on the first solvation shell; the coordination numbers of long-chain n-alkanes in the first solvation shell are rather low.

Characterization of the nucleation precursor (H2SO4–(CH3)2NH) complex: intra-cluster interactions and atmospheric relevance

Amines have been proposed to participate in the nucleation process, but the electron density analysis and the determination of a temperature dependence of the clusters are still lacking.

Serine and Cysteine π-Interactions in Nature: A Comparison of the Frequency, Structure, and Stability of Contacts Involving Oxygen and Sulfur

Despite many DNA–protein π-interactions in high-resolution crystal structures, only four X–H···π or X···π interactions were found between serine (Ser) or cysteine (Cys) and DNA nucleobase π-systems in over 100 DNA–protein complexes (where X = O for Ser and X = S for Cys). Nevertheless, 126 non-covalent contacts occur between Ser or Cys and the aromatic amino acids in many binding arrangements within proteins. Furthermore, Ser and Cys protein–protein π-interactions occur with similar frequencies and strengths. Most importantly, due to the great stability that can be provided to biological macromolecules (up to –20 kJ mol–1 for neutral π-systems or –40 kJ mol–1 for cationic π-systems), Ser and Cys π-interactions should be considered when analyzing protein stability and function.

What are preferred water-aromatic interactions in proteins and crystal structures of small molecules?

The distribution of water molecules around aromatic rings in the protein structures and crystal structures of small molecules shows quite a small number of the strongest OH-π interactions, a larger number of parallel interactions, and the largest number of the weakest CH-O interactions.

Parallel water/aromatic interactions of non-coordinated and coordinated water

The parallel interactions of non-coordinated and coordinated water molecules with an aromatic ring were studied by analyzing data in the Cambridge structural database (CSD) and by using quantum chemical calculations. The CSD data show that water/aromatic contacts prefer parallel to OH/π interactions, which indicates the importance of parallel interactions. The results reveal the influence of water coordination to a metal ion; the interactions of aqua complexes are stronger. Coordinated water molecules prefer a parallel-down orientation in which one OH bond is parallel to the aromatic ring, whereas the other OH bond points to the plane of the ring. The interactions of aqua complexes with parallel-down water/benzene orientation are as strong as the much better known OH/π orientations. The strongest calculated interaction energy is -14.89 kcal mol(-1) . The large number of parallel contacts in crystal structures and the quite strong interactions indicate the importance of parallel orientation in water/benzene interactions.

Structural flexibility of 4,4'-methylene diphenyl diisocyanate (4,4'-MDI): evidence from first principles calculations

A reactant used globally in the production of polyurethane is the molecule 4,4′-methylene diphenyl diisocyanate (4,4′-MDI). The structural flexibility of 4,4′-MDI is one of the most important molecular properties influencing the polymerization process and this property was therefore modeled using density functional theory (DFT) calculations and Car-Parrinello molecular dynamics (MD) simulations. Global and local minima structures were found and confirmed by vibrational analysis. The energy barriers related to rotation of the aromatic rings were estimated by DFT calculations. The stability of global and local minima was verified by Car-Parrinello (MD) runs at finite temperature. The presence of weak C–H⋯π hydrogen bonds was confirmed by atoms in molecules analysis and found to be responsible for the low energy barriers.