The Impact of Interdisciplinary, Gender and Geographic Distributions on the Citation Patterns of the Journal of Chemical Information and Modeling.

There has been a growing recognition of the need for diversity and inclusion in scientific fields. This trend is reflected in the Journal of Chemical Information and Modeling (JCIM), where there has been a gradual increase in the number of papers that embrace this diversity. In this viewpoint, we analyze the evolution of the profile of papers published in JCIM from 1996 to 2022 addressing three diversity criteria, namely interdisciplinarity, geographic and gender distributions, and their impact on citation patterns. We used natural language processing tools for the classification of main areas and gender, as well as metadata, to analyze a total of 7384 articles published in the categories of research articles, reviews, and brief reports. Our analyses reveal that the relative number of articles and citation patterns are similar across the main areas within the scope of JCIM, and international collaboration and publications encompassing two to three research areas attract more citations. The percentage of female authors has increased from 1996 (less than 20%) to 2022 (more than 32%), indicating a positive trend toward gender diversity in almost all geographic regions, although the percentage of publications by single female authors remains lower than 20%. Most JCIM citations come from Europe and the Americas, with a tendency for JCIM papers to cite articles from the same continent. Furthermore, there is a correlation between the gender of the authors, as JCIM manuscripts authored by females are more likely to be cited by other JCIM manuscripts authored by females.

Bias Amplification in Gender, Gender Identity, and Geographical Affiliation.

In the quest for greater equity in science, individual attitudes and institutional policies should also embrace greater diversity and inclusion of minority groups. This viewpoint calls for a broader definition of gender bias in STEM to include gender identity and for increased attention to the issue of bias amplification due to geographic affiliation in the field of computational chemistry and chemoinformatics. It briefly discusses some active interventions to tackle bias on gender, gender identity, and geographic affiliation in STEM.

The Effect of Surface Composition on the Selective Capture of Atmospheric CO2 by ZIF Nanoparticles: The Case of ZIF-8.

We performed theoretical studies of CO2 capture in atmospheric conditions by the zeolitic imidazolate framework-8 (ZIF-8) via classical Monte Carlo (MC) simulations with Metropolis sampling and classical molecular dynamics (MD) simulations in the NVT and NPT ensembles and different thermodynamic conditions. The ZIF-8 framework was described by varying unit cell dimensions in the presence of pure gases of CO2, N2, O2, Ar, and H2O steam as well as binary mixtures of CO2:N2 and CO2:H2O in s 1:1 concentration. Different chemical compositions of the framework surface was considered to provide an accurate treatment of charge and charge distribution in the nanoparticle. Hence, surface groups were represented as unsaturated zinc atom (Zn+2), 2-methylimidazole (mImH), and deprotonated 2-methylimidazole (mIm-). Force field reparameterization of the surface sites was required to reproduce the interactions of the gas molecules with the ZIF-8 surface consistent with quantum mechanics (QM) calculations and Born-Oppenheimer molecular dynamics (BOMD). It was observed that ZIF-8 selectively captures CO2 due to the negligible concentrations of N2, O2, Ar, and H2O. These molecules spontaneously migrate to the inner pores of the framework. At the surface, there is a competitive interaction between H2O, CO2, and N2, for the positively charged ZIF-8 nanoparticle with a large binding energy advantage for water molecules (on average -62, -15, and -8 kcal/mol respectively). For the neutral ZIF-8 nanoparticle, the water molecules dominate the interactions due to the occurrence of hydrogen bond with the imidazolate groups at the surface. Simulations of binary mixtures of CO2/water steam and CO2/N2 were performed to investigate binding competition between these molecules for the framework positively charged and neutral surfaces. It was found that water molecules drastically block the interaction between CO2 molecules and the framework surface, decreasing CO2 capture in the central pore, and CO2 molecules fully block the interaction between N2 molecules and the framework. These findings show that CO2 capture by ZIF-8 is possible in atmospheric environments only upon dehydration of the atmospheric gas. It further shows that ZIF-8 capture of CO2 from the atmospheric environment is dependent on thermodynamic conditions and can be increased by decreasing temperature and/or increasing pressure.

Surface Assessment via Grid Evaluation (SuAVE) for Every Surface Curvature and Cavity Shape.

