“Mirror”-like Protein Dimers Stabilized by Local ... - ACS Publications

The protein dimers form “mirror”-like orientations with some degree of rotation. Such signature orientations are interpreted based on the predomin...
2 downloads 0 Views 1MB Size
Subscriber access provided by UNIV AUTONOMA DE COAHUILA UADEC

B: Biophysics; Physical Chemistry of Biological Systems and Biomolecules

“Mirror”-like Protein Dimers Stabilized by Local Heterogeneity at Protein Surfaces Baofu Qiao, Luis G. Lopez, and Monica Olvera de la Cruz J. Phys. Chem. B, Just Accepted Manuscript • DOI: 10.1021/acs.jpcb.9b01394 • Publication Date (Web): 18 Apr 2019 Downloaded from http://pubs.acs.org on April 19, 2019

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 20 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

“Mirror”-like Protein Dimers Stabilized by Local Heterogeneity at Protein Surfaces Baofu Qiao, † Luis G. Lopez † and Monica Olvera de la Cruz *,†,‡,§ †Department of Materials Science and Engineering, ‡Department of Chemistry, and § Department of Physics and Astronomy, Northwestern University, Evanston, IL 60208

Abstract: Protein aggregation has been observed inside cells and holds true for membraneless organelles. The precise understanding of protein dimerization is a prerequisite for manipulating protein aggregation, which is promising for elevating enzyme concentration to enhance their catalytic performance. Here, the dimerization of two industrially important enzymes of cytochrome P450 (P450) and organophosphorus hydrolase (OPH) are investigated using all-atoms explicit solvent molecular dynamics (MD) simulations, umbrella sampling, and protein-protein docking calculations. The calculated potentials of mean force of dimer-monomer dissociation demonstrate that the dimeric forms are more stable with the free energy barrier of around 60 kJ/mol for P450 and 101 kJ/mol for OPH. The docking calculations on the OPH dimer evidence the uniqueness of the native orientation. The protein dimers form “mirror”-like orientations with some degree of rotation. Such signature orientations are interpreted based on the predominant polar amino acids in the contact regime. In the dimer conformations the active sites are exposed. This work highlights the crucial roles of the polar and nonpolar protein surface domains to form enzymatically active protein dimer aggregates. Our work will potentially aid the design of molecules that can deliver and protect native protein function in various environments.

ACS Paragon Plus Environment

1

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 20

INTRODUCTION Proteins are one of the most important components in biology. Recently, protein phase separation has been discovered inside cells.1-3 Protein aggregation is believed to be of relevance to numerous diseases (e.g., cataract,4 neurodegenerative diseases5). It also drives the formation of spatiotemporal aggregates termed membraneless organelles, 6-7 which enable cellular processes to take place. For example, Nott et al. found that membraneless organelles made of Ddx4N1 proteins selectively stabilized single-stranded nucleic acids.8 Recently, heteropolymer-protein hybrid membraneless organelles were assembled9 by exploiting the local heterogeneity of polar and nonpolar amino acids at protein surfaces and by designing random heteropolymers that form enzymatically active protein-polymer complexes. Interestingly, proteins with given distributions of surface domains were found to selectively interact with amphiphilic random copolymers with similar monomer sequence.10 That is, knowledge of the protein surface domains is essential to manipulate the aggregation behavior of proteins. Another extremely important biochemical process where the surface domain heterogeneity is critically important is protein aggregation into dimers. Protein dimers are common in biology for regulating enzyme activation, for the functions of membrane proteins and for DNA binding.11 The understanding of dimerization, for example, will aid us in manipulating enzyme oligomerization,11 which is promising for advancing catalytic performance by increasing enzyme concentrations. Here, we investigate the homodimerization of two industrially important enzymes, namely P45012 and organophosphorus hydrolase (OPH)13-14 (Figure 1). The dimerization has been reported to affect the catalytic activities of P45015 and OPH.16 The cytochrome P450s (P450s or CYPs for short) are a superfamily of one of nature’s most versatile biocatalysts.17 P450s have been found to play a pivotal role in human steroid metabolism, xenobiotic detoxification and metabolism. These enzymes can catalyze a broad variety of chemical reactions of different substrates (e.g., aromatic hydrocarbon hydroxylation and alkene epoxidation).18 Around 75% of drugs are metabolized by P450s.19 P450s show remarkable versatility with thousands of P450s reported in various organisms.17 Among them, P450BM3 (CYP102A1) is one of the most significant as it shows higher activity and catalytic self-sufficiency, since the HEME and diflavin reductase domains are fused together in one single molecule. 20-22 The aggregation of P450 has been recognized in both crystal and solution.23 Oligomers up to tetramer coexist for CYP126A1.24 The probability of larger oligomers increases with the increase in protein concentration from 10 to 20 µM, . Ligands can also influence the dissociation of CYP126A124 and CYP12125 dimers. The dominant distributions of P450BM3 dimers were supported by a variety of experimental techniques. 26-27 Remarkably, the P450BM3 dimeric form is argued ACS Paragon Plus Environment

2

Page 3 of 20 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

to be catalytically active for fatty acid hydroxylation, where electron transfer occurs between the reductase domain of one subunit and the HEME domain of the other.15

Figure 1. Snapshots of protein dimers (A) P450 (pdb code:1BU7) and (B) OPH (pdb code:1HZY).

In contrast to P450, OPH is much more specific. It is unique in hydrolyzing organophosphates by cleaving their P-X (X = S, O, F) covalent bonds.13-14 Being extremely toxic chemicals, organophosphates are the primary components of pesticides and chemical warfare agents (VX, sarin, soman, etc.). As such, these chemicals pose a huge threat to public safety. For instance, globally, around 38% of pesticides are composed of organophosphorus compounds, which result in around 3 million poisonings and 200,000 deaths annually.14 Therefore, the development of highly efficient methods for degrading organophosphate-based chemicals for on-site bioremediation is critical. 9 The oligomerization state of OPH is in debate. The dimeric form has been extensively supported by X-ray crystallization,28-29 ultracentrifugation28 and gel filtration measurement.16 The monomeric form was claimed in gel filtration measurements at varying ionic strengths (0 to 0.2 M KCl) with up to 2 mM Zn2+ and Co2+ ions.30 Note that the OPH dimer has been shown to be active, which can partially unfold to an intermediate, inactive dimeric form, and then to the dissociated monomeric form. 16 By coupling classical all-atom explicit solvent molecular dynamics (MD) simulations with steered MD simulations and protein-protein docking methods, we investigated the aggregation behavior of P450 and OPH (i.e., the dimer-monomer dissociation). The free energy calculations support the dimeric forms to be energetically stable for both proteins under the conditions investigated. The subunits in the dimeric forms are forming “mirror”-like structures. The protein orientations are explained on the basis of the local heterogeneity of polar and nonpolar amino acids at protein surfaces. These findings improve our ACS Paragon Plus Environment

