Weak Self-association Research Articles

A scaling relationship between thermodynamic and hydrodynamic interactions in protein solutions

Weak protein interactions are associated with a broad array of biological functions and are often implicated in molecular dysfunction accompanying human disease. In addition, these interactions are a critical determinant in the effective manufacturing, stability, and administration of biotherapeutic proteins. Despite their prominence, much remains unknown about how molecular attributes influence the hydrodynamic and thermodynamic contributions to the overall interaction mechanism. To systematically probe these contributions, we have evaluated self-interaction in a diverse set of proteins that demonstrate a broad range of behaviors from attractive to repulsive. Analysis of the composite trending in the data provides a convenient interconversion among interaction parameters measured from the concentration dependence of the molecular weight, diffusion coefficient, and sedimentation coefficient, as well as insight into the relationship between thermodynamic and hydrodynamic interactions. We find relatively good agreement between our data and a model for interacting hard spheres in the range of weak self-association. In addition, we propose an empirically derived, general scaling relationship applicable across a broad range of self-association and repulsive behaviors.

Characterization of interactions within the Igα/Igβ transmembrane domains of the human B-cell receptor provides insights into receptor assembly

The B-cell receptor (BCR), a complex comprised of a membrane-associated immunoglobulin and the Igα/β heterodimer, is one of the most important immune receptors in humans and controls B-cell development, activity, selection, and death. BCR signaling plays key roles in autoimmune diseases and lymphoproliferative disorders, yet, despite the clinical significance of this protein complex, key regions (i.e., the transmembrane domains) have yet to be structurally characterized. The mechanism for BCR signaling also remains unclear and has been variously described by the mutually exclusive cross-linking and dissociation activation models. Common to these models is the significance of local plasma membrane composition, which implies that interactions between BCR transmembrane domains (TMDs) play a role in receptor functionality. Here we used an in vivo assay of TMD oligomerization called GALLEX alongside spectroscopic and computational methods to characterize the structures and interactions of human Igα and Igβ TMDs in detergent micelles and natural membranes. We observed weak self-association of the Igβ TMD and strong self-association of the Igα TMD, which scanning mutagenesis revealed was entirely stabilized by an E–X10–P motif. We also demonstrated strong heterotypic interactions between the Igα and Igβ TMDs both in vitro and in vivo, which scanning mutagenesis and computational models suggest is multiconfigurational but can accommodate distinct interaction sites for self-interactions and heterotypic interactions of the Igα TMD. Taken together, these results demonstrate that the TMDs of the human BCR are sites of strong protein–protein interactions that may direct BCR assembly, endoplasmic reticulum retention, and immune signaling.

Polymer solubility in ionic liquids: dominated by hydrogen bonding.

Polymer solubility in ionic liquids (ILs) cannot be predicted by the solubility parameter approach based on the "like dissolves like" principle. According to the Kamlet-Abraham-Taft (KAT) multi-parameter polarity scale, ILs can be categorized on the basis of hydrogen-bond acidity or basicity ones. The experimental observations, that acidic ILs easily dissolve basic polymers and basic ILs dissolve acidic polymers, reflect the complementary nature of hydrogen-bonding interactions. A quantitative hydrogen-bonding analysis is proposed for predicting the solubility by taking the product of ΔαΔβ as an indicator of the competition between cross-association and self-association hydrogen bonding (H-bonding), where Δα is the difference of acidity parameters between the polymer and IL, and Δβ is the difference of basicity. This solubility criterion has been validated by the solubility data of 19 polymers (11 acidic and 8 basic) in 11 ILs (7 acidic and 4 basic). These principles based on KAT parameters can be applied to other systems dominated by hydrogen bonding.

Biochemical and structural analyses reveal that the tumor suppressor neurofibromin (NF1) forms a high-affinity dimer

Biochemical and structural analyses reveal that the tumor suppressor neurofibromin (NF1) forms a high-affinity dimer

Measuring Ultra-Weak Protein Self-Association by Non-ideal Sedimentation Velocity.

