Geometric Complementarity Research Articles

A Combined Experimental and Computational Exploration of Heteroleptic cis-Pd2L2L'2 Coordination Cages through Geometric Complementarity.

Heteroleptic (mixed-ligand) coordination cages are of interest as host systems with more structurally and functionally complex cavities than homoleptic architectures. The design of heteroleptic cages, however, is far from trivial. In this work, we experimentally probed the self-assembly of Pd(II) ions with binary ligand combinations in a combinatorial fashion to search for new cis-Pd2L2L'2 heteroleptic cages. A hierarchy of computational analyses was then applied to these systems with the aim of elucidating key factors for rationalising self-assembly outcomes. Simple and inexpensive geometric analyses were shown to be effective in identifying complementary ligand pairs. Preliminary results demonstrated the viability of relatively rapid semi-empirical calculations for predicting the topology of thermodynamically favoured assemblies with rigid ligands, whilst more flexible systems proved challenging. Stemming from this, key challenges were identified for future work developing effective computational forecasting tools for self-assembled metallo-supramolecular systems.

In-cell bioluminescence resonance energy transfer (BRET)-based assay uncovers ceritinib and CA-074 as SARS-CoV-2 papain-like protease inhibitors

Papain-like protease (PLpro) is an attractive anti-coronavirus target. The development of PLpro inhibitors, however, is hampered by the limitations of the existing PLpro assay and the scarcity of validated active compounds. We developed a novel in-cell PLpro assay based on BRET and used it to evaluate and discover SARS-CoV-2 PLpro inhibitors. The developed assay demonstrated remarkable sensitivity for detecting the reduction of intracellular PLpro activity while presenting high reliability and performance for inhibitor evaluation and high-throughput screening. Using this assay, three protease inhibitors were identified as novel PLpro inhibitors that are structurally disparate from those previously known. Subsequent enzymatic assays and ligand-protein interaction analysis based on molecular docking revealed that ceritinib directly inhibited PLpro, showing high geometric complementarity with the substrate-binding pocket in PLpro, whereas CA-074 methyl ester underwent intracellular hydrolysis, exposing a free carboxyhydroxyl group essential for hydrogen bonding with G266 in the BL2 groove, resulting in PLpro inhibition.

Biaryl carboxamide-based peptidomimetics analogs as potential pancreatic lipase inhibitors for treating obesity.

A series of 1,1'-biphenyl-3-carboxamide and furan-phenyl-carboxamide analogs were synthesized using an optimized scheme and confirmed by 1H and 13C nuclear magnetic resonance and high-resolution mass spectrometry techniques. The synthesized peptidomimetics analogs were screened in vitro to understand the inhibitory potential of pancreatic lipase (PL). Analogs were assessed for the PL inhibitory activity based on interactions, geometric complementarity, and docking score. Among the synthesized analogs, 9, 29, and 24 were found to have the most potent PL inhibitory activity with IC50 values of 3.87, 4.95, and 5.34 µM, respectively, compared to that of the standard drug, that is, orlistat, which inhibits PL with an IC50 value of 0.99 µM. The most potent analog, 9, exhibited a competitive-type inhibition with an inhibition constant (Ki) of 2.72 µM. In silico molecular docking of analog 9 with the PL (PDB ID:1LPB) showed a docking score of -11.00 kcal/mol. Analog 9 formed crucial hydrogen bond interaction with Ser152, His263, π-cation interaction with Asp79, Arg256, and π-π stacking with Phe77, Tyr114 at the protein's active site. The molecular dynamic simulation confirmed that analog 9 forms stable interactions with PL at the end of 200 ns with root mean square deviation values of 2.5 and 6 Å. No toxicity was observed for analog 9 (concentration range of 1-20 µM) when tested by MTT assay in RAW 264.7 cells.

De novo design of peptides that bind specific conformers of α-synuclein.