3

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 20

understanding of condensing enzymes and they provide guidelines for developing a potential method for increasing enzyme concentration, and consequently enhancing their catalytic performance. COMPUTATIONAL DETAILS All-Atom Explicit Solvent MD Simulations Classical MD simulations were performed at the all-atom resolution using the package GROMACS (version 2016.4).31 The CHARMM36m force field 32 was used along with the recommended CHARMM TIP3P water model33 whose structures were constrained using the SETTLE algorithm.34 The initial structures of P450 and OPH were downloaded from the protein data bank webpage (www.rcsb.org) with the PDB ID of 1BU7 35 for P450 dimer (the HEME domain of P450BM3) and 1HZY 29

for OPH dimer. For both proteins, the ligands required for crystallization were removed. Some

mutations were performed for OPH to be consistent with the experimental sequence from our collaborator (unpublished): (a) Lys169 was carboxylated30, 36 to be Kcx169 by including the atoms of formic acids; (b) Glu159 was mutated to be Lys159; (c) the oxygen atoms which were coordinating two Zn2+ ions simultaneously were modified to be hydroxide ions (OH-);30, 36 (c) the capping amino acids AGSIGTG were added before Asp35. Each protein dimer was dissolved in a water box with the edge length of around 16 nm. The net charge of each monomer is -15e for P450 and +3e for OPH, where e stands for the elementary charge. Consequently, 30 Na+ for P450 and 6 Cl- for OPH were added to neutralize the protein dimers. Then 247 Na+ and Cl- ions were added, corresponding to a salt concentration of approximately 0.1 M. The potential energy of the system was minimized using the steepest descent algorithm, followed by an equilibration of 50 ns. The three-dimensional periodic boundary conditions were applied. The neighbor searching was done for a distance of 1.2 nm using the Verlet particle-based cut-off approach and was updated every 20 time steps. The short-range Lennard-Jones (LJ) 12-6 interactions were switched off from 1 to 1.2 nm. The short-range Coulomb interactions was truncated at 1.2 nm with the long-range interactions calculated by means of the Particle Mesh Ewald algorithm.37-38 By constraining all the covalent bonds using the LINCS algorithm,39-40 a time step of 2 fs was employed. The NTP ensemble (constant number of particles, temperature and pressure) was applied with the temperature coupled via the Nosé-Hover algorithm (reference temperature 298K, characteristic time 1 ps) and the isotropic Parrinello-Rahman barostat (reference pressure 1 bar, characteristic time 4 ps and compressibility 4.5×10-5 bar-1). The protein dimers were stable during the simulation duration of 50 ns based on the calculated density, the potential energy and the root-mean-square-deviation of protein backbone atoms ACS Paragon Plus Environment

4

Page 5 of 20 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(data not presented). In addition to the steerer MD and umbrella sampling calculations below, an unbiased control simulation of 800 ns also supported the stability of the protein dimer. In this control simulation, an OPH dimer was solvated in a smaller simulation box with the edge length of around 12 nm than the length of around 16 nm in the advanced sampling calculations to decrease the computational cost. Steered MD Simulations and Umbrella Sampling The steered simulations were subsequently employed to pull the protein dimers apart in mimicking the dimer-monomer dissociation process. Based on the generated snapshots, the umbrella sampling was applied to calculate the potential of mean force using the Weighted Histogram Analysis Method (WHAM) method. 41-42 In the steered simulations, the atoms, except hydrogen, of one subunit were restrained using a force constant of 1000 kJ/mol/nm2. The other subunit was pulled away as a function of the protein-protein center-of-mass distance, which approached half of the simulation box (~ 8 nm). A pulling velocity of 5 nm/ns was employed. In the steered simulations, snapshots were saved with an interval of around 1 Å within the distance range of 5.0 - 7.8 nm for P450 dimer, and 4.2 - 7.9 nm for OPH dimer. The umbrella sampling simulations were then performed for all the snapshots. Each system was first equilibrated 100 ps, followed by a production simulation of 20 ns with the first 1 ns abandoned for the subsequent data analysis. The WHAM method was employed to calculate the potential of mean force using the program g_wham.

42

The histogram distributions of umbrella sampling are presented in Figure S1, and they support the convergence of the calculations. The umbrella sampling was computationally very expensive due to the large simulations box (~ 16 nm in each dimension). Protein-Protein Docking Calculations for OPH Dimer In order to obtain a free energy landscape of the protein dimers, docking calculations were performed. Package Hex 8.0

43

was used to generate protein-protein complexes, and program FastContact

44

was

used to estimate both the electrostatic and desolvation energy of each complex. The construction of such a landscape was not possible for the P450 dimer due to the highly lumpy surface as well as the highly negatively charged feature of P450 (i.e., -15e per monomer). Consequently, in this work, the docking analyses were exclusively performed for the OPH dimer. Before the docking calculations, the capping amino acids from Ala28 to Gly34 were deleted from the OPH dimer, since these tails are incompatible with rigid body calculations. Kcx169 was replaced with Lys169 in both subunits due to the absence of the atomistic potential in the programs used. 10,000 ACS Paragon Plus Environment

5

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 20

protein dimer complexes were generated by means of a rigid docking of OPH subunits A and B. The only criteria to sort the 10,000 “artificial” dimers was the geometric complementarity (also known as surface or shape complementarity). Both the position and orientation of the subunit A (the reference subunit) were kept fixed, while all possible positions and orientations of the subunit B (the neighbor subunit), always in close contact with the subunit A, were scanned. The native OPH complex was used as the reference orientation to specify the relative orientation of the “artificial” complexes. In the reference orientation, the geometric center of the subunit A was at (0, 0, 0) in the Cartesian coordinates, while the geometric center of the subunit B was at (0, 0, 4.2 nm), where 4.2 nm was the average distance between the centers of the two proteins. The energies of the 10,000 complexes were minimized using the program Minimize.x of the Tinker package (version 8.4).45 The free energy of each complex was estimated using FastContact.