Ultra-weak self-association can govern the macroscopic solution behavior of concentrated macromolecular solutions ranging from food products to pharmaceutical formulations and the cytosol. For example, it can promote dynamic assembly of multi-protein signaling complexes, lead to intracellular liquid–liquid phase transitions, and seed crystallization or pathological aggregates. Unfortunately, weak self-association is technically extremely difficult to study, as it requires very high protein concentrations where short intermolecular distances cause strongly correlated particle motion. Additionally, protein samples near their solubility limit in vitro frequently show some degree of polydispersity. Here we exploit the strong mass-dependent separation of assemblies in the centrifugal field to study ultra-weak binding, using a sedimentation velocity technique that allows us to determine particle size distributions while accounting for colloidal hydrodynamic interactions and thermodynamic non-ideality (Chaturvedi, S. K.; et al. Nat. Commun.2018, 9, 4415; DOI: 10.1038/s41467-018-06902-x). We show that this approach, applied to self-associating proteins, can reveal a time-average association state for rapidly reversible self-associations from which the free energy of binding can be derived. The method is label-free and allows studying mid-sized proteins at millimolar protein concentrations in a wide range of solution conditions. We examine the performance of this method with hen egg lysozyme as a model system, reproducing its well-known ionic-strength-dependent weak self-association. The application to chicken γS-crystallin reveals weak monomer–dimer self-association with KD = 24 mM, corresponding to a standard free energy change of approximately −9 kJ/mol, which is a large contribution to the delicate balance of forces ensuring eye lens transparency.

Self-Assembly Controls Reactivity with Nitric Oxide: Implications for Fluorescence Sensing.

Three molecules containing the fluorophore 4-amino-1,8-naphthalimide (ANI) and showing different tendencies to self-assembly in aqueous environment have been prepared and fully characterized. The fluorescence emissions of two of these compounds in aqueous solutions are efficiently quenched in the presence of nitric oxide (NO) in aerated medium. Nuclear magnetic resonance and mass spectrometry techniques indicate that NO/O2 induces deamination of the ANI fluorophore, resulting in nonemissive 1,8-naphtalimide derivatives. It is found that the reactivity toward NO/O2 is regulated by the different aggregation modes presented by the molecules in aqueous medium. In this way, the molecules displaying fluorescence response toward NO/O2 are those with weak self-association properties whereas the compound with a high hydrophobic character (self-assembling into large nanoparticles) is insensitive to this species. Ultimately, the results described here could not only set the basis for the design of fluorescent bioprobes for NO/O2 based on ANI derivatives or other monoamino compounds but also could raise awareness about the importance of supramolecular interactions for the design of chemosensors.

Effects of Doxorubicin on the Liquid-Liquid Phase Change Properties of Elastin-Like Polypeptides

The lower critical solution temperature (LCST) of the thermo-responsive engineered elastin-like polypeptide (ELP) biopolymer is being exploited for the thermal targeted delivery of doxorubicin (Dox) to solid tumors. We examine the impact of Dox labeling on the thermodynamic and hydrodynamic behavior of an ELP drug carrier and how Dox influences the liquid-liquid phase separation (LLPS). Turbidity, dynamic light scattering (DLS), and differential scanning calorimetry measurements show that ELP undergoes a cooperative liquid-liquid phase separation from a soluble to insoluble coacervated state that is enhanced by Dox labeling. Circular dichroism measurements show that below the LCST ELP consists of both random coils and temperature-dependent β-turn structures. Labeling with Dox further enhances β-turn formation. DLS measurements reveal a significant increase in the hydrodynamic radius of ELP below the LCST consistent with weak self-association. Dox-labeled SynB1-ELP1 (Dox-ELP) has a significant increase in the hydrodynamic radius by DLS measurements that is consistent with stable oligomers and, at high Dox-ELP concentrations, micelle structures. Enhanced association by Dox-ELP is confirmed by sedimentation velocity analytical ultracentrifugation measurements. Both ELP self-association and the ELP inverse phase transition are entropically driven with positive changes in enthalpy and entropy. We show by turbidity and DLS that the ELP phase transition is monophasic, whereas mixtures of ELP and Dox-ELP are biphasic, with Dox-labeled ELP phase changing first and unlabeled ELP partitioning into the coacervate as the temperature is raised. DLS reveals a complex growth in droplet sizes consistent with coalescence and fusion of liquid droplets. Differential scanning calorimetry measurements show a −11 kcal/mol change in enthalpy for Dox-ELP coacervation relative to the unlabeled ELP, consistent with droplet formation being stabilized by favorable enthalpic interactions. We propose that the ELP phase change is initiated by ELP self-association, enhanced by increased Dox-ELP oligomer and micelle formation and stabilized by favorable enthalpic interactions in the liquid droplets.

