Langevin Dipoles Research Articles

Fast and Effective Prediction of the Absolute Binding Free Energies of Covalent Inhibitors of SARS-CoV-2 Main Protease and 20S Proteasome.

The COVID-19 pandemic has been a public health emergency with continuously evolving deadly variants around the globe. Among many preventive and therapeutic strategies, the design of covalent inhibitors targeting the main protease (Mpro) of SARS-CoV-2 that causes COVID-19 has been one of the hotly pursued areas. Currently, about 30% of marketed drugs that target enzymes are covalent inhibitors. Such inhibitors have been shown in recent years to have many advantages that counteract past reservation of their potential off-target activities, which can be minimized by modulation of the electrophilic warhead and simultaneous optimization of nearby noncovalent interactions. This process can be greatly accelerated by exploration of binding affinities using computational models, which are not well-established yet due to the requirement of capturing the chemical nature of covalent bond formation. Here, we present a robust computational method for effective prediction of absolute binding free energies (ABFEs) of covalent inhibitors. This is done by integrating the protein dipoles Langevin dipoles method (in the PDLD/S-LRA/β version) with quantum mechanical calculations of the energetics of the reaction of the warhead and its amino acid target, in water. This approach evaluates the combined effects of the covalent and noncovalent contributions. The applicability of the method is illustrated by predicting the ABFEs of covalent inhibitors of SARS-CoV-2 Mpro and the 20S proteasome. Our results are found to be reliable in predicting ABFEs for cases where the warheads are significantly different. This computational protocol might be a powerful tool for designing effective covalent inhibitors especially for SARS-CoV-2 Mpro and for targeted protein degradation.

Revisiting the protomotive vectorial motion of F0-ATPase

The elucidation of the detailed mechanism used by F0 to convert proton gradient to torque and rotational motion presents a major puzzle despite significant biophysical and structural progress. Although the conceptual model has advanced our understanding of the working principles of such systems, it is crucial to explore the actual mechanism using structure-based models that actually reproduce a unidirectional proton-driven rotation. Our previous work used a coarse-grained (CG) model to simulate the action of F0 However, the simulations were based on a very tentative structural model of the interaction between subunit a and subunit c. Here, we again use a CG model but with a recent cryo-EM structure of cF1F0 and also explore the proton path using our water flooding and protein dipole Langevin dipole semimacroscopic formalism with its linear response approximation version (PDLD/S-LRA) approaches. The simulations are done in the combined space defined by the rotational coordinate and the proton transport coordinate. The study reproduced the effect of the protomotive force on the rotation of the F0 while establishing the electrostatic origin of this effect. Our landscape reproduces the correct unidirectionality of the synthetic direction of the F0 rotation and shows that it reflects the combined electrostatic coupling between the proton transport path and the c-ring conformational change. This work provides guidance for further studies in other proton-driven mechanochemical systems and should lead (when combined with studies of F1) to a complete energy transduction picture of the F0F1-ATPase system.

Solvation thermodynamics of neutral and charged solutes using the solvation‐layer interface condition continuum dielectric model

AbstractWe demonstrate that the solvation‐layer interface condition (SLIC) continuum dielectric model for molecular electrostatics, combined with a simple solvent‐accessible‐surface‐area (SASA)‐proportional model for nonpolar solvent effects, accurately predicts solvation entropies of neutral and charged small molecules. The SLIC/SASA model has only seven fitting parameters in total and achieves this accuracy using a training set with only 20 compounds. Despite this simplicity, solvation free energies and entropies are nearly as accurate as those predicted by the more sophisticated Langevin dipoles solvation model. Surprisingly, the model automatically reproduces the negligible contribution of electrostatics to the solvation of hydrophobic compounds. Opportunities for improvement include nonpolar solvation, anion solvation entropies, and heat capacities. More molecular realism may be needed for these quantities. To enable a future, explicit‐solvent‐based assessment of the SLIC/SASA implicit‐solvent model, we predict solvation entropies for the Mobley test set, which are available as Supporting Information.

Nucleic Bases Alkylation with Acrylonitrile and Cyanoethylene Oxide: A Computational Study.

