SARS-CoV-2 Spike Protein Interaction Space.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a +sense single-strand RNA virus. The virus has four major surface proteins: spike (S), envelope (E), membrane (M), and nucleocapsid (N), respectively. The constitutive proteins present a high grade of symmetry. Identifying a binding site is difficult. The virion is approximately 50-200 nm in diameter. Angiotensin-converting enzyme 2 (ACE2) acts as the cell receptor for the virus. SARS-CoV-2 has an increased affinity to human ACE2 compared with the original SAR strain. Topological space, and its symmetry, is a critical component in molecular interactions. By exploring this space, a suitable ligand space can be characterized accordingly. A spike protein (S) computational model in a complex with ACE 2 was generated using silica methods. Topological spaces were probed using high computational throughput screening techniques to identify and characterize the topological space of both SARS and SARS-CoV-2 spike protein and its ligand space. In order to identify the symmetry clusters, computational analysis techniques, together with statistical analysis, were utilized. The computations are based on crystallographic protein data bank PDB-based models of constitutive proteins. Cartesian coordinates of component atoms and some cluster maps were generated and analyzed. Dihedral angles were used in order to compute a topological receptor space. This computational study uses a multimodal representation of spike protein interactions with some fragment proteins. The chemical space of the receptors (a dimensional volume) suggests the relevance of the receptor as a drug target. The spike protein S of SARS and SARS-CoV-2 is analyzed and compared. The results suggest a mirror symmetry of SARS and SARS-CoV-2 spike proteins. The results show thatSARS-CoV-2 space is variable and has a distinct topology. In conclusion, surface proteins grant virion variability and symmetry in interactions with a potential complementary target (protein, antibody, ligand). The mirror symmetry of dihedral angle clusters determines a high specificity of the receptor space.

Atoms-In-Molecules' Faces of Chemical Hardness by Conceptual Density Functional Theory.

The chemical hardness concept and its realization within the conceptual density functional theory is approached with innovative perspectives, such as the electronegativity and hardness equalization of atoms in molecules connected with the softness kernel, in order to examine the structure-reactivity equalization ansatz between the electronic sharing index and the charge transfer either in the additive or geometrical mean picture of bonding. On the other hand, the maximum hardness principle presents a relation with the chemical stability of the hardness concept. In light of the inverse relation between hardness and polarizability, the minimum polarizability principle has been proposed. Additionally, this review includes important applications of the chemical hardness concept to solid-state chemistry. The mentioned applications support the validity of the electronic structure principles regarding chemical hardness and polarizability in solid-state chemistry.

Three-Body Excitations in Fock-Space Coupled-Cluster: Fourth Order Perturbation Correction to Electron Affinity and Its Relation to Bondonic Formalism.

In this paper, we present a formulation of highly correlated Fock-space multi-reference coupled-cluster (FSMRCC) methods, including approximate triples on top of the FSMRCC with singles and doubles, which correct the electron affinities by at least at third and up to the fourth order in perturbation. We discuss various partial fourth-order schemes, which are reliable and yet computationally more efficient than the full fourth-order triples scheme. The third-order scheme is called MRCCSD+T*(3). We present two approximate fourth-order schemes, MRCCSD+T*-a(4) and MRCCSD+T*(4). The results that are presented allow one to choose an appropriate fourth-order scheme, which is less expensive and right for the problem. All these schemes are based on the effective Hamiltonian scheme, and provide a direct calculation of the vertical electron affinities. We apply these schemes to a prototype Li2 molecule, using four different basis sets, as well as BeO and CH+. We have calculated the vertical electron affinities of Li2 at the geometry of the neutral Li2 molecule. We also present the vertical ionization potentials of the Li2 anion at the geometry of the anion ground state. We have also shown how to calculate adiabatic electron affinity, though in that case we lose the advantages of direct calculation. BeO has been examined in two basis sets. For CH+, four different basis sets have been used. We have presented the partial fourth-order schemes to the EA in all the basis sets. The results are analyzed to illustrate the importance of triples, as well as highlight computationally efficient partial fourth-order schemes. The choice of the basis set on the electron affinity calculation is also emphasized. Comparisons with available experimental and theoretical results are presented. The general fourth-order schemes, which are conceptually equivalent with the Fock-space multi-reference coupled-cluster singles, doubles, and triplets (MRCCSD+T) methods, based on bondonic formalism, are also presented here in a composed way, for quantum electronic affinity.

