Genetic determinants of IgG antibody response to COVID-19 vaccination.

Human humoral immune responses to SARS-CoV-2 vaccines exhibit substantial inter-individual variability and have been linked to vaccine efficacy. To elucidate the underlying mechanism behind this variability, we conducted a genome-wide association study (GWAS) on the anti-spike IgG serostatus of UK Biobank participants who were previously uninfected by SARS-CoV-2 and had received either the first dose (n = 54,066) or the second dose (n = 46,232) of COVID-19 vaccines. Our analysis revealed significant genome-wide associations between the IgG antibody serostatus following the initial vaccine and human leukocyte antigen (HLA) class II alleles. Specifically, the HLA-DRB1∗13:02 allele (MAF = 4.0%, OR = 0.75, p = 2.34e-16) demonstrated the most statistically significant protective effect against IgG seronegativity. This protective effect was driven by an alteration from arginine (Arg) to glutamic acid (Glu) at position 71 on HLA-DRß1 (p = 1.88e-25), leading to a change in the electrostatic potential of pocket 4 of the peptide binding groove. Notably, the impact of HLA alleles on IgG responses was cell type specific, and we observed a shared genetic predisposition between IgG status and susceptibility/severity of COVID-19. These results were replicated within independent cohorts where IgG serostatus was assayed by two different antibody serology tests. Our findings provide insights into the biological mechanism underlying individual variation in responses to COVID-19 vaccines and highlight the need to consider the influence of constitutive genetics when designing vaccination strategies for optimizing protection and control of infectious disease across diverse populations.

A novel approach to study multi-domain motions in JAK1's activation mechanism based on energy landscape.

The family of Janus Kinases (JAKs) associated with the JAK-signal transducers and activators of transcription signaling pathway plays a vital role in the regulation of various cellular processes. The conformational change of JAKs is the fundamental steps for activation, affecting multiple intracellular signaling pathways. However, the transitional process from inactive to active kinase is still a mystery. This study is aimed at investigating the electrostatic properties and transitional states of JAK1 to a fully activation to a catalytically active enzyme. To achieve this goal, structures of the inhibited/activated full-length JAK1 were modelled and the energies of JAK1 with Tyrosine Kinase (TK) domain at different positions were calculated, and Dijkstra's method was applied to find the energetically smoothest path. Through a comparison of the energetically smoothest paths of kinase inactivating P733L and S703I mutations, an evaluation of the reasons why these mutations lead to negative or positive regulation of JAK1 are provided. Our energy analysis suggests that activation of JAK1 is thermodynamically spontaneous, with the inhibition resulting from an energy barrier at the initial steps of activation, specifically the release of the TK domain from the inhibited Four-point-one, Ezrin, Radixin, Moesin-PK cavity. Overall, this work provides insights into the potential pathway for TK translocation and the activation mechanism of JAK1.

Engineering of substrate specificity in a plant cell-wall modifying enzyme through alterations of carboxyl-terminal amino acid residues.

Structural determinants of substrate recognition remain inadequately defined in broad specific cell-wall modifying enzymes, termed xyloglucan xyloglucosyl transferases (XETs). Here, we investigate the Tropaeolum majus seed TmXET6.3 isoform, a member of the GH16_20 subfamily of the GH16 network. This enzyme recognises xyloglucan (XG)-derived donors and acceptors, and a wide spectrum of other chiefly saccharide substrates, although it lacks the activity with homogalacturonan (pectin) fragments. We focus on defining the functionality of carboxyl-terminal residues in TmXET6.3, which extend acceptor binding regions in the GH16_20 subfamily but are absent in the related GH16_21 subfamily. Site-directed mutagenesis using double to quintuple mutants in the carboxyl-terminal region - substitutions emulated on barley XETs recognising the XG/penta-galacturonide acceptor substrate pair - demonstrated that this activity could be gained in TmXET6.3. We demonstrate the roles of semi-conserved Arg238 and Lys237 residues, introducing a net positive charge in the carboxyl-terminal region (which complements a negative charge of the acidic penta-galacturonide) for the transfer of xyloglucan fragments. Experimental data, supported by molecular modelling of TmXET6.3 with the XG oligosaccharide donor and penta-galacturonide acceptor substrates, indicated that they could be accommodated in the active site. Our findings support the conclusion on the significance of positively charged residues at the carboxyl terminus of TmXET6.3 and suggest that a broad specificity could be engineered via modifications of an acceptor binding site. The definition of substrate specificity in XETs should prove invaluable for defining the structure, dynamics, and function of plant cell walls, and their metabolism; these data could be applicable in various biotechnologies.