The surface assessment via grid evaluation (SuAVE) software was developed to account for the effect of curvature in the calculations of structural properties of chemical interfaces regardless of the chemical composition, asymmetry, and level of atom coarseness. It employs differential geometry techniques, enabling the representation of chemical surfaces as fully differentiable. In this article, we present novel developments of SuAVE to treat closed surfaces and complex cavity shapes. These developments expand the repertoire of curvature-dependent analyses already available in the previous version of SuAVE (e.g., area per lipid, density profiles, membrane thickness, deuterium-order parameters, volume per lipid, and surface curvature angle) to include new functionalities applicable to soft matter (e.g., sphericity, average radius, principal moment of inertia, and roundness) and crystalline porous materials (e.g., pore diameter, internal void volume, total area, and the total void volume of the unit cell structure). SuAVE can accurately handle chemical systems with high and low atom density as demonstrated for two distinct chemical systems: the lipid A vesicle and a set of selected metal-organic frameworks. The SuAVE software v2.0 is fully parallel and benefits from a compiler that supports OpenMP. SuAVE is freely available from https://github.com/SuAVE-Software/source and https://www.biomatsite.net/.

The tug of war between Al3+ and Na+ for order-disorder transitions in lipid-A membranes.

Cations play a critical role in the stability and morphology of lipid-A aggregates by neutralizing, hydrating and cross-linking these glycolipid molecules. Monophosphorylated lipid-A is the major immunostimulatory principle in commercially available adjuvants containing Al3+ such as adjuvant system 04 (AS04). The antagonist/agonist immunomodulatory properties of lipid-A are associated with chemical variations (e.g. the number of acyl chains and phosphate groups) and their aggregate arrangements (e.g. lamellar, nonlamellar or mixed). Therefore, the identification of the active form of lipid-A can provide valuable guidance in the development of vaccine adjuvants capable of boosting the immune system with decreased reactogenicity. Although the effect of mono and divalent cations on the structural polymorphism and endotoxicity of LPS has been previously investigated, much less is known about the effect of trivalent cations. We have investigated the effect of NaCl and AlCl3 salt solutions on the structural dynamics and stability of mono and diphosphorylated lipid-A membranes via atomistic MD simulations. The Al3+ ion exerts two major effects on the structural dynamics of lipid-A membranes. It acts as an efficient cross-linker of mono or diphosphorylated lipid-A molecules, thus stabilizing the lamellar arrangement of these glycolipids. It also alters the lipid-A packing and membrane fluidity, inducing disorder → order structural transitions of the membrane. This effect is promptly reversed upon the addition of NaCl solution, which promotes a nearly threefold increase in the amount of water in the carbohydrate moiety of the Al3+-containing lipid-A membranes. The exchange dynamics and residence times of cation-coordinated water molecules in these membranes provide insights into the molecular mechanism for the Na+-induced transition from a densely packed ordered phase to a disordered one. Al3+ counter-ions favor ordered lamellar aggregates, which has been previously associated with the lack of endotoxic activity and cytokine-inducing action. The resulting microscopic understanding of the structure and dynamics of lipid-A aggregates in the presence of Al3+ and Na+ salts can provide valuable guidance in the development of vaccine adjuvants capable of boosting the immune system with decreased reactogenicity.

SuAVE: A Tool for Analyzing Curvature-Dependent Properties in Chemical Interfaces.

Curvature is an intrinsic feature of biological membranes underlying vital cellular processes such as endocytosis, membrane fusion-fission, trafficking, and remodeling. The continuous expansion of the spatiotemporal scales accessible to computational simulations nowadays makes possible quasi-atomistic molecular dynamics simulations of these processes. In despite of that, computation of the shapes and curvatures associated with the dynamics of biological membranes remains challenging. For this reason, the effect of curvature is often neglected in the analysis of quantities essential for the accurate description of membrane properties (e.g., area and volume per lipid, density profiles, membrane thickness). We propose an algorithm for surface assessment via grid evaluation (SuAVE) that relies on the application of a radial base function to interpolate points scattered across an interface of any shape. This enables the representation of the chemical interface as fully differentiable so that related geometrical properties can be calculated through the straightforward employment of well-established differential geometry techniques. Hence, the effect of different types or degrees of curvature can be accurately taken into account in the calculations of structural properties of any interfaces regardless of chemical composition, asymmetry, and level of atom coarseness. The main functionalities implemented in SuAVE are featured for a number of tetraacylated and hexaacylated Lipid-A membranes of distinct curvatures and a surfactant micelle. We show that the properties calculated for moderately to highly curved membranes differ significantly between curvature-dependent and -independent algorithms. The SuAVE software is freely available from www.biomatsite.net/suave-software .