Acrylonitrile (AN) is widely used in the manufacture of resins, plastics, and polymers, where workers are exposed to it during its production, transportation, and application. After intake a portion of AN is converted to cyanoethylene oxide (CEO) by cytochrome P450 2E1. Both AN and CEO represent possible chemical carcinogens leading to DNA damage mainly in the form of the major 7-(2-oxoethyl)deoxyguanosine adduct. A kinetic model for its formation was devised and a corresponding second-order rate constant obtained from the experimental data on the reaction with CEO. A series of ab initio, density functional theory, and semiempirical calculations of activation free energies was then performed on the alkylation of nucleic bases with both CEO and AN. The combination of Hartree-Fock level of theory with the flexible 6-311++G(d,p) basis set and Langevin dipoles implicit solvation model gave the best agreement with the experimental activation barrier. It also predicted relative reactivities of all four nucleobases that are in agreement with the experimentally reported adduct yields. Moreover, this combination predicted higher reactivity of CEO than AN with all four nucleobases corroborating the experimental hypothesis that SN2 substitution of CEO rather than direct Michael addition of AN is responsible for the genotoxic properties of AN. In a broader context this paper points to the applicability of quantum chemical methods to the studies of carcinogenesis.

Aggregation of photosensitizers: the role of dispersion and solvation on dimer formation energetics

A detailed account for the dimer formation energetics is computed for two commonly used photosensitizers: metal-free phthalocyanine (Pc) and methylene blue (MB), a cationic phenothiazine. Gas-phase calculations were performed using the density functional theory augmented with dispersion corrections using both the dispersion-correcting atom-centred potentials and the D2/D3 methods. SCS-MP2/CBS results were also provided for comparison. Hydration free energies were determined using continuum methodologies: namely the conductor-like polarizable continuum model, the solvation model density and a reparametrized version of the polarizable continuum model. For methylene blue, the Langevin dipole model was also used. Our results showed that the nature of aggregation is remarkably different for these two molecules, being dispersion driven for Pc while solvent driven for MB. For water and ethanol our results are in good agreement with the experimental data available for MB. We also provide quantitative estimates to the popular explanation put out by Mukerjee and Ghosh more than 40 years ago that aggregation arises from solvent entropic and inter-species enthalpic contributions.

Quantum mechanical analysis of nonenzymatic nucleotidyl transfer reactions: kinetic and thermodynamic effects of β-γ bridging groups of dNTP substrates.

Rate (k) and equilibrium (K) constants for the reaction of tetrahydrofuranol with a series of Mg2+ complexes of methyl triphosphate analogues, CH3O-P(O2)-O-P(O2)-X-PO34–, X = O, CH2, CHCH3, C(CH3)2, CFCH3, CHF, CHCl, CHBr, CFCl, CF2, CCl2, and CBr2, forming phosphate diester and pyrophosphate or bisphosphonate in aqueous solution were evaluated by B3LYP/TZVP//HF/6-31G* quantum chemical calculations and Langevin dipoles and polarized continuum solvation models. The calculated log k and log K values were found to depend linearly on the experimental pKa4 of the conjugate acid of the corresponding pyrophosphate or bisphosphonate leaving group. The calculated slopes of these Brønsted linear free energy relationships were βlg = −0.89 and βeq = −0.93, respectively. The studied compounds also followed the linear relationship Δlog k = 0.8Δlog K, which became less steep, Δlog k = 0.6Δlog K, after the range of studied compounds was extended to include analogues that were doubly protonated on γ-phosphate, CH3O-P(O2)-O-P(O2)-X-PO3H22–. The scissile Pα–Olg bond length in studied methyl triphosphate analogues slightly increases with decreasing pKa of the leaving group; concomitantly, the CH3OPα(O2) moiety becomes more positive. These structural effects indicate that substituents with low pKa can facilitate both Pα–Olg bond breaking and the Pα–Onuc bond forming process, thus explaining the large negative βlg calculated for the transition state geometry that has significantly longer Pα–Onuc distance than the Pα–Olg distance.

Semimicroscopic investigation of active site pK a values in peptidylarginine deiminase 4