The Density Functional Theory Account of Interplaying Long-Range Exchange and Dispersion Effects in Supramolecular Assemblies of Aromatic Hydrocarbons with Spin.

Aromatic hydrocarbons with fused benzene rings and regular triangular shapes, called n-triangulenes according to the number of rings on one edge, form groundstates with n-1 unpaired spins because of topological reasons. Here, we focus on methodological aspects emerging from the density functional theory (DFT) treatments of dimer models of the n = 2 triangulene (called also phenalenyl), observing that it poses interesting new problems to the issue of long-range corrections. Namely, the interaction comprises simultaneous spincoupling and van der Waals effects, i.e., a technical conjuncture not considered explicitly in the benchmarks calibrating long-range corrections for the DFT account of supramolecular systems. The academic side of considering dimer models for calculations and related analysis is well mirrored in experimental aspects, and synthetic literature revealed many compounds consisting of stacked phenalenyl cores, with intriguing properties, assignable to their long-range spin coupling. Thus, one may speculate that a thorough study assessing the performance of state-of-the-art DFT procedures has relevance for potential applications in spintronics based on organic compounds.

Chemical Bonding by the Chemical Orthogonal Space of Reactivity.

The fashionable Parr-Pearson (PP) atoms-in-molecule/bonding (AIM/AIB) approach for determining the exchanged charge necessary for acquiring an equalized electronegativity within a chemical bond is refined and generalized here by introducing the concepts of chemical power within the chemical orthogonal space (COS) in terms of electronegativity and chemical hardness. Electronegativity and chemical hardness are conceptually orthogonal, since there are opposite tendencies in bonding, i.e., reactivity vs. stability or the HOMO-LUMO middy level vs. the HOMO-LUMO interval (gap). Thus, atoms-in-molecule/bond electronegativity and chemical hardness are provided for in orthogonal space (COS), along with a generalized analytical expression of the exchanged electrons in bonding. Moreover, the present formalism surpasses the earlier Parr-Pearson limitation to the context of hetero-bonding molecules so as to also include the important case of covalent homo-bonding. The connections of the present COS analysis with PP formalism is analytically revealed, while a numerical illustration regarding the patterning and fragmentation of chemical benchmarking bondings is also presented and fundamental open questions are critically discussed.

Fock-Space Coupled Cluster Theory: Systematic Study of Partial Fourth Order Triples Schemes for Ionization Potential and Comparison with Bondonic Formalism.

In this paper, we have made a systematic study of partial fourth order perturbative schemes due to triples to compute the ionization potential within Fock-space multi-reference coupled-cluster theory. In particular, we have obtained computationally less expensive correlation schemes due to fourth order triples. Prototype examples have been considered to explore the efficacy of the approximate methods mentioned, while the bondonic formalism supporting the bonding phenomenology is also respectively for the first time here advanced.

QSAR by Minimal Topological Difference[s]: Post-Modern Perspectives.

In the context of reconsidering the Quantitative Structure-Activity Relationship (QSAR) methods at the economical level, namely the optimization rules of OECD, the present review unfolds the key features of Minimal Sterical, Monte-Carlo and Minimal Topological Difference (MTD) methods, developed for quantitative treatment of the relations between biological activity of organic chemical compounds (drugs, pesticides, and so on) and their structures. The initial Minimal Steric Difference (MSD) is completed by the three-dimensional variant of the MTD method, being the last one referred to here, while the main principles of validating and guiding a viable QSAR method verified by the analytical-automated MTD, thus enlarging the perspectives of understanding the chemical-biological interaction at the level of ligand-receptor sites, cavity, and walls, with a true service to the future adaptive molecular design.

Neurotoxicity of Pesticides: The Roadmap for the Cubic Mode of Action.