Rotational Profiler: A Fast, Automated, and Interactive Server to Derive Torsional Dihedral Potentials for Classical Molecular Simulations.

Rotational Profiler provides an analytical algorithm to compute sets of classical torsional dihedral parameters by fitting an empirical energy profile to a reference one that can be obtained experimentally or by quantum-mechanical methods. The resulting profiles are compatible with the functional forms in the most widely used biomolecular force fields (e.g., GROMOS, AMBER, OPLS, and CHARMM). The linear least-squares regression method is used to generate sets of parameters that best satisfy the fitting. Rotational Profiler is free to use, analytical, and force field/package independent. The formalism is herein described, and its usage, in an interactive and automated manner, is made available as a Web server at http://rotprof.lncc.br.

Discovery of 1,2,4-oxadiazole derivatives as a novel class of noncompetitive inhibitors of 3-hydroxykynurenine transaminase (HKT) from Aedes aegypti.

The mosquito Aedes aegypti is the vector of arboviruses such as Zika, Chikungunya, dengue and yellow fever. These infectious diseases have a major impact on public health. The unavailability of effective vaccines or drugs to prevent or treat most of these diseases makes vector control the main form of prevention. One strategy to promote mosquito population control is the use of synthetic insecticides to inhibit key enzymes in the metabolic pathway of these insects, particularly during larval stages. One of the main targets of the kynurenine detoxification pathway in mosquitoes is the enzyme 3-hydroxykynurenine transaminase (HKT), which catalyzes the conversion of 3-hydroxykynurenine (3-HK) into xanthurenic acid (XA). In this work, we report eleven newly synthesized oxadiazole derivatives and demonstrate that these compounds are potent noncompetitive inhibitors of HKT from Ae. aegypti. The present data provide direct evidence that HKT can be explored as a molecular target for the discovery of novel larvicides against Ae. aegypti. More importantly, it ensures that structural information derived from the HKT 3D-structure can be used to guide the development of more potent inhibitors.

Out of Sight, Out of Mind: The Effect of the Equilibration Protocol on the Structural Ensembles of Charged Glycolipid Bilayers.

Molecular dynamics (MD) simulations represent an essential tool in the toolbox of modern chemistry, enabling the prediction of experimental observables for a variety of chemical systems and processes and majorly impacting the study of biological membranes. However, the chemical diversity of complex lipids beyond phospholipids brings new challenges to well-established protocols used in MD simulations of soft matter and requires continuous assessment to ensure simulation reproducibility and minimize unphysical behavior. Lipopolysaccharides (LPS) are highly charged glycolipids whose aggregation in a lamellar arrangement requires the binding of numerous cations to oppositely charged groups deep inside the membrane. The delicate balance between the fully hydrated carbohydrate region and the smaller hydrophobic core makes LPS membranes very sensitive to the choice of equilibration protocol. In this work, we show that the protocol successfully used to equilibrate phospholipid bilayers when applied to complex lipopolysaccharide membranes occasionally leads to a small expansion of the simulation box very early in the equilibration phase. Although the use of a barostat algorithm controls the system dimension and particle distances according to the target pressure, fluctuation in the fleeting pressure occasionally enables a few water molecules to trickle into the hydrophobic region of the membrane, with spurious solvent buildup. We show that this effect stems from the initial steps of NPT equilibration, where initial pressure can be fairly high. This can be solved with the use of a stepwise-thermalization NVT/NPT protocol, as demonstrated for atomistic MD simulations of LPS/DPPE and lipid-A membranes in the presence of different salts using an extension of the GROMOS forcefield within the GROMACS software. This equilibration protocol should be standard procedure for the generation of consistent structural ensembles of charged glycolipids starting from atomic coordinates not previously pre-equilibrated. Although different ways to deal with this issue can be envisioned, we investigated one alternative that could be readily available in major MD engines with general users in mind.