Insoluble amyloids rich in cross-β fibrils are observed in a number of neurodegenerative diseases. Depending on the clinicopathology, the amyloids can adopt distinct supramolecular assemblies, termed conformational strains. However, rapid methods to study amyloids in a conformationally specific manner are lacking. We introduce a novel computational method for de novo design of peptides that tile the surface of α-synuclein fibrils in a conformationally specific manner. Our method begins by identifying surfaces that are unique to the conformational strain of interest, which becomes a "target backbone" for the design of a peptide binder. Next, we interrogate structures in the PDB with high geometric complementarity to the target. Then, we identify secondary structural motifs that interact with this target backbone in a favorable, highly occurring geometry. This method produces monomeric helical motifs with a favorable geometry for interaction with the strands of the underlying amyloid. Each motif is then symmetrically replicated to form a monolayer that tiles the amyloid surface. Finally, amino acid sequences of the peptide binders are computed to provide a sequence with high geometric and physicochemical complementarity to the target amyloid. This method was applied to a conformational strain of α-synuclein fibrils, resulting in a peptide with high specificity for the target relative to other amyloids formed by α-synuclein, tau, or Aβ40. This designed peptide also markedly slowed the formation of α-synuclein amyloids. Overall, this method offers a new tool for examining conformational strains of amyloid proteins.

An efficient computational protocol for template-based design of peptides that inhibit interactions involving SARS-CoV-2 proteins.

The RNA-dependent RNA polymerase (RdRp) complex of SARS-CoV-2 lies at the core of its replication and transcription processes. The interfaces between holo-RdRp subunits are highly conserved, facilitating the design of inhibitors with high affinity for the interaction interface hotspots. We, therefore, take this as a model protein complex for the application of a structural bioinformatics protocol to design peptides that inhibit RdRp complexation by preferential binding at the interface of its core subunit nonstructural protein, nsp12, with accessory factor nsp7. Here, the interaction hotspots of the nsp7-nsp12 subunit of RdRp, determined from a long molecular dynamics trajectory, are used as a template. A large library of peptide sequences constructed from multiple hotspot motifs of nsp12 is screened in-silico to determine sequences with high geometric complementarity and interaction specificity for the binding interface of nsp7 (target) in the complex. Two lead designed peptides are extensively characterized using orthogonal bioanalytical methods to determine their suitability for inhibition of RdRp complexation. Binding affinity of these peptides to accessory factor nsp7, determined using a surface plasmon resonance (SPR) assay, is slightly better than that of nsp12: dissociation constant of 133nM and 167nM, respectively, compared to 473nM for nsp12. A competitive ELISA is used to quantify inhibition of nsp7-nsp12 complexation, with one of the lead peptides giving an IC50 of 25μM . Cell penetrability and cytotoxicity are characterized using a cargo delivery assay and MTT cytotoxicity assay, respectively. Overall, this work presents a proof-of-concept of an approach for rational discovery of peptide inhibitors of SARS-CoV-2 protein-protein interactions.

Size of the protein-protein energy funnel in crowded environment.

Association of proteins to a significant extent is determined by their geometric complementarity. Large-scale recognition factors, which directly relate to the funnel-like intermolecular energy landscape, provide important insights into the basic rules of protein recognition. Previously, we showed that simple energy functions and coarse-grained models reveal major characteristics of the energy landscape. As new computational approaches increasingly address structural modeling of a whole cell at the molecular level, it becomes important to account for the crowded environment inside the cell. The crowded environment drastically changes protein recognition properties, and thus significantly alters the underlying energy landscape. In this study, we addressed the effect of crowding on the protein binding funnel, focusing on the size of the funnel. As crowders occupy the funnel volume, they make it less accessible to the ligands. Thus, the funnel size, which can be defined by ligand occupancy, is generally reduced with the increase of the crowders concentration. This study quantifies this reduction for different concentration of crowders and correlates this dependence with the structural details of the interacting proteins. The results provide a better understanding of the rules of protein association in the crowded environment.

Crescent-shaped meta-substituted benzene derivatives as a new class of non-nucleoside ribonuclease A inhibitors

Ribonuclease A is used as a model enzyme system for the design of RNase inhibitors. Previous studies have established that the geometric nature of the active site cleft is an important feature for the accommodation of crescent-shaped compounds in the active site of RNase A. In the current research, benzene-based triazolylated semicircular hybrid molecules carrying different polar functionalities were synthesized and screened for their RNase A inhibitory potency. An additional carboxylic acid group at the C1-position of the 1,3,5-trisubstituted benzene ring enhanced the inhibitory properties significantly. Furthermore, the studies revealed that the reduced arm lengths of 3,5-substituents result in a better geometric complementarity that makes the molecules fit favorably in the semicircular cavity of the active site as visualized by docking studies. In a series of ten such new compounds, the 3,5-bis[4-(sulfomethyl)-1H-1,2,3-triazol-1-yl]benzoic acid exhibited, the highest inhibition efficiency with a Ki value of 12 ± 0.9 µM. This study identifies a new class of non-nucleoside inhibitors which are competitive inhibitors of the ribonucleolytic activity of RNase A.