Unraveling pnicogen bonding cooperativity: Insights from molecular electrostatic potential analysis.

A theoretical investigation on the cooperativity of a series of binary, ternary, and quaternary complexes interconnected by pnicogen bonds has been conducted using calculations at the M06-2X/aug-cc-pVTZ level of density functional theory. By measuring changes in the molecular electrostatic potential (MESP) at the nucleus of interacting atoms in all of the complexes, it is possible to quantify the substantial reorganization of the electron density triggered by the formation of pnicogen bonds. The positive change in MESP, indicating a loss of electron density from the donor molecule in a dimer, facilitates the acceptance of electron density from a third molecule, resulting in the formation of a ternary complex with a stronger pnicogen bond compared to the one present in the binary complex. Similarly, the acceptor molecule in a dimer with a negative change in MESP showed an enhanced tendency to donate electron density to an electron-deficient third molecule. The MESP analysis provided valuable insights into the donor/acceptor characteristics of pnicogen bonds within the quaternary complexes. The proposed MESP hypotheses are consistent with the positive cooperativity observed in the pnicogen-bonded clusters. To quantify the changes in MESP, both at the donor atom (ΔVdonor ) and the acceptor atom (ΔVacceptor ), for all pnicogen bonds in the cluster, the total change in MESP (ΔΔVn ) was measured as ΔΔVn = ∑(ΔVdonor )-∑(ΔVacceptor ). Remarkably, ΔΔVn exhibited a strong linear relationship with the sum of the bond energies of the pnicogen bonds in the cluster. This establishes the MESP analysis as a robust approach for understanding the strength and cooperative behavior of pnicogen-bonded clusters. Additionally, the MESP features provided clear evidence of pnicogen bond formation, further supporting the reliability of this approach.

Constructing Directional Electrostatic Potential Difference via Gradient Nitrogen Doping for Efficient Oxygen Reduction Reaction.

Nitrogen doping has been recognized as an important strategy to enhance the oxygen reduction reaction (ORR) activity of carbon-encapsulated transition metal catalysts (TM@C). However, previous reports on nitrogen doping have tended to result in a random distribution of nitrogen atoms, which leads to disordered electrostatic potential differences on the surface of carbon layers, limiting further control over the materials' electronic structure. Herein, a gradient nitrogen doping strategy to prepare nitrogen-deficient graphene and nitrogen-rich carbon nanotubes encapsulated cobalt nanoparticles catalysts (Co@CNTs@NG) is proposed. The unique gradient nitrogen doping leads to a gradual increase in the electrostatic potential of the carbon layer from the nitrogen-rich region to the nitrogen-deficient region, facilitating the directed electron transfer within these layers and ultimately optimizing the charge distribution of the material. Therefore, this strategy effectively regulates the density of state and work function of the material, further optimizing the adsorption of oxygen-containing intermediates and enhancing ORR activity. Theoretical and experimental results show that under controlled gradient nitrogen doping, Co@CNTs@NG exhibits significantly ORR performance (Eonset = 0.96 V, E1/2 = 0.86 V). At the same time, Co@CNTs@NG also displays excellent performance as a cathode material for Zn-air batteries, with peak power density of 132.65 mA cm-2 and open-circuit voltage (OCV) of 1.51 V. This work provides an effective gradient nitrogen doping strategy to optimize the ORR performance.

A Critical Review of Short Antimicrobial Peptides from Scorpion Venoms, Their Physicochemical Attributes, and Potential for the Development of New Drugs.