Peptidylarginine deiminase 4 (PAD4), also known as protein arginine deiminase 4, performs a post-translational deimination that converts arginine to citrulline. The dysregulation of PAD4 has been implicated in a number of diseases, including rheumatoid arthritis (RA) and cancer. This makes PAD4 an important therapeutic target. To develop small-molecule inhibitors as potential treatments, it is advantageous if the catalytic mechanism is well understood. The protonation states of the active site residues, which have long been under controversy, have a direct impact on the catalytic mechanism. Two competing mechanisms are under investigation in the current literature. The first is a reverse protonation mechanism that depends on the active site histidine and cysteine existing as an ion pair. The second is a substrate-assisted mechanism that depends on the active site histidine and cysteine being neutral. This study uses the semimicroscopic protein dipoles Langevin dipoles (PDLD/S) linear response approximation method in the MOLARIS software package to calculate the change in solvation energy of moving the residue from water to the protein interior, and then using that information to assess the protonation states of the active site residues of PAD4. Results from these calculations suggest that in the enzyme–substrate complex of PAD4, the cysteine and histidine are protonated and deprotonated, respectively, and are therefore both neutral, analogous to the proposed protonation states of the active site residues in the Michaelis complex in the substrate-assisted mechanism.

Can the chemical reactivity of an ultimate carcinogen be related to its carcinogenicity? An application to propylene oxide

In this article we report calculations of the activation free energy for a chemical reaction between propylene oxide and DNA, in particular with the guanine at the N7 position. Calculations were performed using Hartree–Fock and MP2 methods in conjunction with flexible basis sets. The effects of solvation were considered using the Langevin dipoles and solvent reaction field methods. The calculated activation free energies are in good agreement with the experimental value of 25.36kcal/mol. Moreover, the reactivities of a series of ultimate carcinogens of the epoxy type toward DNA are shown to be related to their carcinogenicities. Higher reactivity is generally associated with higher carcinogenicity, although transport properties, reactions with ultimate carcinogen scavengers, and the DNA correction mechanism are also very important. It is very likely that the window of reactivity rather than a high reactivity value is the relevant measure of carcinogenicity, since highly reactive ultimate carcinogens interact with water and proteins before they reach DNA.

Validating the vitality strategy for fighting drug resistance

The current challenge in designing effective drugs against HIV-1 is to find novel candidates with high potency, but with a lower susceptibility to mutations associated with drug resistance. Trying to address this challenge, we developed in our previous study (Ishikita and Warshel, Angew Chem Int Ed Engl 2008; 47:697-700) a novel computational strategy for fighting drug resistance by predicting the likely moves of the virus through constraints on binding and catalysis. This has been based on calculating the ratio between the vitality values ((K(i) k(cat)/K(M))(mutant)/(K(i) k(cat)/K(M))(wild-type)) and using it as a guide for predicting the moves of the virus. The corresponding calculations of the binding affinity, K(i), were carried out using the semi-macroscopic version of the protein dipole Langevin dipole (PDLD/S) in its linear response approximation (LRA) in its β version (PDLD/S-LRA/β). We also calculate the proteolytic efficiency, k(cat)/K(M), by evaluating the transition state (TS) binding free energies using the PDLD/S-LRA/β method. Here we provide an extensive validation of our strategy by calculating the vitality of six existing clinical and experimental drug candidates. It is found that the computationally determined vitalities correlate reasonably well with those derived from the corresponding experimental data. This indicates that the calculated vitality may be used to identify mutations that would be most effective for the survival of the virus. Thus, it should be possible to use our approach in screening for mutations that would provide the most effective resistance to any proposed antiviral drug. This ability should be very useful in guiding the design of drug molecules that will lead to the slowest resistance.

Chemical Reaction of Soybean Flavonoids with DNA: A Computational Study Using the Implicit Solvent Model

Genistein, daidzein, glycitein and quercetin are flavonoids present in soybean and other vegetables in high amounts. These flavonoids can be metabolically converted to more active forms, which may react with guanine in the DNA to form complexes and can lead to DNA depurination. We assumed two ultimate carcinogen forms of each of these flavonoids, diol epoxide form and diketone form. Density functional theory (DFT) and Hartree-Fock (HF) methods were used to study the reaction thermodynamics between active forms of flavonoids and DNA guanine. Solvent reaction field method of Tomasi and co-workers and the Langevin dipoles method of Florian and Warshel were used to calculate the hydration free energies. Activation free energy for each reaction was estimated using the linear free energy relation. Our calculations show that diol epoxide forms of flavonoids are more reactive than the corresponding diketone forms and are hence more likely flavonoid ultimate carcinogens. Genistein, daidzein and glycitein show comparable reactivity while quercetin is less reactive toward DNA.