Conformational Dynamics and Responsiveness of Weak and Strong Polyelectrolyte Brushes: Atomistic Simulations of Poly(dimethyl aminoethyl methacrylate) and Poly(2-(methacryloyloxy)ethyl trimethylammonium chloride).

The complex solution behavior of polymer brushes is key to control their properties, including for biomedical applications and catalysis. The swelling behavior of poly(dimethyl aminoethyl methacrylate) (PDMAEMA) and poly(2-(methacryloyloxy)ethyl trimethylammonium chloride) (PMETAC) in response to changes in pH, solvent, and salt types has been investigated using atomistic molecular dynamics simulations. PDMAEMA and PMETAC have been selected as canonical models for weak and strong polyelectrolytes whose complex conformational behavior is particularly challenging for the development and validation of atomistic models. The GROMOS-derived atomic parameters reproduce the experimental swelling coefficients obtained from ellipsometry measurements for brushes of 5-15 nm thickness. The present atomistic models capture the protonated morphology of PDMAEMA, the swollen and collapsed conformations of PDMAEMA and PMETAC in good and bad solvents, and the salt-selective response of PMETAC. The modular nature of the molecular models allows for the simple extension of atomic parameters to a variety of polymers or copolymers.

Polymyxin Binding to the Bacterial Outer Membrane Reveals Cation Displacement and Increasing Membrane Curvature in Susceptible but Not in Resistant Lipopolysaccharide Chemotypes.

Lipid-A is the causative agent of Gram-negative sepsis and is responsible for an increasingly high mortality rate among hospitalized patients. Compounds that bind Lipid-A can limit this inflammatory process. The cationic antimicrobial peptide polymyxin B (Pmx-B) is one of the simplest molecules capable of selectively binding to Lipid-A and may serve as a model for further development of Lipid-A binding agents. Gram-negative bacteria resistance to Pmx-B relies on the upregulation of a number of regulatory systems, which promote chemical modifications of the lipopolysaccharide (LPS) structure and leads to major changes in the physical-chemical properties of the outer membrane. A detailed understanding of how the chemical structure of the LPS modulates macroscopic properties of the outer membrane is paramount for the design and optimization of novel drugs targeting clinically relevant strains. We have performed a systematic investigation of Pmx-B binding to outer membrane models composed of distinct LPS chemotypes experimentally shown to be either resistant or susceptible to the peptide. Molecular dynamics simulations were carried out for Pmx-B bound to the penta- and hexa-acylated forms of Lipid-A (more susceptible) and Lipid-A modified with 4-amino-4-deoxy-l-arabinose (resistant) as well as the penta-acylated form of LPS Re (less susceptible). The present simulations show that upon binding to the bacterial outer membrane surface, Pmx-B promotes cation displacement and structural changes in membrane curvature and integrity as a function of the LPS chemotype susceptibility or resistance to the antimicrobial peptide.

Single-Amino Acid Modifications Reveal Additional Controls on the Proton Pathway of [FeFe]-Hydrogenase.

The proton pathway of [FeFe]-hydrogenase is essential for enzymatic H2 production and oxidation and is composed of four residues and a water molecule. A computational analysis of this pathway in the [FeFe]-hydrogenase from Clostridium pasteurianum revealed that the solvent-exposed residue of the pathway (Glu282) forms hydrogen bonds to two residues outside of the pathway (Arg286 and Ser320), implying that these residues could function in regulating proton transfer. In this study, we show that substituting Arg286 with leucine eliminates hydrogen bonding with Glu282 and results in an ∼3-fold enhancement of H2 production activity when methyl viologen is used as an electron donor, suggesting that Arg286 may help control the rate of proton delivery. In contrast, substitution of Ser320 with alanine reduces the rate ∼5-fold, implying that it either acts as a member of the pathway or influences Glu282 to permit proton transfer. Interestingly, quantum mechanics/molecular mechanics and molecular dynamics calculations indicate that Ser320 does not play a structural role or indirectly influence the barrier for proton movement at the entrance of the channel. Rather, it may act as an additional proton acceptor for the pathway or serve in a regulatory role. While further studies are needed to elucidate the role of Ser320, collectively these data provide insights into the complex proton transport process.