Scorpion venoms have proven to be excellent sources of antimicrobial agents. However, although many of them have been functionally characterized, they remain underutilized as pharmacological agents, despite their evident therapeutic potential. In this review, we discuss the physicochemical properties of short scorpion venom antimicrobial peptides (ssAMPs). Being generally short (13-25 aa) and amidated, their proven antimicrobial activity is generally explained by parameters such as their net charge, the hydrophobic moment, or the degree of helicity. However, for a complete understanding of their biological activities, also considering the properties of the target membranes is of great relevance. Here, with an extensive analysis of the physicochemical, structural, and thermodynamic parameters associated with these biomolecules, we propose a theoretical framework for the rational design of new antimicrobial drugs. Through a comparison of these physicochemical properties with the bioactivity of ssAMPs in pathogenic bacteria such as Staphylococcus aureus or Acinetobacter baumannii, it is evident that in addition to the net charge, the hydrophobic moment, electrostatic energy, or intrinsic flexibility are determining parameters to understand their performance. Although the correlation between these parameters is very complex, the consensus of our analysis suggests that there is a delicate balance between them and that modifying one affects the rest. Understanding the contribution of lipid composition to their bioactivities is also underestimated, which suggests that for each peptide, there is a physiological context to consider for the rational design of new drugs.

Electrostatic Surface Potentials and Chalcogen-Bonding Motifs of Substituted 2,1,3-Benzoselenadiazoles Probed via 77Se Solid-State NMR Spectroscopy.

Chalcogen bonds (ChB) are moderately strong, directional, and specific non-covalent interactions that have garnered substantial interest over the last decades. However, ChB applications are currently hampered by a lack of methods to characterize and control chalcogen bonds. We report on the influence of various substituents (halogens, cyano, and methyl groups) on the observed self-complementary ChB networks of 2,1,3-benzoselenadiazoles. From molecular electrostatic potential calculations, we show that the electrostatic surface potentials (ESP) of the σ-holes on selenium are largely influenced by the electron-withdrawing character of these substituents. Structural analyses via X-ray diffraction reveal a variety of ChB geometries and binding modes that are rationalized via the computed ESP maps, although the structure of 5,6-dimethyl-2,1,3-benzoselenadiazole also demonstrates the influence of steric interactions. 77Se solid-state magic-angle spinning NMR spectroscopy, in particular the analysis of the selenium chemical shift tensors, is found to be an effective probe able to characterize both structural and electrostatic features of these self-complementary ChB systems. We find a positive correlation between the value of the ESP maxima at the σ-holes and the experimentally measured 77Se isotropic chemical shift, while the skew of the chemical shift tensor is established as a metric which is reflective of the ChB binding motif.

Halogen Bonding to the π-Systems of Polycyclic Aromatics.

The propensity of the π-electron system lying above a polycyclic aromatic system to engage in a halogen bond is examined by DFT calculations. Prototype Lewis acid CF3I is placed above the planes of benzene, naphthalene, anthracene, phenanthrene, naphthacene, chrysene, triphenyl, pyrene, and coronene. The I atom positions itself some 3.3-3.4 Å above the polycyclic plane, and the associated interaction energy is about 4 kcal/mol. This quantity is a little smaller for benzene, but is roughly equal for the larger polycyclics. The energy only oscillates a little as the Lewis acid slides across the face of the polycyclic, preferring regions of higher π-electron density over minima of the electrostatic potential. The binding is dominated by dispersion which contributes half of the total interaction energy.

Radial Behavior of Electrostatic Potentials of Atoms and Ions Revisited: Isotropy and Anisotropy.

In this paper we revisit earlier work relating to monoatomic atoms and ions published by pioneers in the area of electrostatic potentials. We include plots of the radial distributions of the electrostatic potentials for spherically symmetric atoms and cations, and for singly, doubly and triply negative anions. For atoms with anisotropy in their densities and electrostatic potentials, such as the halonium cations, it is shown how the molecular surface approach for plotting electrostatic potentials complements that achieved by directional radial distributions.

Investigation on the Influence of Solvents Environment on the Optoelectronic Properties of the Fluorescent Probe 1,3,4-Oxadiazole Analogues: A Combined Theoretical and Experimental Study.