44

This

program rapidly estimates both the electrostatic and the desolvation contributions of the free energy. The former is the standard intermolecular Coulomb electrostatic potential, calculated using a distancedependent dielectric constant of 4.5r,

46

and the latter is based on an empirical contact

potential developed using a database of structures from the protein databank.47 The desolvation component takes into account the hydrophobic interactions, the change in energy due to the desolvation of polar and charged atoms, and the conformational entropy loss of the side chains. The total free energy is the sum of both contributions. Subsequently, all the protein-protein complexes with repulsive interactions (positive free energy) were discarded. For each position of protein B around protein A, only the complex with the strongest attractive interaction was kept, regardless of its orientation. 409 out of 10,000 complexes remained after this process, and these complexes were used to create the free energy landscape of the protein-protein interaction. RESULTS AND DISCUSSION Free Energy of Protein Dimerization To justify the stability of the protein dimers, we calculated the free energy of protein dimerization. In this respect, the protein-protein docking calculations were performed first, which provided the threedimensional free energy map of the possible protein dimers. Moreover, the potential of mean force of protein dimer-monomer dissociation was calculated for the native complex as a function of the proteinprotein distance. These two calculations combined provide us with a better understanding of the free energy involved in the dimerization. ACS Paragon Plus Environment

6

Page 7 of 20 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 2. (A) Top and (B) bottom hemispheres of the free energy profile of protein-protein docking for the OPH dimer. Only the attractive interactions (ΔG < 0) are displayed. The center of the reference subunit (in cyan in the top insert) is located at (0, 0, 0), with the neighbor subunit colored in silver. The three strongest orientations are labeled in panel (A), with their free energy values listed in the bottom inset.

Presented in Figure 2 are the three-dimensional distributions of the protein-protein orientations with the attractive free energies projected on the XY plane. It is shown that the free energy attractions are relatively weak (ΔG > -50 kJ/mol) for most of the orientations. The stronger attractions (ΔG < -100 kJ/mol) only occur at some orientations. The three orientations with the strongest free energy attractions are labeled with a, b and c in Figure 2, respectively. Their free energy values are listed in the bottom inset. The orientation at (0, 0, 4.2 nm) shows the strongest attractive free energy, which supports the favorable feature of the native orientation. In complement to the three-dimensional protein dimerization profile in Figure 2, the potential of mean force (PMF, ΔG(ξ)) demonstrates the one-dimensional free energy profile along the protein-protein distance (ξ) (Figure 3A). In what follows, the protein-protein distance refers to the distance between the center-of-mass of the subunits, which include the cofactors (HEME for P450 and two Zn2+ ions for OPH). The calculated results have been shifted so that the average free energy in the distance range of 6.5 ≤ ξ ≤ 7.5 nm is considered as the reference (zero). Consequently, the PMFs of the dimer-monomer dissociation were obtained, with ΔG = 60.9 ± 0.9 kJ/mol for P450 and 101.1 ± 0.7 kJ/mol for OPH, where the error bars refer to the standard deviations obtained from the corresponding PMF values within 6.5 ≤ ξ ≤ 7.5 nm. These results support that the dimeric form is more stable for both proteins. The calculated

ACS Paragon Plus Environment

7

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 20

free energy of 101 kJ/mol for OPH dimer reasonably agrees with the docking results (Figure 2). The dimeric form of OPH is energetically more stable than that of P450, which is expected to result from the less net charge of +3e per OPH monomer in comparison with -15e per P450 monomer. The obtained free energies are close to those of insulin dimer (30,48 40,49 50,50 126 kJ/mol 51), trypsin-pancreatic trypsin inhibitor complex (188 kJ/mol) 51 and hemoglobin αβ dimer (> 159 kJ/mol).51 Figure 3. (A) Potential of mean force of the dimer-monomer dissociation for P450 and OPH as a function of the protein-protein distance. Insets are the corresponding simulation snapshots for associated and A

B

dissociated configurations. (B) Solvent accessible surface area of protein dimers of P450 and OPH as a function of the protein-protein distance. The error bars in (B) stand for the standard deviations.

Note that some uncertainties exist in calculating the PMFs, e.g., the accuracy of the force field, the convergence of the umbrella sampling, the numerical integration of the WHAM calculation. In terms of the force field, the CHARMM36m potential is one of the most extensively employed and calibrated allatom potential for protein simulations.32 The simulation duration in the umbrella sampling calculations does affect the obtained values in a quantitative manner. For instance, a shorter simulation of 5 ns for each sampling window results in a slightly weaker value of around 50 kJ/mol for P450 and 91 kJ/mol for OPH, in comparison with the values of 60.9 kJ/mol for P450 and 101.1 kJ/mol for OPH with 20 ns for each sampling window. The convergence of the numerical integration of the WHAM calculations was supported by the overlap of the neighbor occurrence profiles (Figure S1). Additionally, some other approaches were proposed for protein-protein dissociation calculations.52 To further understand the energetic driving force of protein dimerization, we calculated the solvent accessible surface areas (SASA). As demonstrated in Figure 3B the SASAs increase as a function of the protein-protein distance. As the protein dimers dissociate, waters will be adsorbed at the protein surface ACS Paragon Plus Environment

8

Page 9 of 20 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

which was previously buried by the neighbor subunit. The corresponding entropy decreases for the water molecules because of a loss in the degrees of freedom. This kind of solvent entropy contribution is approximately 10 kJ/mol per 1 nm2 of protein surface area.51,

53

Therefore, the solvent entropy

contribution (-TΔS) is estimated to be around 190 kJ/mol for P450 and 310 kJ/mol for OPH for the dimermonomer dissociation. Accordingly, the enthalpy change (ΔH = ΔG + TΔS) is estimated to be -140 kJ/mol for P450 and -210 kJ/mol for OPH by ignoring other contributions (e.g., entropy of the protein). It could be suggested that for both proteins, the solvent entropy plays a determinant role in protein dimerization.