Reactivity of bisphenol A-3,4-quinone with DNA. A quantum chemical study

Bisphenol A is a primary component of polycarbonate plastics found in many different products – from reusable drink bottles to cell phones. It is also an important component of epoxy resins, which serve as a protective layer inside food and drink cans. Chronic exposure to bisphenol A from food, drink and other sources is the reason that this chemical shows up at low levels in the urine of nearly everyone. Bisphenol A is a known endocrine disruptor with a variety of other effects, including genotoxicity. These facts have created a lot of concern about how toxic it is. To investigate the possible genotoxic mechanisms of bisphenol A we calculated the chemical reactivity of the bisphenol A metabolite bisphenol A-3,4-quinone with deoxyguanosine by using density functional theory in conjunction with Langevin dipoles solvation model. The calculated activation free energy of 23.1kcal/mol shows that bisphenol A-3,4-quinone could form an adduct with deoxyguanosine and could be a mutagen. The subsequent depurination reaction was also studied with the same methodology. The calculated activation energy was 22.5kcal/mol, providing evidence that the rate-limiting step essential to causing genotoxicity is nucleophilic addition of the N7-deoxyguanosine to bisphenol A-3,4-quinone rather than depurination.

Impact of Proximal and Distal Pocket Site-Directed Mutations on the Ferric/Ferrous Heme Redox Potential of Yeast Cytochrome-c-Peroxidase.

Cytochrome-c-peroxidase (CCP) contains a five-coordinate heme active site. The reduction potential for the ferric to ferrous couple in CCP is anomalously low and pH dependent (Eo = ~-180 mV vs. S.H.E. at pH 7). The contribution of the protein environment to the tuning of the redox potential of this couple is evaluated using site directed mutants of several amino acid residues in the environment of the heme. These include proximal pocket mutation to residues Asp-235, Trp-191, Phe-202 and His-175, distal pocket mutation to residues Trp-51, His-52, and Arg-48; and a heme edge mutation to Ala-147. Where unknown, the structural changes resulting from the amino acid substitution have been studied by X-ray crystallography. In most cases, ostensibly polar or charged residues are replaced by large hydrophobic groups or alternatively by Ala or Gly. These latter have been shown to generate large, solvent filled cavities. Reduction potentials are measured as a function of pH by spectroelectrochemistry. Starting with the X-ray derived structures of CCP and the mutants, or with predicted structures generated by Molecular Dynamics (MD), predictions of redox potential changes are modeled using the Protein Dipoles Langevin Dipoles (PDLD) method. These calculations serve to model an electrostatic assessment of the redox potential change with simplified assumptions about heme iron chemistry, with the balance of the experimentally observed shifts in redox potential being thence attributed to changes in the ligand set and heme coordination chemistry, and/or other changes in the structure not directly evident in the X-ray structures (e.g. ionization states, specific roles played by solvent species, or conformationally flexible portions of the protein). Agreement between theory and experiment is good for all mutant proteins with the exception of the mutation Arg 48 to Ala, and His 52 to Ala. In the former case, the influence of phosphate buffer is adduced to account for the discrepancy, and measurements made in a bis-tris propane/2-(N-morpholino)ethanesulfonic acid buffer system agree well with theory. For the latter case, an unknown structural element relevant to His-52, and/or solvent influence in the mutant akin to anion binding in the distal pocket (though lacking proof that it is) manifests in this mutant. The use of exogenous (sixth) ligands in dissecting the contributions to control of redox potential are also explored as a pathway for model building.

Theoretical Comparison ofp-Nitrophenyl Phosphate and Sulfate Hydrolysis in Aqueous Solution: Implications for Enzyme-Catalyzed Sulfuryl Transfer