This study investigates the photophysical properties of a nitrobenzene-substituted 1,3,4-oxadiazole derivative (OX-NO) using both theoretical and experimental methods. The impact of the solvent on OX-NO absorption and fluorescence spectra, as well as its fluorescence quantum yield, was initially studied. A noticeable bathochromic shift in the Stokes shift indicated a π→ π* transition within the molecules. Solute-solvent interactions were analysed using Catalan parameters, distinguishing between specific and nonspecific interactions. Excited state dipole moments were derived from Lippert's, Bakshiev's, and Chamma Viallet's equations, showing increased polarity in the excited state compared to the ground state. Ground state dipole moments were determined via solvatochromic shift methods and ab initio techniques. Additionally, detailed analyses of bond length, angles, dihedral angles, Mulliken charge distribution, and HOMO-LUMO energy gap were conducted using the DFT-B3LYP-6-311G basis set in Gaussian-09 W. The energy band gap values obtained from theoretical calculations and experimental methods (cyclic voltammetry and UV-Visible spectroscopy) exhibited excellent agreement. Reactive sites such as electrophilic and nucleophilic regions were identified through total electron density, electrostatic maps, molecular electrostatic potential, and 3D plots using DFT computational analysis. Global descriptors were employed to characterize the compounds' chemical reactivity comprehensively. The observed photophysical attributes underscore the potential of these fluorophores in various applications like organic light-emitting diodes, solar cells, and chemosensors. This study contributes crucial insights into the optoelectronic properties of nitrobenzene-substituted 1,3,4-oxadiazole derivatives, paving the way for their future integration in advanced technological domains.

DFT investigation of coupling constant anomalies in substituted ß-lactams.

ß-lactams are a chemically diverse group of molecules with a wide range of biological activities. Having recently observed curious trends in 2JHH coupling values in studies on this structural class, we sought to obtain a more comprehensive understanding of these diagnostic NMR parameters, specifically interrogating 1JCH, 2JCH, and 2JHH, to differentiate 3- and 4-monosubstituted ß-lactams. Further investigation using computational chemistry methods was employed to explore the geometric and electronic origins for the observed and calculated differences between the two substitution patterns.

A comprehensive study of experimental and theoretical characterization and in silico toxicity analysis of new molecules.

In this study, for the first time in the literature, a 2-(3-methoxyphenylamino)-2-oxoethyl acrylate (3MPAEA) molecule was synthesized in two steps, and a 2-chloro-N-(3-methoxyphenyl)acetamide (m-acetamide) was obtained in the first step. Experimental results were obtained using FTIR, 1H, and 13C NMR spectroscopy methods for m-acetamide and 3MPAEA compounds created in the laboratory environment and compared with theoretical results. Band gap (BG) energy, chemical hardness, electronegativity, chemical potential, and electrophilicity index were calculated. With vibration spectroscopic analysis, atom-molecule vibrations of the theoretical and experimental peaks of the spectrum were observed. The locations of C and H atoms were determined by nuclear magnetic resonance spectroscopy. The green, blue, and red regions of the potential energy map (MEP) map were examined. Some observed that the energy thermal, heat capacity, and entropy graphs increased in direct proportion to increasing the temperature in Kelvin, which is known as thermochemistry. The changes in the rotation, translation, and vibration of the molecule as its temperature increased were examined. When the thermochemistry surface map was examined, some observed that the temperature was high in the middle binding site of the molecules. Covalent interactions were graphed using the non-covalent interactions (NCIs) calculation method. In silico toxicity studies were carried out for m-acetamide and 3MPAEA molecules: fathead minnow LC50 (96 h), Daphnia magna LC50 (48 h), Tetrahymena pyriformis IGC50 (48 h), oral rat LD50, water solubility, bioconcentration factor, developmental toxicity, mutation, normal boiling point, flash point, melting point, density, thermal conductivity, viscosity, vapor pressure, etc. parameters were investigated.

Structure-Activity Relationship Models to Predict Properties of the Dielectric Fluids for Transformer Insulation System.

Mineral oils and synthetic and natural esters are the predominant insulating liquids in electrical equipment. Structure-activity relationship models to predict the key properties of pure insulating liquids, including pulse breakdown strengths, AC breakdown voltages, dielectric constants, flash points, and kinematic viscosities, have been proposed for the first time. Dependence of the specific properties on the molecular structures has been illustrated quantitatively in terms of surface area, statistical total variance, and average deviation of positive and negative electrostatic potentials, as augmented by molecular weight, volume, and ovality. Moreover, the individual contribution of the functional groups to viscosity has been revealed by an additive approach. The predicted properties are in good agreement with the experimental data. The present theoretical work provides new insights on the development of novel dielectric fluids.

Novel rod-like [Cu(phen)2(OAc)]·PF6 complex for high-performance visible-light-driven photocatalytic degradation of hazardous organic dyes: DFT approach, Hirshfeld and fingerprint plot analysis.