51

Moreover, the distributions of the SASA and the protein-protein distance in the whole

umbrella sampling process (Fig. S2) also support that the protein dimers are the dominant forms. “Mirror”-Like Protein Dimer Orientation Stabilized by Inter-Protein Hydrogen Bonds The previous section demonstrates that the native orientations of protein dimers are energetically favorable. Thus, a follow-up question is: what structural uniqueness is driving the protein dimer native orientation? To understand the mechanism that is driving the inter-protein orientation, we first calculated the interprotein potential energies. As listed in Table S1, the inter-protein interactions between polar amino acids outcompete the others, evidencing the significant role of hydrophilic interactions in protein dimerization. Moreover, it is shown that the Coulomb interactions dominate between polar amino acids, whereas the LJ interactions dominate when the nonpolar amino acids are involved. In these calculations, a cutoff distance of 1.2 nm was employed, which is smaller than the radius of the P450 monomer of around 2.5 nm and about 2.1 nm for the OPH monomer (Fig. 3). Therefore, these calculations represent the interprotein interactions on the contact regime and adjacent to the contact regime. We then focus on the contact regime exclusively where the two protein subunits are in direct contact with each other. A visualization supports a large amount of inter-protein H-bonds (Figure 4A). In the dimeric form, the P450 proteins are stabilized by around 13 ± 2 inter-protein H-bonds, and 8 ± 1 H-bonds for OPH. The number of H-bonds drops quickly with the increase of the protein-protein distance and disappears beyond around 6 nm. Here, the structural criteria of H-bonds are employed, that is, the donoracceptor distance rDA ≤ 0.35 nm, and the H−donor−acceptor angle θHDA ≤ 30°.54-55

ACS Paragon Plus Environment

9

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 20

Figure 4. (A) The number of inter-protein H-bonds as a function of the protein-protein distance. Snapshots of (B) P450 and (C) OPH dimers. Highlighted in the insets of (B, C) are the amino acids which are forming inter-protein H-bonds (yellow rods). The two subunits are colored in silver and cyan, with the O/N/H atoms colored in red/blue/white, respectively.

The energy barrier in breaking one amide H-bond has been estimated to be approximately 17 – 21 kJ/mol.51 Consequently, the energy barrier in breaking all the inter-protein H-bonds is around 221 – 273 kJ/mol for P450 and 136 – 168 kJ/mol for OPH. As demonstrated below, the inter-protein H-bonds play a crucial role in driving the relative orientations of the protein dimers. Based on the last snapshots in the protein dimer simulations, we obtained a list of the inter-protein Hbonds, which are illustrated in Figure 4B and listed in Table 1 for P450, and Figure 4C and Table 2 for OPH. It is found that the two subunits are roughly the “mirror” to each other. For instance, for P450 dimer Asp121 on subunit A forms an H-bond with Arg132 on subunit B; meanwhile, Arg132 on protein A forms an H-bond with Asp121 on protein B. Similar observations have been obtained for the OPH dimer. For example, Ser280 and Ser35 on subunit A form H-bonds with Ser110 on subunit B; meanwhile, Ser280 and Ser35 on protein B form H-bonds with Ser110 on protein A. The majority of the contacting amino acids (6 out of 8 for P450; 2 out of 3 for OPH) on protein A are shown to have their “mirror” amino acids on protein B located in the contact regime. Such “mirror”-like structures strongly demonstrate the uniqueness of the contact regime in protein dimerization. ACS Paragon Plus Environment

10

Page 11 of 20 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Also listed in Tables 1 and 2 are the types of the amino acids (P0: polar and neutral; Pn: polar and negatively charged; Pp: polar and positively charged) which are forming inter-protein H-bonds. The contacting amino acids are exclusively polar amino acids (neutral or charged). It is thus concluded that the local hydrophilic interactions between the neighbor proteins are playing a predominant role in the orientation of protein neighbors.

Table 1. H-bonds Between the Two Subunits of P450 Dimer in Figure 4B. Subunit A Subunit B number a type b

Asp 121 Arg 132 1 Pn - Pp

Arg Asn Arg Asp Gln Asp Gln Asn 132 134 161 168 169 168 125 163 Asp Arg Asn Gln Asp Asp Gln Arg 121 161 134 169 168 168 128 132 3 2 2 1 1 2 1 1 Pp - Pn P0 - Pp Pp - P0 Pn - P0 P0 - Pn Pn - Pn P0 - P0 P0 - Pp

a) Number of H-bonds b) Pn: polar and negatively charged; Pp: polar and positively charged; P0: polar and neutral

Table 2. H-bonds Between the Two Subunits of OPH Dimer in Figure 4C. Subunit A Subunit B number a type b

Ser Ser Ser 35 110 280 Ser Ser Ser 110 35 110 2 2 1 P0 - P0 P0 - P0 P0 - P0

Ser 110 Ser 280 1 P0 - P0

Lys 132 Glu 44 1 Pp - Pn

a) Number of H-bonds b) Pn: polar and negatively charged; Pp: polar and positively charged; P0: polar and neutral

In fact, it has been shown that proteins rarely work alone. Around 50% of the reported structures in the Protein Data Bank are protein complexes, most of which are composed of four subunits or less, and are homomeric.56 The homodimeric protein complexes can form symmetric57-59 or asymmetric60-61 structures. To explore if our findings above are universal, another protein dimer is explored, too. Remarkably, the “mirror”-like orientation is also observed for the XoxF dimer (Table S2). XoxF is an enzyme that efficiently catalyzes methanol oxidation in the exclusive presence of light lanthanide metal ions,62-63 and is therefore an important enzyme for applications in creating low-carbon chemicals. Condensing XoxF ACS Paragon Plus Environment

11

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 20

in aqueous solution provides a promising method to efficiently separate lanthanides, and importantly, an alternative means for recycling rare earth elements. 64-66 Local Heterogeneity at Protein Surface Drives Protein Dimer Orientation

Figure 5. Electrostatic potential at OPH dimer surface: (A) lower and (B) upper subunits. Negative and positive potentials are colored in red/blue, respectively. The snapshots are rotated to display the active sites (cave regions, blue ellipses) of the two subunits.