Both phosphoryl and sulfuryl transfers are ubiquitous in biology, being involved in a wide range of processes, ranging from cell division to apoptosis. Additionally, it is becoming increasingly clear that enzymes that can catalyze phosphoryl transfer can often cross-catalyze sulfuryl transfer (and vice versa). However, while there have been extensive experimental and theoretical studies performed on phosphoryl transfer, the body of available research on sulfuryl transfer is comparatively much smaller. The present work presents a direct theoretical comparison of p-nitrophenyl phosphate and sulfate monoester hydrolysis, both of which are considered prototype systems for probing phosphoryl and sulfuryl transfer, respectively. Specifically, free energy surfaces have been generated using density functional theory, by initial geometry optimization in PCM using the 6-31+G* basis set and the B3LYP density functional, followed by single-point calculations using the larger 6-311+G** basis set and the COSMO continuum model. The resulting surfaces have been then used to identify the relevant transition states, either by further unconstrained geometry optimization or from the surface itself where possible. Additionally, configurational entropies were evaluated using a combination of the quasiharmonic approximation and the restraint release approach and added to the calculated activation barriers as a correction. Finally, the overall activation entropy was estimated by approximating the solvent contribution to the total activation entropy using the Langevin dipoles solvation model. We have reproduced both the experimentally observed activation barriers and the observed trend in the activation entropies with reasonable accuracy, as well as providing a comparison of calculated and observed (15)N and (18)O kinetic isotope effects. We demonstrate that, counterintuitively, the hydrolysis of the p-nitrophenyl sulfate proceeds through a more expansive pathway than its phosphate analogue. Additionally, we show that the solvation effects upon moving from the ground state to the transition state are quite different for both reactions, suggesting that the enzymes that catalyze these reactions would need active sites with quite different electrostatic preorganization for the efficient catalysis of either reaction (despite which many enzymes can catalyze both phosphoryl and sulfuryl transfer). We believe that such a comparative study is an important foundation for understanding the molecular basis for phosphate-sulfate cross-promiscuity within members of the alkaline phosphatase superfamily.

Langevin Poisson-Boltzmann equation: point-like ions and water dipoles near a charged surface

Water ordering near a charged membrane surface is important for many biological processes such as binding of ligands to a membrane or transport of ions across it. In this work, the mean-field Poisson-Boltzmann theory for point-like ions, describing an electrolyte solution in contact with a planar charged surface, is modified by including the orientational ordering of water. Water molecules are considered as Langevin dipoles, while the number density of water is assumed to be constant everywhere in the electrolyte solution. It is shown that the dielectric permittivity of an electrolyte close to a charged surface is decreased due to the increased orientational ordering of water dipoles. The dielectric permittivity close to the charged surface is additionally decreased due to the finite size of ions and dipoles.

Excluded volume effect and orientational ordering near charged surface in solution of ions and Langevin dipoles

The influence of a finite volume of ions and orientational ordering of water Langevin dipoles on the dielectric permittivity profile in the vicinity of charged surface is studied theoretically via a numerical solution of the modified Poisson–Boltzmann equation. It is shown that the dielectric permittivity profile close to the charged surface is mainly determined by two mechanisms; specifically, the depletion of dipoles at the charged surface due to accumulated counterions and the increased orientational ordering of the water dipoles.

Absolute binding free energy calculations: On the accuracy of computational scoring of protein–ligand interactions

Calculating the absolute binding free energies is a challenging task. Reliable estimates of binding free energies should provide a guide for rational drug design. It should also provide us with deeper understanding of the correlation between protein structure and its function. Further applications may include identifying novel molecular scaffolds and optimizing lead compounds in computer-aided drug design. Available options to evaluate the absolute binding free energies range from the rigorous but expensive free energy perturbation to the microscopic linear response approximation (LRA/beta version) and related approaches including the linear interaction energy (LIE) to the more approximated and considerably faster scaled protein dipoles Langevin dipoles (PDLD/S-LRA version) as well as the less rigorous molecular mechanics Poisson-Boltzmann/surface area (MM/PBSA) and generalized born/surface area (MM/GBSA) to the less accurate scoring functions. There is a need for an assessment of the performance of different approaches in terms of computer time and reliability. We present a comparative study of the LRA/beta, the LIE, the PDLD/S-LRA/beta, and the more widely used MM/PBSA and assess their abilities to estimate the absolute binding energies. The LRA and LIE methods perform reasonably well but require specialized parameterization for the nonelectrostatic term. The PDLD/S-LRA/beta performs effectively without the need of reparameterization. Our assessment of the MM/PBSA is less optimistic. This approach appears to provide erroneous estimates of the absolute binding energies because of its incorrect entropies and the problematic treatment of electrostatic energies. Overall, the PDLD/S-LRA/beta appears to offer an appealing option for the final stages of massive screening approaches.