A novel octahedral distorted coordination complex was formed from a copper transition metal with a bidentate ligand (1,10-Phenanthroline) and characterized by Ultraviolet-visible spectroscopy, Ultraviolet-visible diffuse reflectance spectroscopy, Fourier-transform infrared spectroscopy, Brunauer-Emmett-Teller, Field emission scanning electron microscopy, and Single-crystal X-ray diffraction. The Hirshfeld surface and fingerprint plot analyses were conducted to determine the interactions between atoms in the Cu(II) complex. DFT calculations showed that the central copper ion and its coordinated atoms have an octahedral geometry. The Molecular electrostatic potential (MEP) map indicated that the copper (II) complex is an electrophilic compound that can interact with negatively charged macromolecules. The HOMO-LUMO analysis demonstrated the π nature charge transfer from acetate to phenanthroline. The band gap of [Cu(phen)2(OAc)]·PF6 photocatalyst was estimated to be 2.88 eV, confirming that this complex is suitable for environmental remediation. The photocatalytic degradation of erythrosine, malachite green, methylene blue, and Eriochrome Black T as model organic pollutants using the prepared complex was investigated under visible light. The [Cu(phen)2(OAc)]·PF6 photocatalyst exhibited degradation 94.7, 90.1, 82.7, and 74.3 % of malachite green, methylene blue, erythrosine, and Eriochrome Black T, respectively, under visible illumination within 70 min. The results from the Langmuir-Hinshelwood kinetic analysis demonstrated that the Cu(II) complex has a higher efficiency for the degradation of cationic pollutants than the anionic ones. This was attributed to surface charge attraction between photocatalyst and cationic dyes promoting removal efficiency. The reusability test indicated that the photocatalyst could be utilized in seven consecutive photocatalytic degradation cycles with an insignificant decrease in efficiency.

An Innovative Vortex-Assisted Liquid-Liquid Microextraction Approach Using Deep Eutectic Solvent: Application for the Spectrofluorometric Determination of Rhodamine B in Water, Food and Cosmetic Samples.

A new green and highly sensitive method for the determination of rhodamine B (RhB) by deep eutectic solvent-based vortex-assisted liquid-liquid microextraction with fluorescence detection (DES-VALLME-FLD) was developed. The extraction efficiency of conventional solvents and different deep eutectic solvent (DES) systems composed of tetrabutylammonium bromide (TBAB) and an alcohol (hexanol, octanol, or decanol) in different ratios were compared. DFT calculations of intermolecular electrostatic and non-covalent interactions of the most stable RhB forms with DES and water explain the experimental DESs' extraction efficiency. Semiempirical PM7 computations were used to obtain Hansen solubility parameters, which supported the good solubility of the monocationic RhB form in selected DESs. The dependence of the linear calibration of microextraction into 100 µL DES was observed in the RhB calibration range from 0.2 to 10.0 µg L-1 with a correlation coefficient of R2 = 0.9991. The LOD value was calculated to be 0.023 µg L-1. The accuracy and precision of the proposed method were verified over two days with RSD values of 2.9 to 4.1% and recovery of 94.6 to 103.7%. The developed method was applied to the determination of RhB in real samples (tap water, energy drink, and lipstick).

Cocrystal Prediction of Nifedipine Based on the Graph Neural Network and Molecular Electrostatic Potential Surface.

Nifedipine (NIF) is a dihydropyridine calcium channel blocker primarily used to treat conditions such as hypertension and angina. However, its low solubility and low bioavailability limit its effectiveness in clinical practice. Here, we developed a cocrystal prediction model based on Graph Neural Networks (CocrystalGNN) for the screening of cocrystals with NIF. And scoring 50 coformers using CocrystalGNN. To validate the reliability of the model, we used another prediction method, Molecular Electrostatic Potential Surface (MEPS), to verify the prediction results. Subsequently, we performed a second validation using experiments. The results indicate that our model achieved high performance. Ultimately, cocrystals of NIF were successfully obtained and all cocrystals exhibited better solubility and dissolution characteristics compared to the parent drug. This study lays a solid foundation for combining virtual prediction with experimental screening to discover novel water-insoluble drug cocrystals.