To understand the mechanism driving protein dimerization, it is a prerequisite to calibrate the local heterogeneity at protein surface. We first calculated the electrostatic potential at the surface of the OPH dimer. In these calculations, the lower and the upper subunits were saved separately, both of which were converted to the .pqr format using the PDB2PQR program.67 The recommended parameters were employed: the implicit solvent PARSE force field and the PROPKA program68 at pH = 7 for the ionization states of amino acids. The Adaptive Poisson–Boltzmann Solver (APBS)69 plugin of the package VMD 70 was then employed for the electrostatic potential calculation and visualization. The obtained results are demonstrated in Figure 5. First of all, it is demonstrated that the active sites are not buried by the neighbor protein. Such organization is essential for preserving their enzymatic activity when dimerized. This finding holds true for the P450 dimer. It may thus be feasible to synthesize enzyme dimer aggregates for advanced catalytic performance by elevating their concentrations. The negative (red) and positive (blue) electrostatic potentials coexist on both the lower and the upper subunits. Figure 5 also shows that the two subunits rotate against each other by some degree. The rotation ACS Paragon Plus Environment

12

Page 13 of 20 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

angle was calculated to be 103 ± 6 degree for the P450 dimer and 60 ± 3 for the OPH dimer. These stand for the angles between the normal vectors of two planes. The plane of each subunit was created using three coordinates: the metals in the active site (Fe in P450, one Zn in OPH), the geometric center of the subunit backbone atoms, and the corresponding center of the neighbor subunit. The error bars refer to the standard deviations calculated over the simulation time. The small error bars indicate that the orientations of the protein dimers are quite stable.

Figure 6. Distribution of surface domains for proteins P450 (top) and OPH (bottom). The distributions of the amino acids in the contact regimes (Tables 1 and 2) are labeled by yellow crosses.

To further understand the uniqueness of the proteins’ contact regime, we calculated the distributions of the polar and nonpolar protein surface domains (Figure 6). The spherical coordinates were employed. To display the protein-protein contact regime at the centers of the plots, the reference subunit was translated to the center of the simulations box, the dimer was rotated so that the neighbor subunit was aligned in the opposite of the X-dimension, as shown in the upper inset. The translation and rotation of the protein dimers were removed over the simulation duration using the first configuration. These calculations were done using an in-house script, which is based on our previous work

9-10

and was

improved to handle the lumpy surface of protein P450. See the Supporting Information for more details. It is demonstrated in Figure 6 that small polar surface domains coexist with small nonpolar surface domains at the surfaces of both proteins, similar to what we previously found for other proteins.9-10 Also presented in Figure 6 using yellow crosses are the locations of the amino acids which are distributed in the contact regime (Tables 1 and 2). These contacting amino acids are exclusively located on (charged ACS Paragon Plus Environment

13

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 20

or neutral) polar surface domains. They are distributed within the region of 50⁰ < θ < 130⁰ and 150⁰ < φ < 220⁰. Such broad distributions favor multiple inter-protein contacts, which cooperatively drive the specific orientations of protein dimers. Note that for some proteins the hydrophobic interactions are argued to be crucial in protein dimerization, where the contacting regions are occupied by the nonpolar amino acids.71 For instance, the aromatic side groups on phenylalanine, tyrosine and tryptophan amino acids are abundant in the contact region of STAT5 protein.72 Conclusions Within the context of quantifying protein aggregation, we calculated the free energy profile at the protein surface using protein-protein docking, and the potential of mean force of protein dimer-monomer dissociation using the steered MD and umbrella sampling approaches. The dimeric forms are found to be more energetically favorable for proteins P450 and OPH. Calculations of the solvent accessible surface areas supported that the solvent entropy contribution plays a determinant role in stabilizing the protein dimers. Moreover, the protein-protein docking results for the OPH dimer evidenced the uniqueness of the native orientation of protein dimers. The inter-protein H-bonds were shown to stabilize the relative orientation of the protein dimers. In particular, there exist around 10 inter-protein H-bonds for P450 and OPH dimers. Inspection of the contact regime demonstrated that the protein dimers are forming “mirror”-like structures with some degree of rotation. The protein surface heterogeneity is shown to play a crucial role in stabilizing the “mirror”-like protein dimer orientation. These findings support that the local heterogeneity of polar and nonpolar amino acids on the protein surface, though small in size (nm), is strong enough to drive the recognition of ligands. In our previous work, the local heterogeneity at a protein surface was shown to form stable protein-polymer complexes, which consequently preserved the enzymatic activity in non-native environments.9-10 Moreover, the unique fingerprint of protein-lipid interactions73 are likely related to the local heterogeneity at a protein surface. In the dimeric forms the enzyme active sites are exposed, suggesting that it is feasible to fabricate P450 and OPH dimer aggregates for elevated catalytic performance by increasing the concentrations of the enzymes, for example, by using suitable random heteropolymers

9

or by designing specific

polyampholytes. 74 The aggregates will be structurally similar to membraneless organelles. 8 That is, the P450 and OPH dimers, which are the active forms, will be employed as the building units to design the

ACS Paragon Plus Environment

14

Page 15 of 20 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

molecules. The molecules can connect neighbor building units, while preserving the stability and the enzymatic activity of P450 and OPH dimers in aqueous and organic conditions.

ASSOCIATED CONTENT Supporting Information The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/XXX. Supporting figures and tables, and details in calculating protein surface domains.

AUTHOR INFORMATION Corresponding Authors *E-mail: [email protected], Phone: +1-8474917801 ORCID Baofu Qiao: 0000-0001-8870-5985 Luis G. Lopez: 0000-0002-1587-6043 Notes The authors declare no competing financial interest. ACKNOWLEDGMENTS. The research was supported by Department of Energy Award DE-FG0208ER46539 and the Sherman Fairchild Foundation.

REFERENCES 1. Mikhaleva, S.; Lemke, E. A. Beyond the Transport Function of Import Receptors: What's All the FUS About? Cell 2018, 173, 549-553. 2. Dolgin, E. What Lava Lamps and Vinaigrette Can Teach Us About Cell Biology. Nature 2018, 555, 300-302. 3. Maharana, S.; Wang, J.; Papadopoulos, D. K.; Richter, D.; Pozniakovsky, A.; Poser, I.; Bickle, M.; Rizk, S.; Guillén-Boixet, J.; Franzmann, T., et al. RNA Buffers the Phase Separation Behavior of Prion Like RNA Binding Proteins. Science 2018, 360, 918-921.