Coarse-Grained Modeling of Antibodies from Small-Angle Scattering Profiles.

Predicting the concentrated solution behavior for monoclonal antibodies requires developing and using minimal models to describe their shape and interaction potential. Toward this end, the small-angle X-ray scattering (SAXS) profiles for a monoclonal antibody (COE-03) have been measured under solution conditions chosen to produce weak self-association. The experiments are complemented with molecular simulations of a three-bead antibody model with and without interbead attraction. The scattering profile is extracted directly from the molecular simulation to avoid using the decoupling approximation. We examine the ability of the three-bead model to capture features of the scattering profile and the dependence of compressibilty on protein concentration. The three-bead model is able to reproduce generic features of the experimental structure factor as a function of wave vector S(k) including a well-defined shoulder, which is a consequence of the planar structure of the antibody, and a well-defined minimum in S(k) at k ∼ 0.025 Å-1. We also show the decoupling approximation is incapable of accounting for highly anisotropic shapes. The best-fit parameters obtained from matching spherical models to simulated scattering profiles are protein concentration dependent, which limits their applicability for predicting thermodynamic properties. Nevertheless, the experimental compressibility curves can be accurately reproduced by an appropriate parametrization of the Baxter adhesive model, indicating the model provides a semiempirical equation of state for the antibody. The results provide insights into how equations of state can be improved for antibodies by accounting for their anisotropic shapes.

Weak self-association of cytochrome c peroxidase molecules observed by paramagnetic NMR.

There is growing experimental evidence that many proteins exhibit a tendency for (ultra)weak homo- or hetero- oligomerization interactions. With the development of paramagnetic relaxation enhancement NMR spectroscopy it has become possible to characterize weak complexes experimentally and even detect complexes with affinities in the 1–25 mM range. We present evidence for a weak complex between cytochrome c peroxidase (CcP) molecules. In a previous study, we attached nitroxide based spin labels at three positions on CcP with the intent of observing intramolecular PRE effects. However, several intermolecular PRE effects were also observed suggesting a weak self-association between CcP molecules. The CcP–CcP complex was characterized using paramagnetic NMR and protein docking. The interaction occurs between the surface that is also part of the stereo-specific binding site for its physiological partner, cytochrome c (Cc), and several small, positively charged patches on the “back” of CcP. The CcP–CcP complex is not a stereo-specific complex. It is a dynamic ensemble of orientations, characteristic of an encounter state. The contact areas resemble those observed for CcP molecules in crystals. The CcP–CcP complex formation competes with that of the CcP-Cc complex. However, the affinity for Cc is much larger and thus it is expected that, under physiological conditions, auto-inhibition will be limited.Graphical A weak self-association between cytochrome c peroxidase molecules was characterized using paramagnetic NMR.Electronic supplementary materialThe online version of this article (doi:10.1007/s10858-016-0035-z) contains supplementary material, which is available to authorized users.

The Ras G Domain Lacks the Intrinsic Propensity to Form Dimers

Ras GTPase is a molecular switch controlling a number of cellular pathways including growth, proliferation, differentiation, and apoptosis. Recent reports indicated that Ras undergoes dimerization at the membrane surface through protein-protein interactions. If firmly established this property of Ras would require profound reassessment of a large amount of published data and modification of the Ras signaling paradigm. One proposed mechanism of dimerization involves formation of salt bridges between the two GTPase domains (G domains) leading to formation of a compact dimer as observed in Ras crystal structures. In this work, we interrogated the intrinsic ability of Ras to self-associate in solution by creating conditions of high local concentration through irreversibly tethering the two G domains together at their unstructured C-terminal tails. We evaluated possible self-association in this inverted tandem conjugate via analysis of the time-domain fluorescence anisotropy and NMR chemical shift perturbations. We did not observe the increased rotational correlation time expected for the G domain dimer. Variation of the ionic strength (to modulate stability of the salt bridges) did not affect the rotational correlation time in the tandem further supporting independent rotational diffusion of two G domains. In a parallel line of experiments to detect and map weak self-association of the G domains, we analyzed NMR chemical shifts perturbations at a number of sites near the crystallographic dimer interface. The nearly complete lack of chemical shift perturbations in the tandem construct supported a simple model with the independent G domains repelled from each other by their overall negative charge. These results lead us to the conclusion that self-association of the G domains cannot be responsible for homodimerization of Ras reported in the literature.