ESIPT-influenced C3-symmetry, disk-shaped fluorescence turn-on probe for Zn2+ based on melamine-naphthyl moiety with high affinity towards Cu2+ in CH3CN

A C3-symmetry turn-on probe for Zn2+ based on the melamine-naphthyl moiety (H) was synthesized via a simple one-step Schiff base reaction mechanism. The photophysical properties were investigated using spectroscopic methods of UV–Vis and fluorescence. The chemosensing properties of H were investigated via colorimetric methods, using heavy metal cations as well as biological anions, in aqueous soluble acetonitrile (CH3CN). Subsequently, the probe demonstrated that it could discriminate few cations (Zn2+, Co2+, Ni2+ and Cu2+) and anions (CN−, F−, AcO− and OH−) via colorimetric methods, observable by naked eye. The colorimetric method was further complemented by spectroscopic methods of UV–Vis, where spectral shifts were observed upon adding molar quantities of these cations and anions, evidencing that, chemical interactions has taken place, between the host and guest. Interestingly, fluorescence analysis carried out, demonstrated that the ESIPT-steered probe was highly selective to Zn2+ only, among the cations and anions used, turning the initially non-fluorescent H into a dual emissive H-Zn complex upon excited at 380 nm. The selectivity of H towards Zn2+ was based on the geometrical complementarity through coordination induced charge transfer, upon binding. The DFT predictions were carried out to complement the experimental data.

SARS-CoV-2 Variants, RBD Mutations, Binding Affinity, and Antibody Escape.

Since 2020, the receptor-binding domain (RBD) of the spike protein of the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been constantly mutating, producing most of the notable missense mutations in the context of “variants of concern”, probably in response to the vaccine-driven alteration of immune profiles of the human population. The Delta variant, in particular, has become the most prevalent variant of the epidemic, and it is spreading in countries with the highest vaccination rates, causing the world to face the risk of a new wave of the contagion. Understanding the physical mechanism responsible for the mutation-induced changes in the RBD’s binding affinity, its transmissibility, and its capacity to escape vaccine-induced immunity is the “urgent challenge” in the development of preventive measures, vaccines, and therapeutic antibodies against the coronavirus disease 2019 (COVID-19) pandemic. In this study, entropy–enthalpy compensation and the Gibbs free energy change were used to analyze the impact of the RBD mutations on the binding affinity of SARS-CoV-2 variants with the receptor angiotensin converting enzyme 2 (ACE2) and existing antibodies. Through the analysis, we found that the existing mutations have already covered almost all possible detrimental mutations that could result in an increase of transmissibility, and that a possible mutation in amino-acid position 498 of the RBD can potentially enhance its binding affinity. A new calculation method for the binding energies of protein–protein complexes is proposed based on the entropy–enthalpy compensation rule. All known structures of RBD–antibody complexes and the RBD–ACE2 complex comply with the entropy–enthalpy compensation rule in providing the driving force behind the spontaneous protein–protein docking. The variant-induced risk of breakthrough infections in vaccinated people is attributed to the L452R mutation’s reduction of the binding affinity of many antibodies. Mutations reversing the hydrophobic or hydrophilic performance of residues in the spike RBD potentially cause breakthrough infections of coronaviruses due to the changes in geometric complementarity in the entropy–enthalpy compensations between antibodies and the virus at the binding sites.

A Computational Study of Asymmetric Hydrogenation of 2‐Phenyl Acrylic Acids Catalyzed by a Rh(I) Catalyst with Ferrocenyl Chiral Bisphosphorus Ligand: The Role of Ion‐Pair Interaction†