ACS Paragon Plus Environment

15

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 20

4. Wang, Y.; Lomakin, A.; McManus, J. J.; Ogun, O.; Benedek, G. B. Phase Behavior of Mixtures of Human Lens Proteins Gamma D and Beta B1. Proc. Natl. Acad. Sci. U. S. A. 2010, 107, 1328213287. 5. Levine, Z. A.; Larini, L.; LaPointe, N. E.; Feinstein, S. C.; Shea, J.-E. Regulation and Aggregation of Intrinsically Disordered Peptides. Proc. Natl. Acad. Sci. U. S. A. 2015, 112, 2758-2763. 6. Falahati, H.; Wieschaus, E. Independent Active and Thermodynamic Processes Govern the Nucleolus Assembly in Vivo. Proc. Natl. Acad. Sci. U. S. A. 2017, 114, 1335-1340. 7. Boeynaems, S.; Alberti, S.; Fawzi, N. L.; Mittag, T.; Polymenidou, M.; Rousseau, F.; Schymkowitz, J.; Shorter, J.; Wolozin, B.; Van Den Bosch, L., et al. Protein Phase Separation: A New Phase in Cell Biology. Trends Cell Biol. 2018, 28, 420-435. 8. Nott, T. J.; Craggs, T. D.; Baldwin, A. J. Membraneless Organelles Can Melt Nucleic Acid Duplexes and Act as Biomolecular Filters. Nat. Chem. 2016, 8, 569–575. 9. Panganiban, B.; Qiao, B.; Jiang, T.; DelRe, C.; Obadia, M. M.; Nguyen, T. D.; Smith, A. A. A.; Hall, A.; Sit, I.; Crosby, M. G., et al. Random Heteropolymers Preserve Protein Function in Foreign Environments. Science 2018, 359, 1239-1243. 10. Nguyen, T. D.; Qiao, B.; Olvera de la Cruz, M. Efficient Encapsulation of Proteins with Random Copolymers. Proc. Natl. Acad. Sci. U. S. A. 2018, 115, 6578-6583. 11. Marianayagam, N. J.; Sunde, M.; Matthews, J. M. The Power of Two: Protein Dimerization in Biology. Trends Biochem. Sci. 2004, 29, 618-625. 12. Arnold, F. H. Directed Evolution: Bringing New Chemistry to Life. Angew. Chem. Int. Ed. 2018, 57, 4143-4148. 13. Simonian, A. L.; Good, T. A.; Wang, S. S.; Wild, J. R. Nanoparticle Based Optical Biosensors for the Direct Detection of Organophosphate Chemical Warfare Agents and Pesticides. Anal. Chim. Acta 2005, 534, 69-77. 14. Singh, B. K.; Walker, A. Microbial Degradation of Organophosphorus Compounds. FEMS Microbiol. Rev. 2006, 30, 428-471. 15. Neeli, R.; Girvan, H. M.; Lawrence, A.; Warren, M. J.; Leys, D.; Scrutton, N. S.; Munro, A. W. The Dimeric Form of Flavocytochrome P450 BM3 Is Catalytically Functional as a Fatty Acid Hydroxylase. FEBS Lett. 2005, 579, 5582-5588. 16. Grimsley, J. K.; Scholtz, J. M.; Pace, C. N.; Wild, J. R. Organophosphorus Hydrolase Is a Remarkably Stable Enzyme That Unfolds through a Homodimeric Intermediate. Biochemistry 1997, 36, 14366-14374. 17. Ortiz de Montellano, P. R. Hydrocarbon Hydroxylation by Cytochrome P450 Enzymes. Chem. Rev. 2010, 110, 932-948. 18. Munro, A. W.; Girvan, H. M.; Mason, A. E.; Dunford, A. J.; McLean, K. J. What Makes a P450 Tick? Trends Biochem. Sci. 2013, 38, 140-150. 19. Guengerich, F. P. Cytochrome P450 and Chemical Toxicology. Chem. Res. Toxicol. 2008, 21, 70-83. 20. Munro, A. W.; Leys, D. G.; McLean, K. J.; Marshall, K. R.; Ost, T. W. B.; Daff, S.; Miles, C. S.; Chapman, S. K.; Lysek, D. A.; Moser, C. C., et al. P450 BM3: The Very Model of a Modern Flavocytochrome. Trends Biochem. Sci. 2002, 27, 250-257. 21. Fasan, R.; Chen, M. M.; Crook, N. C.; Arnold, F. H. Engineered Alkane Hydroxylating Cytochrome P450BM3 Exhibiting Nativelike Catalytic Properties. Angew. Chem. Int. Ed. 2007, 46, 8414-8418. 22. Whitehouse, C. J. C.; Bell, S. G.; Wong, L.-L. P450BM3 (CYP102A1): Connecting the Dots. Chem. Soc. Rev. 2012, 41, 1218-1260. 23. Scott, E. E.; He, Y. A.; Wester, M. R.; White, M. A.; Chin, C. C.; Halpert, J. R.; Johnson, E. F.; Stout, C. D. An Open Conformation of Mammalian Cytochrome P450 2B4 at 1.6 Å Resolution. Proc. Natl. Acad. Sci. U. S. A. 2003, 100, 13196-13201. ACS Paragon Plus Environment

16

Page 17 of 20 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