Pesticides are used today on a planetary-wide scale. The rising need for substances with this biological activity due to an increasing consumption of agricultural and animal products and to the development of urban areas makes the chemical industry to constantly investigate new molecules or to improve the physicochemical characteristics, increase the biological activities and improve the toxicity profiles of the already known ones. Molecular databases are increasingly accessible for in vitro and in vivo bioavailability studies. In this context, structure-activity studies, by their in silico - in cerebro methods, are used to precede in vitro and in vivo studies in plants and experimental animals because they can indicate trends by statistical methods or biological activity models expressed as mathematical equations or graphical correlations, so a direction of study can be developed or another can be abandoned, saving financial resources, time and laboratory animals. Following this line of research the present paper reviews the Structure-Activity Relationship (SAR) studies and proposes a correlation between a topological connectivity index and the biological activity or toxicity made as a result of a study performed on 11 molecules of organophosphate compounds, randomly chosen, with a basic structure including a Phosphorus atom double bounded to an Oxygen atom or to a Sulfur one and having three other simple covalent bonds with two alkoxy (-methoxy or -ethoxy) groups and to another functional group different from the alkoxy groups. The molecules were packed on a cubic structure consisting of three adjacent cubes, respecting a principle of topological efficiency, that of occupying a minimal space in that cubic structure, a method that was called the Clef Method. The central topological index selected for correlation was the Wiener index, since it was possible this way to discuss different adjacencies between the nodes in the graphs corresponding to the organophosphate compounds molecules packed on the cubic structure; accordingly, "three dimensional" variants of these connectivity indices could be considered and further used for studying the qualitative-quantitative relationships for the specific molecule-enzyme interaction complexes, including correlation between the Wiener weights (nodal specific contributions to the total Wiener index of the molecular graph) and the biochemical reactivity of some of the atoms. Finally, when passing from SAR to Q(uantitative)-SAR studies, especially by the present advanced method of the cubic molecule (Clef Method) and its good assessment of the (neuro)toxicity of the studied molecules and of their inhibitory effect on the target enzyme - acetylcholinesterase, it can be seen that a predictability of the toxicity and activity of different analogue compounds can be ensured, facilitating the in vivo experiments or improving the usage of pesticides.

Atoms in Generalized Orbital Configurations: Towards Atom-Dedicated Density Functionals.

Deriving a practical formula for the atomic body with generalized shell occupations, we perform a detective analysis of the radial distribution in the exchange energy, hinting at ideas about new types of density functionals, dedicated to the specifics of the electronic structure of atoms, exploiting the intrinsic spherical symmetry.

Ligand Shaping in Induced Fit Docking of MraY Inhibitors. Polynomial Discriminant and Laplacian Operator as Biological Activity Descriptors.

Docking-i.e., interaction of a small molecule (ligand) with a proteic structure (receptor)-represents the ground of drug action mechanism of the vast majority of bioactive chemicals. Ligand and receptor accommodate their geometry and energy, within this interaction, in the benefit of receptor-ligand complex. In an induced fit docking, the structure of ligand is most susceptible to changes in topology and energy, comparative to the receptor. These changes can be described by manifold hypersurfaces, in terms of polynomial discriminant and Laplacian operator. Such topological surfaces were represented for each MraY (phospho-MurNAc-pentapeptide translocase) inhibitor, studied before and after docking with MraY. Binding affinities of all ligands were calculated by this procedure. For each ligand, Laplacian and polynomial discriminant were correlated with the ligand minimum inhibitory concentration (MIC) retrieved from literature. It was observed that MIC is correlated with Laplacian and polynomial discriminant.

Molecular Dynamic Studies of the Complex Polyethylenimine and Glucose Oxidase.