Main observation and conclusionAsymmetric hydrogenation reaction is one of the most efficient synthetic methods to form useful chiral compounds for synthetic chemistry, medicinal chemistry and material chemistry. Generally, the enantioselectivity of many hydrogenation reactions is controlled by steric hindrance between the chiral ligand and substrate. Recently, Zhang group developed a highly asymmetric hydrogenation of 2‐aryl and 2‐alkyl acrylic acids catalyzed by a Rh(I) catalyst with a chiral Wudaphos ligand. The excellent enantioselectivity of this asymmetric reaction was proposed to be controlled by ion‐pair interaction between the substrate and chiral ligand. In this study, a systematic density functional theory study has been carried out to investigate the reaction mechanism and origin of the enantioselectivity. Our computational results suggest that this reaction follows the classic mechanism involving oxidative addition of H2, migratory insertion and reductive elimination. Different from the C=C coordination to the metal in the common oxidative addition step, our study found that the chelation of the carboxyl group of the substrate to the cationic Rh(I) metal is more favorable in this oxidative addition step. The high enantioselectivity is proposed to be dictated by a better catalyst/substrate geometric complementarity in the major pathway to have less distortion of the catalyst for a strong ion‐pair interaction.

Molecular and Structural Basis of Cross-Reactivity in M. tuberculosis Toxin-Antitoxin Systems.

Mycobacterium tuberculosis genome encodes over 80 toxin–antitoxin (TA) systems. While each toxin interacts with its cognate antitoxin, the abundance of TA systems presents an opportunity for potential non-cognate interactions. TA systems mediate manifold interactions to manage pathogenicity and stress response network of the cell and non-cognate interactions may play vital roles as well. To address if non-cognate and heterologous interactions are feasible and to understand the structural basis of their interactions, we have performed comprehensive computational analyses on the available 3D structures and generated structural models of paralogous M. tuberculosis VapBC and MazEF TA systems. For a majority of the TA systems, we show that non-cognate toxin–antitoxin interactions are structurally incompatible except for complexes like VapBC15 and VapBC11, which show similar interfaces and potential for cross-reactivity. For TA systems which have been experimentally shown earlier to disfavor non-cognate interactions, we demonstrate that they are structurally and stereo-chemically incompatible. For selected TA systems, our detailed structural analysis identifies specificity conferring residues. Thus, our work improves the current understanding of TA interfaces and generates a hypothesis based on congenial binding site, geometric complementarity, and chemical nature of interfaces. Overall, our work offers a structure-based explanation for non-cognate toxin-antitoxin interactions in M. tuberculosis.

Metallosupramolecules of pillar[5]-bis-trithiacrown including a mercury(ii) iodide ion-triplet complex.

A combination of pillar[5]-bis-trithiacrown (L) and mercury(ii) halides afforded a monomer complex (Cl--form), a 1-D coordination polymer (Br--form) and a supramolecular ion-triplet complex [(I·Hg·I)@L] (I--form). In the ion-triplet complex, the host encapsulates the (I--Hg2+-I-) entity via Hg2+π and C-HI- interactions, reflecting geometrical complementarity.

Hydrogen-Bond-Assisted Symmetry Breaking in a Network of Chiral Metal-Organic Assemblies.

Herein we elucidate the interplay of chiral, chelate, solvent, and hydrogen-bonding information in the self-assembly of a series of new three-dimensional metal-organic architectures. Enantiopure ligands, each containing H-bond donors and acceptors, form different structures, depending on the ratio in which they are combined: enantiopure components form M4L4 assemblies, whereas racemic mixtures form M3L3 stacks. Chiral amplification within M3L3 enantiomers was observed when a 2:1 ratio of R and S subcomponent enantiomers was employed. Simply switching the solvent (from MeCN to MeOH) or chelating unit (from bidentate to tridentate) increased the diversity of structures that can be generated from these building blocks, leading to the selective formation of novel M2L2 and M3L2 assemblies. The addition of achiral ligand building blocks resulted in the formation of further structures: When an achiral subcomponent was combined with its R and S chiral congeners, a three-layer heteroleptic architecture was generated, with the achiral unit sitting at the top of the stack. When combined with the S enantiomer only, however, the achiral unit assembled in the center of the structure, thus demonstrating the selective placement of achiral units within chiral systems. Further sorting experiments revealed that combining R and S stereocenters within a single ligand led to diastereoselective product generation. These results show how geometric complementarity between different ligands impacts upon the degree of hydrogen-bonding within the assembly, stabilizing specific low-symmetry architectures from among many possible structural outcomes.