Acyclic cucurbit[n]uril-type molecular containers: influence of aromatic walls on their function as solubilizing excipients for insoluble drugs.

We studied the influence of the aromatic sidewalls on the ability of acyclic CB[n]-type molecular containers (1a–1e) to act as solubilizing agents for 19 insoluble drugs including the developmental anticancer agent PBS-1086. All five containers exhibit good water solubility and weak self-association (Ks ≤ 624 M–1). We constructed phase solubility diagrams to extract Krel and Ka values for the container·drug complexes. The acyclic CB[n]-type containers generally display significantly higher Ka values than HP-β-CD toward drugs. Containers 1a–1e bind the steroidal ring system and aromatic moieties of insoluble drugs. Compound 1b displays highest affinity toward most of the drugs studied. Containers 1a and 1b are broadly applicable and can be used to formulate a wider variety of insoluble drugs than was previously possible with cyclodextrin technology. For drugs that are solubilized by both HP-β-CD and 1a–1e, lower concentrations of 1a–1e are required to achieve identical [drug].

The structure of α-haemoglobin in complex with a haemoglobin-binding domain from Staphylococcus aureus reveals the elusive α-haemoglobin dimerization interface.

Adult haemoglobin (Hb) is made up of two α and two β subunits. Mutations that reduce expression of the α- or β-globin genes lead to the conditions α- or β-thalassaemia, respectively. Whilst both conditions are characterized by anaemia of variable severity, other details of their pathophysiology are different, in part owing to the greater stability of the β chains that is conferred through β self-association. In contrast, α subunits interact weakly, and in the absence of stabilizing quaternary interactions the α chain (α) is prone to haem loss and denaturation. The molecular contacts that confer weak self-association of α have not been determined previously. Here, the first structure of an α2 homodimer is reported in complex with one domain of the Hb receptor from Staphylococcus aureus. The α2 dimer interface has a highly unusual, approximately linear, arrangement of four His side chains within hydrogen-bonding distance of each other. Some interactions present in the α1β1 dimer interface of native Hb are preserved in the α2 dimer. However, a marked asymmetry is observed in the α2 interface, suggesting that steric factors limit the number of stabilizing interactions that can form simultaneously across the interface.

Specific and nonspecific interactions in ultraweak protein-protein associations revealed by solvent paramagnetic relaxation enhancements.

Weak and transient protein–protein interactions underlie numerous biological processes. However, the location of the interaction sites of the specific complexes and the effect of transient, nonspecific protein–protein interactions often remain elusive. We have investigated the weak self-association of human growth hormone (hGH, KD = 0.90 ± 0.03 mM) at neutral pH by the paramagnetic relaxation enhancement (PRE) of the amide protons induced by the soluble paramagnetic relaxation agent, gadodiamide (Gd(DTPA-BMA)). Primarily, it was found that the PREs are in agreement with the general Hwang-Freed model for relaxation by translational diffusion (J. Chem. Phys.1975, 63, 4017–4025), only if crowding effects on the diffusion in the protein solution are taken into account. Second, by measuring the PREs of the amide protons at increasing hGH concentrations and a constant concentration of the relaxation agent, it is shown that a distinction can be made between residues that are affected only by transient, nonspecific protein–protein interactions and residues that are involved in specific protein–protein associations. Thus, the PREs of the former residues increase linearly with the hGH concentration in the entire concentration range because of a reduction of the diffusion caused by the transient, nonspecific protein–protein interactions, while the PREs of the latter residues increase only at the lower hGH concentrations but decrease at the higher concentrations because of specific protein–protein associations that impede the access of gadodiamide to the residues of the interaction surface. Finally, it is found that the ultraweak aggregation of hGH involves several interaction sites that are located in patches covering a large part of the protein surface.