Glucose oxidase (GOx) is an enzyme produced by Aspergillus, Penicillium and other fungi species. It catalyzes the oxidation of ß-d-glucose (by the molecular oxygen or other molecules, like quinones, in a higher oxidation state) to form d-glucono-1,5-lactone, which hydrolyses spontaneously to produce gluconic acid. A coproduct of this enzymatic reaction is hydrogen peroxide (H2O2). GOx has found several commercial applications in chemical and pharmaceutical industries including novel biosensors that use the immobilized enzyme on different nanomaterials and/or polymers such as polyethylenimine (PEI). The problem of GOx immobilization on PEI is retaining the enzyme native activity despite its immobilization onto the polymer surface. Therefore, the molecular dynamic (MD) study of the PEI ligand (C14N8_07_B22) and the GOx enzyme (3QVR) was performed to examine the final complex PEI-GOx stabilization and the affinity of the PEI ligand to the docking sites of the GOx enzyme. The docking procedure showed two places/regions of major interaction of the protein with the polymer PEI: (LIG1) of -5.8 kcal/mol and (LIG2) of -4.5 kcal/mol located inside the enzyme and on its surface, respectively. The values of enthalpy for the PEI-enzyme complex, located inside of the protein (LIG1) and on its surface (LIG2) were computed. Docking also discovered domains of the GOx protein that exhibit no interactions with the ligand or have even repulsive characteristics. The structural data clearly indicate some differences in the ligand PEI behavior bound at the two places/regions of glucose oxidase.

Learning the Relationship between the Primary Structure of HIV Envelope Glycoproteins and Neutralization Activity of Particular Antibodies by Using Artificial Neural Networks.

The dependency between the primary structure of HIV envelope glycoproteins (ENV) and the neutralization data for given antibodies is very complicated and depends on a large number of factors, such as the binding affinity of a given antibody for a given ENV protein, and the intrinsic infection kinetics of the viral strain. This paper presents a first approach to learning these dependencies using an artificial feedforward neural network which is trained to learn from experimental data. The results presented here demonstrate that the trained neural network is able to generalize on new viral strains and to predict reliable values of neutralizing activities of given antibodies against HIV-1.

Hydractinia echinata test-system. IV. Toxic synergism of human pharmaceuticals in mixtures with iodoform.

The structure-toxicity relationships for a series of singular human stomatological pharmaceuticals preparations and in mixture with Iodoform on Hydractinia echinata were obtained and their synergism was analyzed through the Metamorphosis Reduction Concentration (MRC50) within the "Köln model". The differences manifested between the total and individual components of the samples and mixtures, associated with toxic versus non-toxic synergism, are dependent on three essential factors of synthesis (the nature, the concentration and the number) besides manifested isotoxicity of the given components. The method represents a practical alternative useful for the reduction of experimental tests on animal to the lowest possible level, in accordance to the '3Rs' (reduce, reuse and recycle) integrative concept.

Chemical Structure-Biological Activity Models for Pharmacophores' 3D-Interactions.

Within medicinal chemistry nowadays, the so-called pharmaco-dynamics seeks for qualitative (for understanding) and quantitative (for predicting) mechanisms/models by which given chemical structure or series of congeners actively act on biological sites either by focused interaction/therapy or by diffuse/hazardous influence. To this aim, the present review exposes three of the fertile directions in approaching the biological activity by chemical structural causes: the special computing trace of the algebraic structure-activity relationship (SPECTRAL-SAR) offering the full analytical counterpart for multi-variate computational regression, the minimal topological difference (MTD) as the revived precursor for comparative molecular field analyses (CoMFA) and comparative molecular similarity indices analysis (CoMSIA); all of these methods and algorithms were presented, discussed and exemplified on relevant chemical medicinal systems as proton pump inhibitors belonging to the 4-indolyl,2-guanidinothiazole class of derivatives blocking the acid secretion from parietal cells in the stomach, the 1-[(2-hydroxyethoxy)-methyl]-6-(phenylthio)thymine congeners' (HEPT ligands) antiviral activity against Human Immunodeficiency Virus of first type (HIV-1) and new pharmacophores in treating severe genetic disorders (like depression and psychosis), respectively, all involving 3D pharmacophore interactions.

Linear and Branched PEIs (Polyethylenimines) and Their Property Space.