Ring shape-dependent self-sorting of pillar[n]arenes assembled on a surface

Self-sorting, in which multiple components selectively assemble themselves by recognising self from others, is an attractive approach to produce supramolecular assemblies with controlled structures. Lock-and-key type complementary physical interactions are required for self-sorting because selective affinity is necessary to distinguish self from others. Here we show self-sorting behaviour based on a principle of geometrical complementarity by shape during our investigation of assembly of pentagonal pillar[5]arenes and hexagonal pillar[6]arenes on a surface. In the homoassembly systems, anionic pillar[5]arenes and pillar[6]arenes are adsorbed onto positively charged layers of cationic pillar[5]arenes and pillar[6]arenes, respectively, through cationic-anionic electrostatic interactions. In contrast, ionic pillar[5]arenes are adsorbed onto layers constructed from oppositely charged pillar[5]arenes, but ionic pillar[6]arenes are not. Equally, for the reverse combination, ionic pillar[6]arenes are adsorbed onto layers constructed from oppositely charged pillar[6]arenes, but ionic pillar[5]arenes are not. The geometrical complementarity by shape realises effective self-sorting even in non-directional multivalent ionic interactions.

HBX plays a potential role in developing hepatocellular carcinoma in chronic HBV carriers, a series of evidences suggesting plausible oncogenic factor

Background: Chronic Hepatitis caused by HBV is one of the principal causative agent for incidence of hepatocellular carcinoma around the world. In this respect, HBV viral protein HBX is possibly the most probable target that reinforces the prevalence of HBV mediated HCC. In this study HBX is used to delineate the molecular pathway it follows to develop oncogenesis. Methods & Materials: Around 500 complete genome sequences of HBV and respective HBX sequences were extracted from NCBI and UNIPROT database. They were subjected to multiple sequence alignment. Using maximum likelihood method the phylogenetic tree was reconstructed, engaging 1000 bootstrap replicates. Representative sequence from each genotype was selected to develop molecular model of HBX. By means of structural and thermodynamic refinements HBX dimers were constructed on basis of geometrical complementarity. Evidences suggests that longitudinal positioning of residues may lead to variations in the dimerization of HBX among different genotypes and subsequently their related functions. However, the dimer interaction region showed a region which could further be exploited for designing disruptor/small inhibitor. Moreover, selected sequences were predicted for post translational modifications analyses by different kinases to associate their molecular and biological functions. Results: HBX is intrinsically disordered protein structurally, providing justification for subsequent multiple molecular and biological roles including HBV associated carcinoma. HBX contains a Zinc Finger domain, which is an unambiguous evidence of protein-nucleic acid interaction. Structural variance is observed in molecular models inferring varying role in HBV associated hcc. Conformation of dimerization and variations in the post translational modification profile strengthens this fact further. PTM analysis exhibited 14 different kinases involved in the cell cycle regulation, regulate gene expression, signal transduction, apoptosis, cell metabolism, immune response and neuron differentiation tends to alter most of the sequences of HBX pointing towards their role in carcinogenesis. Conclusion: This study delineates the structural and functional polymorphism of HBX and provides enough evidence, suggesting HBX a plausible oncogene. HBX gene is found to be associated with development of hepatocellular carcinoma. This may also augment in understanding the molecular pathway followed by HBX in development of oncogenesis. Amplified with the further studies the present investigation may elucidate potential targets for therapeutic intervention.

An antiparallel double-stranded BODIPY-porphyrin dyad assembled via a self-complementary B-FZn interaction.

An antiparallel double-strand of a BODIPY-zinc-porphyrin dyad was assembled via geometrical complementarity of an unusual B-FZn coordination bonding interaction.

Binding of Mono- and Dianions within Silver Thiazolylurea Tweezers and Capsules.

A silver thiazolylurea complex, [Ag(TUTh)2]+, has been used as a host species for geometrically differently shaped mono- and dianions: trigonal planar (NO3-), tetrahedral (SO42-), and octahedral (SiF62-). In the presence of nitrate a 1:1 near-planar tweezer host-guest species is formed, with poor binding in solution despite excellent geometric complementarity being found between the host and the anion in the solid state. In the presence of either SO42- or SiF62- a 2:1 host-guest species is formed, whereby the guest is held in a capsulelike arrangement stabilized by an array of eight NH hydrogen-bond donors, as confirmed by X-ray crystallographic studies. Solution studies in DMSO-d6 support the host-guest stoichiometry seen in the solid state. The binding constant between SO42- and [Ag(TUTh)2]+ was calculated to be K21 = 2511 M-2 and was shown to be the dominant species in solution, in excellent agreement with the solid-state studies. However, upon the addition of SiF62- ions different speciation is observed: H2·G (capsule), H·G (tweezer), and H·G2 during the course of the study.