24. Chenge, J. T.; Duyet, L. V.; Swami, S.; McLean, K. J.; Kavanagh, M. E.; Coyne, A. G.; Rigby, S. E. J.; Cheesman, M. R.; Girvan, H. M.; Levy, C. W., et al. Structural Characterization and Ligand/Inhibitor Identification Provide Functional Insights into the Mycobacterium Tuberculosis Cytochrome P450 CYP126A1. J. Biol. Chem. 2017, 292, 1310-1329. 25. Duffell, K. M.; Hudson, S. A.; McLean, K. J.; Munro, A. W.; Abell, C.; Matak-Vinković, D. Nanoelectrospray Ionization Mass Spectrometric Study of Mycobacterium Tuberculosis CYP121Ligand Interactions. Anal. Chem. 2013, 85, 5707-5714. 26. Black, S. D.; Martin, S. T. Evidence for Conformational Dynamics and Molecular Aggregation in Cytochrome P450 102 (BM3). Biochemistry 1994, 33, 12056-12062. 27. Zhang, H.; Yokom, A. L.; Cheng, S.; Su, M.; Hollenberg, P. F.; Southworth, D. R.; Osawa, Y. The Full Length Cytochrome P450 Enzyme CYP102A1 Dimerizes at Its Reductase Domains and Has Flexible Heme Domains for Efficient Catalysis. J. Biol. Chem. 2018, 293, 7727-7736. 28. Benning, M. M.; Kuo, J. M.; Raushel, F. M.; Holden, H. M. Three-Dimensional Structure of Phosphotriesterase: An Enzyme Capable of Detoxifying Organophosphate Nerve Agents. Biochemistry 1994, 33, 15001-15007. 29. Benning, M. M.; Shim, H.; Raushel, F. M.; Holden, H. M. High Resolution X-ray Structures of Different Metal Substituted Forms of Phosphotriesterase from Pseudomonas Diminuta. Biochemistry 2001, 40, 2712-2722. 30. Carletti, E.; Jacquamet, L.; Loiodice, M.; Rochu, D.; Masson, P.; Nachon, F. Update on Biochemical Properties of Recombinant Pseudomonas Diminuta Phosphotriesterase. J. Enzyme Inhib. Med. Chem. 2009, 24, 1045-1055. 31. Hess, B.; Kutzner, C.; van der Spoel, D.; Lindahl, E. Gromacs 4: Algorithms for Highly Efficient, Load Balanced, and Scalable Molecular Simulation. J. Chem. Theory Comput. 2008, 4, 435447. 32. Huang, J.; Rauscher, S.; Nawrocki, G.; Ran, T.; Feig, M.; de Groot, B. L.; Grubmuller, H.; MacKerell Jr, A. D. CHARMM36m: An Improved Force Field for Folded and Intrinsically Disordered Proteins. Nat. Meth. 2017, 14, 71-73. 33. MacKerell, A. D.; Bashford, D.; Bellott, M.; Dunbrack, R. L.; Evanseck, J. D.; Field, M. J.; Fischer, S.; Gao, J.; Guo, H.; Ha, S., et al. All-Atom Empirical Potential for Molecular Modeling and Dynamics Studies of Proteins. J. Phys. Chem. B 1998, 102, 3586-3616. 34. Miyamoto, S.; Kollman, P. A. Settle: An Analytical Version of the SHAKE and RATTLE Algorithm for Rigid Water Models. J. Comput. Chem. 1992, 13, 952-962. 35. Sevrioukova, I. F.; Li, H.; Zhang, H.; Peterson, J. A.; Poulos, T. L. Structure of a Cytochrome P450 Redox Partner Electron Transfer Complex. Proc. Natl. Acad. Sci. U. S. A. 1999, 96, 1863-1868. 36. Omburo, G. A.; Kuo, J. M.; Mullins, L. S.; Raushel, F. M. Characterization of the Zinc Binding Site of Bacterial Phosphotriesterase. J. Biol. Chem. 1992, 267, 13278-83. 37. Darden, T.; York, D.; Pedersen, L. Particle Mesh Ewald: An N·Log(N) Method for Ewald Sums in Large Systems. J. Chem. Phys. 1993, 98, 10089-10092. 38. Essmann, U.; Perera, L.; Berkowitz, M. L.; Darden, T.; Lee, H.; Pedersen, L. A Smooth Particle Mesh Ewald Method. J. Chem. Phys. 1995, 103, 8577-93. 39. Hess, B. P-LINCS: A Parallel Linear Constraint Solver for Molecular Simulation. J. Chem. Theory Comput. 2008, 4, 116-122. 40. Hess, B.; Bekker, H.; Berendsen, H. J. C.; Fraaije, J. G. E. M. LINCS: A Linear Constraint Solver for Molecular Simulations. J. Comput. Chem. 1997, 18, 1463-1472. 41. Shankar, K.; M., R. J.; Djamal, B.; H., S. R.; A., K. P. The Weighted Histogram Analysis Method for Free Energy Calculations on Biomolecules. I. The Method. J. Comput. Chem. 1992, 13, 1011-1021.

ACS Paragon Plus Environment

17

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 20