Achieving High Fill Factor in Organic Photovoltaic Cells by Tuning Molecular Electrostatic Potential Fluctuation.

In the field of organic photovoltaics (OPVs), significant progress has been made in tailoring molecular structures to enhance the open-circuit voltage and the short-circuit current density. However, there remains a crucial gap in the development of coordinated material design strategies focused on improving the fill factor (FF). Here, we introduce a molecular design strategy that incorporates electrostatic potential fluctuation to design organic photovoltaic materials. By reducing the fluctuation amplitude of IT-4F, we synthesized a new acceptor named ITOC6-4F. When using PBQx-TF as a donor, the ITOC6-4F-based cell shows a markedly low recombination rate constant of 0.66×10-14 cm3 s-1 and demonstrates an outstanding FF of 0.816, both of which are new records for binary OPV cells. Also, we find that a small fluctuation amplitude could decrease the energetic disorder of OPV cells, reducing energy loss. Finally, the ITOC6-4F-based cell creates the highest efficiency of 16.0 % among medium-gap OPV cells. Our work holds a vital implication for guiding the design of high-performance OPV materials.

A density fitting scheme for the fast evaluation of molecular electrostatic potential.

Molecular electrostatic potential (MEP) is a significant and crucial physical quantity that can be applied to a large number of scenarios, such as the prediction of nucleophilic or electrophilic attacks, fitting atomic charges, σ-hole, and so forth. The computational cost for the MEP has an O(N2 ) scaling with the increase of atoms, which is intractable and laborious for macromolecules. Herein, a density fitting molecular electrostatic potential (DF-MEP) is used to reduce the computational costs for the macromolecular MEP. It is found that the accuracy of DF-MEP is almost identical to the conventional molecular electrostatic potential (Conv-MEP), while the computational costs can be reduced to an O(N) scaling, for example, the computational time of 699,200 grids for the Trp-cage molecule (304 atoms) only takes 16.6 s at the B3LYP-D3(BJ)/def2-SVP level of theory with 16 CPU cores compared with 3060.2 s for the Conv-MEP method.

On hidden anisotropy of formally charged fragments.

The proper and precise reproduction of the molecular electrostatic potential (MEP) is crucial to describe correctly electrostatic interactions in molecular modeling. Most of the classical molecular mechanics force fields for biomolecules and drug-like molecules use the atom-centered point charges to describe MEP. However, it has been systematically pointed out in literature that such an approximation is not always enough, and some groups, like amino group or heavy halogens, require the use of anisotropic model for better description of their MEP. At the same time, the formally charged groups have not been as extensively and systematically studied as their neutral counterparts. In this report, we demonstrate that the anisotropic models for formally charged groups do bring improvements in the reference MEP reproduction, that are comparable in magnitude to those for neutral groups.

Hydrogen bonds of OCNH motif in rings in drugs: A molecular electrostatic potential analysis.

The OCNH unit is one of the most frequently encountered structural motifs in rings in drugs which serves dual role as the proton donor through NH bond and proton acceptor through the CO bond. Here, we predicted the HB strength (Eint ) of OCNH motif with H2 O for commonly observed 37 rings in drugs with DFT method M06L/6-311++G(d,p). The HB strength is rationalized in terms of molecular electrostatic potential (MESP) topology parameters ΔVn(NH) and ΔVn(CO) which describe the relative electron deficient/rich nature of NH and CO, respectively, with respect to the reference formamide. The Eint of formamide is -10.0 kcal/mol whereas the Eint of ring systems is in the range -8.6 to -12.7 kcal/mol-a minor increase/decrease compared to the formamide. The variations in Eint are addressed using the MESP parameters ΔVn(NH) and ΔVn(CO) and proposed the hypothesis that a positive ΔVn(NH) enhances NH…Ow interaction while a negative ΔVn(CO) enhances the CO…Hw interaction. The hypothesis is proved by expressing Eint jointly as ΔVn(NH) and ΔVn(CO) and also verified for twenty FDA approved drugs. The predicted Eint for the drugs using ΔVn(NH) and ΔVn(CO) agreed well with the calculated Eint . The study confirms that even delicate variations in the electronic feature of a molecule can be quantified in terms of MESP parameters and they provide a priori prediction of the HB strength. The MESP topology analysis is recommended to understand the tunability of HB strength in drug motifs.