Structural and Hydrodynamic Analysis of a Novel Drug Delivery Vector: ELP[V5G3A2-150

The therapeutic potential of elastin-like polypeptide (ELP) conjugated to therapeutic compounds is currently being investigated as an approach to target drugs to solid tumors. ELPs are hydrophobic polymers that are soluble at low temperatures and cooperatively aggregate above a transition temperature (TT), allowing for thermal targeting of covalently attached drugs. They have been shown to cooperatively transition from a disordered structure to a repeating type II β-turn structure, forming a β-spiral above the TT. Here we present biophysical measurements of the structural, thermodynamic, and hydrodynamic properties of a specific ELP being investigated for drug delivery, ELP[V5G3A2-150]. We examine the biophysical properties below and above the TT to understand and predict the therapeutic potential of ELP-drug conjugates. We observed that below the TT, ELP[V5G3A2-150] is soluble, with an extended conformation consisting of both random coil and heterogeneous β structures. Sedimentation velocity experiments indicate that ELP[V5G3A2-150] undergoes weak self-association with increasing temperature, and above the TT the hydrophobic effect drives aggregation entropically. These experiments also reveal a previously unreported temperature-dependent critical concentration (Cc) that resembles a solubility constant. Labeling ELP[V5G3A2-150] with fluorescein lowers the TT by 3.5°C at 20 μM, whereas ELP[V5G3A2-150] dissolution in physiological media (fetal bovine serum) increases the TT by ∼2.2°C.

The molecular basis of ATM-dependent dimerization of the Mdc1 DNA damage checkpoint mediator

Mdc1 is a large modular phosphoprotein scaffold that maintains signaling and repair complexes at double-stranded DNA break sites. Mdc1 is anchored to damaged chromatin through interaction of its C-terminal BRCT-repeat domain with the tail of γH2AX following DNA damage, but the role of the N-terminal forkhead-associated (FHA) domain remains unclear. We show that a major binding target of the Mdc1 FHA domain is a previously unidentified DNA damage and ATM-dependent phosphorylation site near the N-terminus of Mdc1 itself. Binding to this motif stabilizes a weak self-association of the FHA domain to form a tight dimer. X-ray structures of free and complexed Mdc1 FHA domain reveal a ‘head-to-tail’ dimerization mechanism that is closely related to that seen in pre-activated forms of the Chk2 DNA damage kinase, and which both positively and negatively influences Mdc1 FHA domain-mediated interactions in human cells prior to and following DNA damage.

Conformational Changes of the HIV-1 Envelope Protein during Membrane Fusion Are Inhibited by the Replacement of Its Membrane-spanning Domain

To help understand the dynamic nature of membrane fusion induced by the human immunodeficiency virus-1 (HIV-1) envelope protein, we developed a new cell-based real-time assay system employing a pair of novel reporter proteins. The reporter proteins consist of a pair of split Renilla luciferase (spRL) fused to split green fluorescent protein (spGFP). The spGFP modules were chosen not only to compensate weak self-association of spRL but also to provide visual reporter signals during membrane fusion. Use of this reporter together with a membrane permeable substrate for Renilla luciferase achieved a simple real-time monitoring of membrane fusion using live cells. We analyzed the HIV-1 envelope mutants whose membrane-spanning domains were replaced with that of glycophorin A or vesicular stomatitis virus G-protein. These mutants showed a slower kinetics of membrane fusion. The analysis of membrane fusion in the presence of fusion inhibitors, soluble CD4 and C34, revealed that these replacements prolonged the period during which the mutants were sensitive to the inhibitors, as compared with the wild type. These results suggest that the mutations within the membrane-spanning domains exerted an allosteric effect on the HIV-1 envelope protein, probably affecting the receptor-induced conformational changes of the ectodomain of the protein.

Synthesis and Self‐Aggregation of Enantiopure and Racemic Molecular Tweezers Based on the Bicyclo[3.3.1]nonane Framework

A pair of molecular tweezers (syn-4) that consists of quinoline and pyrazine units fused to a bicyclic framework is presented. The tweezers were synthesised both as a racemic mixture (rac-4) and an enantiomerically pure form ((R,R,R,R)-4) starting from either racemic or enantiomerically pure bicyclo[3.3.1]nonane-2,6-dione (3). Homochiral dimers were observed in the solid state for rac-4. The self-association of both rac-4 and (R,R,R,R)-4 was studied in solution. A weak self-association constant in CDCl(3) was estimated by (1)H NMR spectroscopic dilution titration experiments in both cases, following several proton resonances. For this purpose, a general normalisation model for the accurate determination of association constants from multiple datasets was developed. In contrast to the solid state, no diastereomeric discrimination was observed for rac-4 in solution.

Macromolecular Interaction of Halichondrin B Analogues Eribulin (E7389) and ER-076349 with Tubulin by Analytical Ultracentrifugation

Halichondrin B is an antimitotic drug that inhibits microtubule assembly. To understand the molecular details of its interaction with tubulin, we investigated the binding of two halichondrin B analogues, eribulin (previously, ER-086526, E7389) and ER-076349, to tubulin by quantitative analytical ultracentrifugation. Eribulin is currently undergoing phase III clinical trials for cancer; ER-076349 is a closely related analogue with C.35 hydroxyl instead of C.35 primary amine [Towle, M. J., et al. (2001) Cancer Res. 61, 1013]. Below the critical concentration for microtubule assembly and in the presence of GDP, tubulin undergoes weak self-association into short curved oligomers. Eribulin inhibits this oligomer formation 4-6-fold, while ER-076349 slightly stimulates oligomer formation by 2-fold. This is in contrast to vinblastine which strongly stimulates large spiral polymers by 1000-fold under these same conditions. Vinblastine-induced spiral formation is strongly inhibited by both eribulin and ER-076349. Colchicine binding to the intradimer interface has no significant effect on small oligomer formation or the inhibitory activity of eribulin on this process. These results suggest that halichondrin B analogues bind to the interdimer interface or to the beta-subunit alone, disrupt polymer stability, and compete with vinblastine-induced spiral formation. Stathmin is known to form a tight 1:2 complex with tubulin. Eribulin strongly inhibits formation of the 1:2 stathmin-tubulin complex (>3.3 kcal/mol), while ER-076349 weakens formation of the 1:2 complex by approximately 1.9 kcal/mol. These results suggest that eribulin is a global inhibitor of tubulin polymer formation, disrupting tubulin-tubulin contacts at the interdimer interface. ER-076349 also perturbs tubulin-tubulin contacts, but in a more polymer specific manner, reflecting adaptability of the interdimer interface to drug and polymer polymorphism. These results suggest halichondrin B analogues exhibit unique tubulin-based activities that may underlie the clinical utility of these compounds.

Implications of the Progressive Self-association of Wild-type Human Factor H for Complement Regulation and Disease

Factor H (FH) is a major regulator of complement alternative pathway activation. It is composed of 20 short complement regulator (SCR) domains and is genetically associated as a risk factor for age-related macular degeneration. Previous studies on FH suggested that it existed in monomeric and dimeric forms. Improved X-ray scattering and analytical ultracentrifugation methodology for wild-type FH permitted a clarification of these oligomeric properties. Data at lower concentrations revealed a dependence of the X-ray radius of gyration values on concentration that corresponded to the weak self-association of FH. Global sedimentation equilibrium fits indicated that a monomer–dimer equilibrium best described the data up to 1.3 mg/ml with a fitted dissociation constant K D of 28 μM and that higher oligomers formed at increased concentrations. The K D showed that about 85–95% of serum FH will be monomeric in the absence of other factors. Size-distribution analyses in sedimentation velocity experiments showed that monomeric FH was the major species but that as many as six oligomeric forms co-existed with it. The data were explained in terms of two weak dimerisation sites recently identified in the SCR-6/8 and SCR-16/20 fragments of FH with similar K D values. These observations indicate a mechanism for the progressive self-association of FH and may be relevant for complement regulation and the formation of drusen deposits that are associated with age-related macular degeneration.

Weak Self-Association in a Carbohydrate System

The physiological importance of weak interactions between biological macromolecules (molar dissociation constants >10 μM) is now well recognized, particularly with regard to cell adhesion and immunological phenomena, and many weak interactions have been measured for proteins. The concomitant importance of carbohydrate-carbohydrate interactions has also been identified, although no weak interaction between pure carbohydrate systems has ever been measured. We now demonstrate for the first time to our knowledge using a powerful probe for weak interactions—sedimentation velocity in the analytical ultracentrifuge—that at least some carbohydrates (from the class of polysaccharides known as heteroxylans and demonstrated here to be biologically active) can show well-defined weak self-interactions of the “monomer-dimer” type frequently found in protein systems. The weak interaction between the heteroxylans is shown from a temperature dependence study to be likely to be hydrophobic in nature.