Geometric Complementarity in Assembly and Guest Recognition of a Bent Heteroleptic cis-[Pd2LA2LB2] Coordination Cage.

Due to the inherent difficulties in achieving a defined and exclusive formation of multicomponent assemblies against entropic predisposition, we present the rational assembly of a heteroleptic [Pd2LA2LB2]4+ coordination cage achieved through the geometric complementarity of two carefully designed ligands, LA and LB. With Pd(II) cations as rigid nodes, the pure distinctly angular components readily form homoleptic cages, a [Pd2LA4]4+ strained helical assembly and a [Pd4LB8]8+ box-like structure, both of which were characterized by X-ray analysis. Combined, however, the two ligands could be used to cleanly assemble a cis-[Pd2LA2LB2]4+ cage with a bent architecture. The same self-sorted product was also obtained by a quantitative cage-to-cage transformation upon mixing of the two homoleptic cages revealing the [Pd2LA2LB2]4+ assembly as the thermodynamic minimum. The structure of the heteroleptic cage was examined by ESI-MS, COSY, DOSY, and NOESY methods, the latter of which pointed toward a cis-conformation of ligands in the assembly. Indeed, DFT calculations revealed that the angular ligands and strict Pd(II) geometry strongly favor the cis-[Pd2LA2LB2]4+ species. The robust nature of the cis-[Pd2LA2LB2]4+ cage allowed us to probe the accessibility of its cavity, which could be utilized for shape recognition toward stereoisomeric guests. The ability to directly combine two different backbones in a controlled manner provides a powerful strategy for increasing complexity in the family of [Pd2L4] cages and opens up possibilities of introducing multiple functionalities into a single self-assembled architecture.

Structural analysis of Fc/FcγR complexes: a blueprint for antibody design.

The number of studies and the quality of the structural data of Fcγ receptors (FcγRs) has rapidly increased in the last few years. Upon critical examination of the literature, we have extracted general conclusions that could explain differences in affinity and selectivity of FcγRs for immunoglobulin G (IgG) based on structural considerations. FcγRs employ a little conserved asymmetric surface of domain D2 composed of two distinct subsites to recognize the well-conserved lower hinge region of IgG1-Fc. The extent of the contact interface with the antibody in subsite 1 of the receptor (but not in subsite 2), the geometrical complementarity between antibody and receptor, and the number of polar interactions contribute decisively toward strengthening the binding affinity of the antibody for the receptor. In addition, the uncertain role of the N-linked glycan of IgG for the binding and effector responses elicited by FcγRs is discussed. The available data suggest that not only the non-covalent interactions between IgG and FcγRs but also their dynamic features are essential for the immune response elicited through these receptors. We believe that the integration of structural, thermodynamic, and kinetic data will be critical for the design and validation of the next generation of therapeutic antibodies with enhanced effector capabilities.

Investigation of the complexation of essential oil components with cyclodextrins

The formation of inclusion complexes of six essential oil (EO) components (β-caryophyllene, cis-ocimene, trans-ocimene, sabinene hydrate (thujanol), γ-terpinene and α-terpineol) with six cyclodextrins (CDs) (α-CD, β-CD, γ-CD, HP-β-CD, RAMEB and CRYSMEB) was investigated by using static headspace-gas chromatography and UV–visible spectroscopy. Retention studies showed that CDs could efficiently reduce the volatility of EO components except for β-caryophyllene with α-CD. In this case, no inclusion complex was detected while for other compounds the formation of 1:1 inclusion complexes was observed. Results revealed that the inclusion stability mainly depends on geometric complementarity between encapsulated molecule and CD's cavity. Molecular modelling was used to investigate the complementarities between host and guest. Thus, CDs could efficiently be regarded as promising encapsulants for EO components leading to improve their application in cosmetic, pharmaceutical and agriculture fields.