A chemical property space defines the adaptability of a molecule to changing conditions and its interaction with other molecular systems determining a pharmacological response. Within a congeneric molecular series (compounds with the same derivatization algorithm and thus the same brute formula) the chemical properties vary in a monotonic manner, i.e., congeneric compounds share the same chemical property space. The chemical property space is a key component in molecular design, where some building blocks are functionalized, i.e., derivatized, and eventually self-assembled in more complex systems, such as enzyme-ligand systems, of which (physico-chemical) properties/bioactivity may be predicted by QSPR/QSAR (quantitative structure-property/activity relationship) studies. The system structure is determined by the binding type (temporal/permanent; electrostatic/covalent) and is reflected in its local electronic (and/or magnetic) properties. Such nano-systems play the role of molecular devices, important in nano-medicine. In the present article, the behavior of polyethylenimine (PEI) macromolecules (linear LPEI and branched BPEI, respectively) with respect to the glucose oxidase enzyme GOx is described in terms of their (interacting) energy, geometry and topology, in an attempt to find the best shape and size of PEIs to be useful for a chosen (nanochemistry) purpose.

Quantitative structure-activity/ecotoxicity relationships (QSAR/QEcoSAR) of a series of phosphonates.

In this paper the structure-toxicity relationship studies were performed for a series of 60 phosphonates. The toxicity of the compounds was determined by two ways: by quantifying the measured toxicity values, Mlog(1/MRIC50) collected by literature, for rodents species; second by using EcoSAR software version 1.11, for calculating the toxicity for fish species, considered as dependent variables and they were related to structural features obtained by molecular and quantum mechanics calculations. The QSAR/QEcoSAR was validated by multiple linear regression (MLR), although the purpose of this work was not to validate the model proposed, but rather to test the influence of structural parameters of the proposed model QSAR/QEcoSAR. The obtained models showed that the toxicity of phosphonates was influenced by steric and molecular geometry which cause inhibition of cholinesterase activity.

Double Variational Binding--(SMILES) Conformational Analysis by Docking Mechanisms for Anti-HIV Pyrimidine Ligands.

Variational quantitative binding-conformational analysis for a series of anti-HIV pyrimidine-based ligands is advanced at the individual molecular level. This was achieved by employing ligand-receptor docking algorithms for each molecule in the 1,3-disubstituted uracil derivative series that was studied. Such computational algorithms were employed for analyzing both genuine molecular cases and their simplified molecular input line entry system (SMILES) transformations, which were created via the controlled breaking of chemical bonds, so as to generate the longest SMILES molecular chain (LoSMoC) and Branching SMILES (BraS) conformations. The study identified the most active anti-HIV molecules, and analyzed their special and relevant bonding fragments (chemical alerts), and the recorded energetic and geometric docking results (i.e., binding and affinity energies, and the surface area and volume of bonding, respectively). Clear computational evidence was also produced concerning the ligand-receptor pocket binding efficacies of the LoSMoc and BraS conformation types, thus confirming their earlier presence (as suggested by variational quantitative structure-activity relationship, variational-QSAR) as active intermediates for the molecule-to-cell transduction process.

Structure-biological function relationship extended to mitotic arrest-deficient 2-like protein Mad2 native and mutants-new opportunity for genetic disorder control.

Overexpression of mitotic arrest-deficient proteins Mad1 and Mad2, two components of spindle assembly checkpoint, is a risk factor for chromosomal instability (CIN) and a trigger of many genetic disorders. Mad2 transition from inactive open (O-Mad2) to active closed (C-Mad2) conformations or Mad2 binding to specific partners (cell-division cycle protein 20 (Cdc20) or Mad1) were targets of previous pharmacogenomics studies. Here, Mad2 binding to Cdc20 and the interconversion rate from open to closed Mad2 were predicted and the molecular features with a critical contribution to these processes were determined by extending the quantitative structure-activity relationship (QSAR) method to large-size proteins such as Mad2. QSAR models were built based on available published data on 23 Mad2 mutants inducing CIN-related functional changes. The most relevant descriptors identified for predicting Mad2 native and mutants action mechanism and their involvement in genetic disorders are the steric (van der Waals area and solvent accessible area and their subdivided) and energetic van der Waals energy descriptors. The reliability of our QSAR models is indicated by significant values of statistical coefficients: Cross-validated correlation q2 (0.53-0.65) and fitted correlation r2 (0.82-0.90). Moreover, based on established QSAR equations, we rationally design and analyze nine de novo Mad2 mutants as possible promoters of CIN.