Calculation of p Ka values of carboxylic acids: Application to bilirubin

One of the major factors in the interactions of proteins with other molecules is the dissociation constant for a weakly acidic or basic group, generally expressed as the p K a. We have critically evaluated various quantum mechanical (QM) methods for calculation of protonation states of carboxylic acids including bilirubin. We calculated p K a values of various carboxylic acids by using quantum-chemical calculations in conjunction with Langevin dipoles and solvent reaction field method of Tomasi and co-workers. We studied formic acid, acetic acid, nicotinic acid and bilirubin. We demonstrated that the calculated p K a values compare well with the experiment if flexible basis sets of at least 6-31+G(d,p) size and Langevin dipoles solvation model are used. Results are discussed in terms of bilirubin transport by bilitranslocase.

Histamine (re)uptake by astrocytes: an experimental and computational study

Astrocytes participate in the clearance of neurotransmitters by their uptake and subsequent enzymatic degradation. Histamine as a polar and/or protonated molecule must use a carrier to be transported across the cell membrane, although a specific histamine transporter has not been elucidated, yet. In this work we upgraded the kinetic studies of histamine uptake into neonatal rat cultured type 1 astrocytes with quantum chemical calculations of histamine pKa values in conjunction with Langevin dipoles solvation model as the first step toward microscopic simulation of transport. Our results indicate that astrocytes transport histamine by at least two carrier mediated processes, a concentration gradient dependent passive and a sodium-dependent and ATP-driven active transport. We also demonstrated that histamine protonation states depend on the polarity of the environment. In conclusion we suggest that histamine, a polar molecule at physiological pH uses at least two different mechanisms for its uptake into astrocytes -an electrodiffusion and Na(+)-dependent and ouabain sensitive active process. We emphasize relevance of knowledge of histamines protonation states at the rate limiting step of its transport for microscopic simulation that will be possible when structure of histamine transporter is known.

A Computational Analysis of ATP Binding of SV40 Large Tumor Antigen Helicase Motor

Simian Virus 40 Large Tumor Antigen (LTag) is an efficient helicase motor that unwinds and translocates DNA. The DNA unwinding and translocation of LTag is powered by ATP binding and hydrolysis at the nucleotide pocket between two adjacent subunits of an LTag hexamer. Based on the set of high-resolution hexameric structures of LTag helicase in different nucleotide binding states, we simulated a conformational transition pathway of the ATP binding process using the targeted molecular dynamics method and calculated the corresponding energy profile using the linear response approximation (LRA) version of the semi-macroscopic Protein Dipoles Langevin Dipoles method (PDLD/S). The simulation results suggest a three-step process for the ATP binding from the initial interaction to the final tight binding at the nucleotide pocket, in which ATP is eventually “locked” by three pairs of charge-charge interactions across the pocket. Such a “cross-locking” ATP binding process is similar to the binding zipper model reported for the F1-ATPase hexameric motor. The simulation also shows a transition mechanism of Mg2+ coordination to form the Mg-ATP complex during ATP binding, which is accompanied by the large conformational changes of LTag. This simulation study of the ATP binding process to an LTag and the accompanying conformational changes in the context of a hexamer leads to a refined cooperative iris model that has been proposed previously.

Effective approach for calculations of absolute stability of proteins using focused dielectric constants

The ability to predict the absolute stability of proteins based on their corresponding sequence and structure is a problem of great fundamental and practical importance. In this work, we report an extensive, refinement and validation of our recent approach (Roca et al., FEBS Lett 2007;581:2065-2071) for predicting absolute values of protein stability DeltaG(fold). This approach employs the semimacroscopic protein dipole Langevin dipole method in its linear response approximation version (PDLD/S-LRA) while using the best fitted values of the dielectric constants epsilon'(p) and epsilon'(eff) for the self energy and charge-charge interactions, respectively. The method is validated on a diverse set of 45 proteins. It is found that the best fitted values of both dielectric constants are around 40. However, the self energy of internal residues and the charge-charge interactions of Lys have to be treated with care, using a somewhat lower values of epsilon'(p) and epsilon'(eff). The predictions of DeltaG(fold) reported here, have an average error of only 1.8 kcal/mole compared to the observed values, making our method very promising for estimating protein stability. It also provides valuable insight into the complex electrostatic phenomena taking place in folded proteins.