42. Hub, J. S.; de Groot, B. L.; van der Spoel, D. g_wham: A Free Weighted Histogram Analysis Implementation Including Robust Error and Autocorrelation Estimates. J. Chem. Theory Comput. 2010, 6, 3713-3720. 43. Ritchie, D. W.; Kozakov, D.; Vajda, S. Accelerating and Focusing Protein-Protein Docking Correlations Using Multi-Dimensional Rotational FFT Generating Functions. Bioinformatics 2008, 24, 1865-1873. 44. Camacho, C. J.; Zhang, C. FastContact: Rapid Estimate of Contact and Binding Free Energies. Bioinformatics 2005, 21, 2534-2536. 45. Rackers, J. A.; Wang, Z.; Lu, C.; Laury, M. L.; Lagardere, L.; Schnieders, M. J.; Piquemal, J.P.; Ren, P.; Ponder, J. W. Tinker 8: Software Tools for Molecular Design. J. Chem. Theory Comput. 2018, 14, 5273-5289. 46. Pickersgill, R. W. A Rapid Method of Calculating Charge-Charge Interaction Energies in Proteins. Protein Eng. 1988, 2, 247-248. 47. Zhang, C.; Vasmatzis, G.; Cornette, J. L.; DeLisi, C. Determination of Atomic Desolvation Energies from the Structures of Crystallized Proteins. J. Mol. Biol. 1997, 267, 707-726. 48. Strazza, S.; Hunter, R.; Walker, E.; Darnall, D. W. The Thermodynamics of Bovine and Porcine Insulin and Proinsulin Association Determined by Concentration Difference Spectroscopy. Arch. Biochem. Biophys. 1985, 238, 30-42. 49. Banerjee, P.; Mondal, S.; Bagchi, B. Insulin Dimer Dissociation in Aqueous Solution: A Computational Study of Free Energy Landscape and Evolving Microscopic Structure Along the Reaction Pathway. J. Chem. Phys. 2018, 149, 114902. 50. Zoete, V.; Meuwly, M.; Karplus, M. Study of the Insulin Dimerization: Binding Free Energy Calculations and Per-Residue Free Energy Decomposition. Proteins: Struct., Funct., Bioinf. 2005, 61, 79-93. 51. Chothia, C.; Janin, J. Principles of Protein-Protein Recognition. Nature 1975, 256, 705-708. 52. Rodriguez, R. A.; Yu, L.; Chen, L. Y. Computing Protein-Protein Association Affinity with Hybrid Steered Molecular Dynamics. J. Chem. Theory Comput. 2015, 11, 4427-4438. 53. Chothia, C. Hydrophobic Bonding and Accessible Surface Area in Proteins. Nature 1974, 248, 338-339. 54. Luzar, A.; Chandler, D. Structure and Hydrogen Bond Dynamics of Water-Dimethyl Sulfoxide Mixtures by Computer Simulations. J. Chem. Phys. 1993, 98, 8160-8173. 55. Qiao, B.; Demars, T.; Olvera de la Cruz, M.; Ellis, R. J. How Hydrogen Bonds Affect the Growth of Reverse Micelles around Coordinating Metal Ions. J. Phys. Chem. Lett. 2014, 5, 1440-1444. 56. Levy, E. D.; Pereira-Leal, J. B.; Chothia, C.; Teichmann, S. A. 3D Complex: A Structural Classification of Protein Complexes. PLoS Comput. Biol. 2006, 2, e155. 57. Goodsell, D. S.; Olson, A. J. Structural Symmetry and Protein Function. Annu. Rev. Biophys. Biomol. Struct. 2000, 29, 105-153. 58. André, I.; Strauss, C. E. M.; Kaplan, D. B.; Bradley, P.; Baker, D. Emergence of Symmetry in Homooligomeric Biological Assemblies. Proc. Natl. Acad. Sci. U. S. A. 2008, 105, 16148-16152. 59. Gross, M. Symmetry and Complexity in Protein Oligomers. Curr. Biol. 2012, 22, R175-R177. 60. Ye, Z.; Needham, P. G.; Estabrooks, S. K.; Whitaker, S. K.; Garcia, B. L.; Misra, S.; Brodsky, J. L.; Camacho, C. J. Symmetry Breaking During Homodimeric Assembly Activates an E3 Ubiquitin Ligase. Sci. Rep. 2017, 7, 1789. 61. Brown, J. H. Breaking Symmetry in Protein Dimers: Designs and Functions. Protein Sci. 2006, 15, 1-13. 62. Pol, A.; Barends, T. R. M.; Dietl, A.; Khadem, A. F.; Eygensteyn, J.; Jetten, M. S. M.; Op den Camp, H. J. M. Rare Earth Metals Are Essential for Methanotrophic Life in Volcanic Mudpots. Environ. Microbiol. 2014, 16, 255-264. 63. Skovran, E.; Martinez-Gomez, N. C. Just Add Lanthanides. Science 2015, 348, 862-863. ACS Paragon Plus Environment

18

Page 19 of 20 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

64. Qiao, B.; Skanthakumar, S.; Soderholm, L. Comparative CHARMM and AMOEBA Simulations of Lanthanide Hydration Energetics and Experimental Aqueous Solution Structures. J. Chem. Theory Comput. 2018, 14, 1781-1790. 65. Qiao, B.; Ferru, G.; Ellis, R. J. Complexation Enhancement Drives Water-to-Oil Ion Transport: A Simulation Study. Chem. Eur. J 2017, 23, 427-436. 66. Ferru, G.; Reinhart, B.; Bera, M. K.; Olvera de la Cruz, M.; Qiao, B.; Ellis, R. J. The Lanthanide Contraction Beyond Coordination Chemistry. Chem. Eur. J. 2016, 22, 6899-6904. 67. Dolinsky, T. J.; Nielsen, J. E.; McCammon, J. A.; Baker, N. A. PDB2PQR: An Automated Pipeline for the Setup of Poisson-Boltzmann Electrostatics Calculations. Nucleic Acids Res. 2004, 32, W665-W667. 68. Olsson, M. H. M.; Søndergaard, C. R.; Rostkowski, M.; Jensen, J. H. Propka3: Consistent Treatment of Internal and Surface Residues in Empirical Pka Predictions. J. Chem. Theory Comput. 2011, 7, 525-537. 69. Baker, N. A.; Sept, D.; Joseph, S.; Holst, M. J.; McCammon, J. A. Electrostatics of Nanosystems: Application to Microtubules and the Ribosome. Proc. Natl. Acad. Sci. U. S. A. 2001, 98, 10037-10041. 70. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual Molecular Dynamics. J. Mol. Graphics 1996, 14, 33-38. 71. Mei, G.; Di Venere, A.; Rosato, N.; Finazzi-Agrò, A. The Importance of Being Dimeric. FEBS J. 2005, 272, 16-27. 72. Fahrenkamp, D.; Li, J.; Ernst, S.; Schmitz-Van de Leur, H.; Chatain, N.; Küster, A.; Koschmieder, S.; Lüscher, B.; Rossetti, G.; Müller-Newen, G. Intramolecular Hydrophobic Interactions Are Critical Mediators of STAT5 Dimerization. Sci. Rep. 2016, 6, 35454. 73. Corradi, V.; Mendez-Villuendas, E.; Ingólfsson, H. I.; Gu, R.-X.; Siuda, I.; Melo, M. N.; Moussatova, A.; DeGagné, L. J.; Sejdiu, B. I.; Singh, G., et al. Lipid-Protein Interactions Are Unique Fingerprints for Membrane Proteins. ACS Cent. Sci. 2018, 4, 709-717. 74. Samanta, R.; Ganesan, V. Influence of Protein Charge Patches on the Structure of ProteinPolyelectrolyte Complexes. Soft Matter 2018, 14, 9475-9488.

ACS Paragon Plus Environment

19

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

TOC Graphic:

ACS Paragon Plus Environment

20

Page 20 of 20