A heterodimerizing leucine zipper coiled coil system for examining the specificity of a position interactions: amino acids I,V, L,N, A,and K

We use a heterodimerizing leucine zipper system to examine the contribution of the interhelical a-a' interaction to dimer stability for six amino acids (A, V, L, I, K, and N). Circular dichroism (CD) spectroscopy monitored the thermal denaturation of 36 heterodimers that generate six homotypic and 30 heterotypic a-a' interactions. Isoleucine (I-I) is the most stable homotypic a-a' interaction, being 9.2 kcal/mol per dimer more stable than the A-A interaction and 4.0 kcal/mol per dimer more stable than either the L-L or V-V interaction, and 7.0 kcal/mol per dimer more stable than the N-N interaction. Only lysine was less stable than alanine. An alanine-based double-mutant thermodynamic cycle calculated coupling energies between the a and a' positions in the heterodimer. The aliphatic amino acids L, V, and I prefer to form homotypic interactions with coupling energies of -0.6 to -0.9 kcal/mol per dimer, but the heterotypic aliphatic interactions have positive coupling energies of <1.0 kcal/mol per dimer. The asparagine homotypic interaction has a coupling energy of -0.5 kcal/mol per dimer, while heterotypic interactions with the aliphatic amino acids produce coupling energies ranging from 2.6 to 4.9 kcal/mol per dimer. The homotypic K-K interaction is 2.9 kcal/mol per dimer less stable than the A-A interaction, but the coupling energy is only 0.3 kcal/mol per dimer. Heterotypic interactions with lysine and either asparagine or aliphatic amino acids produce similar coupling energies ranging from -0.2 to -0.7 kcal/mol per dimer. Thus, of the amino acids that were examined, asparagine contributes the most to dimerization specificity because of the large positive coupling energies in heterotypic interactions with the aliphatic amino acids which results in the N-N homotypic interaction.     

Thermodynamics of amide hydrogen bond formation in polar and apolar solvents

We present the initial findings of a theoretical study of hydrogen bond formation between two formamide molecules in water and in carbon tetrachloride. These systems were chosen as the simplest models for secondary structure formation in the polar environment near the protein surface and the apolar environment of the protein interior. We have employed thermodynamic simulation methods to obtain absolute binding free energies and free energy profiles for the formation of peptide hydrogen bonds in the two solvents. We find that the amide hydrogen bond is stable by 8.4 kcal/mol in CCl4, and by 0.3 kcal/mol in water. Our results indicate also that the hydrogen-bonded dimer is 2.2 kcal/mol more stable in water than it is in CCl4. We compare our results with those from experiment, and discuss their use in interpreting mechanisms of protein folding.     

Stability of a model beta-sheet in water.

We have used molecular dynamics simulations to determine the stability in water of a model beta-sheet formed by two alanine dipeptide molecules with two intermolecular hydrogen bonds in the closely spaced antiparallel arrangement. In this paper we describe our computations of the binding free energy of the model sheet and a portion of the free energy surface as a function of a reaction co-ordinate for sheet formation. We used the free energy surface to identify stable conformations along the reaction co-ordinate. To determine whether or not the model sheet with two hydrogen bonds is more stable than a single amide hydrogen bond in water, we compared the results of the present calculations to results from our earlier study of linear hydrogen bond formation between two formamide molecules (the formamide "dimer"). The free energy surfaces for the sheet and formamide dimer each have two minima corresponding to locally stable hydrogen-bonded and solvent-separated configurations. The binding free energies of the model sheet and the formamide dimer are -5.5 and -0.34 kcal/mol, respectively. Thus, the model sheet with two hydrogen bonds is quite stable while the simple amide hydrogen bond is only marginally stable. To understand the relative stabilities of the model sheet and formamide dimer in terms of solute-solute and solute-water interactions, we decomposed the free energy differences between hydrogen-bonded and solvent-separated conformations into energetic and entropic contributions. The changes in the peptide-peptide energy and the entropy are roughly twice as large for the sheet as they are for the formamide dimer. The magnitude of the peptide-water energy difference for the sheet is less than twice (by about 3.5 kcal/mol) that for the formamide dimer, and this accounts for the stability of the sheet. The presence of the side-chains and/or blocking groups apparently prevents the amide groups in the sheet from being solvated as favorably in the separated arrangement as in the formamide dimer, where the amide groups are completely exposed to the solvent.     

Molecular modeling reveals binding interface of γ‐tubulin with GCP4 and interactions with noscapinoids

The initiation of microtubule assembly within cells is guided by a cone shaped multi‐protein complex, γ‐tubulin ring complex (γTuRC) containing γ‐tubulin and atleast five other γ‐tubulin‐complex proteins (GCPs), i.e., GCP2, GCP3, GCP4, GCP5, and GCP6. The rim of γTuRC is a ring of γ‐tubulin molecules that interacts, via one of its longitudinal interfaces, with GCP2, GCP3, or GCP4 and, via other interface, with α/β?tubulin dimers recruited for the microtubule lattice formation. These interactions however, are not well understood in the absence of crystal structure of functional reconstitution of γTuRC subunits. In this study, we elucidate the atomic interactions between γ‐tubulin and GCP4 through computational techniques. We simulated two complexes of γ‐tubulin‐GCP4 complex (we called dimer1 and dimer2) for 25 ns to obtain a stable complex and calculated the ensemble average of binding free energies of ?158.82 and ?170.19 kcal/mol for dimer1 and ?79.53 and ?101.50 kcal/mol for dimer2 using MM‐PBSA and MM‐GBSA methods, respectively. These highly favourable binding free energy values points to very robust interactions between GCP4 and γ‐tubulin. From the results of the free‐energy decomposition and the computational alanine scanning calculation, we identified the amino acids crucial for the interaction of γ‐tubulin with GCP4, called hotspots. Furthermore, in the endeavour to identify chemical leads that might interact at the interface of γ‐tubulin‐GCP4 complex; we found a class of compounds based on the plant alkaloid, noscapine that binds with high affinity in a cavity close to γ‐tubulin‐GCP4 interface compared with previously reported compounds. All noscapinoids displayed stable interaction throughout the simulation, however, most robust interaction was observed for bromo‐noscapine followed by noscapine and amino‐noscapine. This offers a novel chemical scaffold for γ‐tubulin binding drugs near γ‐tubulin‐GCP4 interface. Proteins 2015; 83:827–843. © 2015 Wiley Periodicals, Inc.     

Strong interaction between disulfide derivatives and aromatic groups in peptides and proteins

The intermolecular interaction energy of complexes of dimethyldisulfide with benzene and cyclohexane, respectively, was computed as function of the relative distance and orientation within each pair of molecules. The energy of the most stable orientation of the dimethyldisulfide-cyclohexane complex is ?2.57 kcal/mol, while that of the most stable orientation of the dimethyldisulfide-benzene complex is ?3.33 kcal/mol. The energy difference of ～0.8 kcal/mol is due to favorable specific nonbonded interactions between the sulfur atoms and the atoms of the aromatic ring. Proper parameterization of empirical interatomic energies, used in computations in this laboratory, accounts for these interactions without the need for a $\text{special}$ sulfur-aromatic potential energy function.     

Structural basis of pesticide detection by enzymatic biosensing: a molecular docking and MD simulation study

Designing of rapid, facile, selective, and cost-effective biosensor technology is a growing area for the detection of various classes of pesticides. The biosensor with these features can be achieved only through the various bio-components using different transducers. This study, therefore, focuses on the usage of molecular docking, specificity tendencies, and capabilities of proteins for the detection of pesticides. Accordingly, the four transducers, acetylcholinesterase (ACH), cytochromes P450 (CYP), glutathione S-transferase (GST), and protein kinase C (PKC) were selected based on their applications including neurotransmitter, metabolism, detoxification enzyme, and protein phosphorylation. Then after molecular docking of the pesticides, fenobucarb, dichlorodiphenyltrichloroethane (DDT), and parathion onto each enzyme, the conformational behavior of the most stable complexes was further analyzed using 50 ns Molecular Dynamics (MD) simulations carried out under explicit water conditions. In the case of protein kinase C (PKC) and cytochrome P450 3A4 enzyme (CYP), the fenobucarb complex showed the most suitable combination of free energy of binding and inhibition constant ?4.42 kcal/mol (573.73 μM) and ?5.1 kcal/mol (183.49 μM), respectively. Parathion dominated for acetylcholinesterase (ACH) with ?4.57 kcal/mol (448.09 μM) and lastly dichlorodiphenyltrichloroethane for glutathione S-transferase (GST), ?5.43 kcal/mol (103.88 μM). The RMSD variations were critical for understanding the impact of pesticides as they distinctively influence the energetic attributes of the proteins. Overall, the outcomes from the extensive analysis provide an insight into the structural features of the proteins studied, thereby highlighting their potential use as a substrate in biorecognition sensing of pesticide compounds.     

Molecular dynamics study of free energy profiles for organic cations in gramicidin A channels.

The free energy profiles for four organic cations in right-handed single-helix gramicidin A dimers were computed by using umbrella sampling molecular dynamics with CHARMM. Ion-water column translocations were facilitated by using a novel "water-tunnel" approach. The overlapping pieces of free energy profile for adjacent windows were selected from three trajectories that differed in initial ion rotation and were aligned by the method of umbrella potential differences. Neglected long-range electrostatic energies from the bulk water and the bilayer were computed with DelPhi and added to the profile. The approach was corroborated for the formamidinium-guanidinium pair by using perturbation dynamics at axial positions 0, 6, 12, and 15 A from the channel center. The barrier to ethylammonium entry was prohibitive at 21 kcal/mol, whereas for methylammonium it was 5.5 kcal/mol, and the profile was quite flat through the channel, roughly consistent with conductance measurements. The profile for formamidinium was very similar to that of methylammonium. Guanidinium had a high entry barrier (deltaF = +8.6 kcal/mol) and a narrow deep central well (deltaF = -2.6 kcal/mol), qualitatively consistent with predictions from voltage-dependent potassium current block measurements. Its deep central well, contrasting with the flat profile for formamidinium, was verified with perturbation dynamics and was correlated with its high propensity to form hydrogen bonds with the channel at the dimer junction (not shared by the other three cations). Analysis of the ensemble average radial forces on the ions demonstrates that all four ions undergo compressive forces in the channel that are at maximum at the center of the monomer and relieved at the dimer junction, illustrating increased flexibility of the channel walls in the center of the channel.     

Ab Initio Quantum Mechanics Analysis of Imidazole C-H…O Water Hydrogen Bonding and a Molecular Mechanics Forcefield Correction

Abstract

While it is well established that classical hydrogen bonds play an important role in enzyme structure, function and dynamics, the role of weaker, but ‘activated’ C-H donor hydrogen bonds is poorly understood. The most important such case involves histidine which often plays a direct role in enzyme catalysis and possesses the most acidic C-H donor group of the standard amino acids. In the present study, we obtained optimized geometries and hydrogen bond interaction energies for C-H…O hydrogen bonded complexes between methane, ethylene, benzene, acetylene, and imidazole with water at the MP2-FC/6-31++G(2d,2p) and MP2-FC/aug-cc-pVDZ//MP2-FC/6-31++G(2d,2p) levels of theory. A strong linear relationship is obtained between the stability of the various hydrogen bonded complexes and both separation distances for H…0 and C—O. In general, these calculations indicate that C-H…0 interactions can be classified as hydrogen bonding interactions, albeit significantly weaker than the classical hydrogen bonds, but significantly stronger than just van der Waals interactions. For instance, while the electronic energy of stabilization at the MP2-FC/aug-cc-pVDZ//MP2-FC/6-31++G(2d,2p) level of theory of a water C-H…O water hydrogen bond is 4.36 kcal/mol more stable than the methane C-H…O water interaction, the water-water hydrogen bond is only 2.06 kcal/mol more stable than the imidazole C_e?H…O water hydrogen bond. Neglecting this latter hydrogen bonding interaction is obviously unacceptable. We next compare the potential energy surfaces for the imidazole C_e?H…O water and imidazole N_d?H…O hydrogen bonded complexes computed at the MP2/6-31++G(2d,2p) level of theory with the potential energy surface computed using the AMBER molecular mechanics program and forcefields. While the Weiner et al and Cornell et al AMBER forcefields reasonably account for the imidazole N-H…O water interaction, these forcefields do not adequately account for the imidazole C_e?H…O water hydrogen bond. A forcefield modification is offered that results in excellent agreement between the ab initio and molecular mechanics geometry and energy for this C-H…O hydrogen bonded complex.     

Chrysophanol anthrone: a theoretical study on the potential energy surface

The potential energy surface (PES) of chrysophanol anthrone, the active component of Goa Powder, was systematically explored and thoroughly scrutinised via density functional theory, in order to gain an understanding of its physicochemical properties. In particular, we focused on the rotations of the two hydroxyl-phenyl dihedral angles. A picture with a four stable rotamers emerged where only the most stable conformer has a planar structure and the less stable conformer has the maximum deviation from planarity. The computed PES shows that the energy barriers for the conformer interconversion are less than 15 kcal/mol. From the analysis of the calculated intramolecular hydrogen bond enthalpy, we conclude that the number of the intramolecular hydrogen bonds governs the conformer stability. Additionally, the conformational equilibrium was pursued by means of an analysis of the energy of OH internal rotation barriers. The total energy changes were decomposed in an electrostatic decomposition scheme in order to gain an insight into the effects governing the torsional barrier and preferred conformations. This analysis shows that the interplay between the repulsive and attractive potentials causes the conformer stability, where the attractive term dominates the conformer stabilisation.     

Human platelet factor 4 subunit association/dissociation thermodynamics and kinetics

Platelet factor 4 (PF4) monomers (7800 daltons) form dimers and tetramers in varying molar ratios under certain solution conditions [Mayo, K. H., & Chen, M. J. (1989) Biochemistry 28, 9469]. The presence of a simplified aromatic region (one Tyr and two His) and resolved monomer, dimer, and tetramer Y60 3,5 ring proton resonances makes study of PF4 aggregate association/dissociation thermodynamics and kinetics possible. PF4 protein subunit association/dissociation equilibrium thermodynamic parameters have been derived by 1H NMR (500MHz) resonance line-fitting analysis of steady-state Y60 3,5 ring proton resonance monomer-dimer-tetramer populations as a function of temperature from 10 to 40 degrees C. Below 10 degrees C and above 40 degrees C, resonance broadening and overlap severely impaired analysis. Enthalpic and entropic contributions to dimer association Gibb's free energy [-5.1 kcal/mol (30 degrees C)] are +2.5 +/- 1 kcal/mol and +26 +/- 7 eu, respectively, and for tetramer association Gibb's free energy [-5.7 kcal/mol (30 degrees C)], they are -7.5 +/- 1 kcal/mol and -7 +/- 3 eu, respectively. These thermodynamic parameters are consistent with low dielectric medium electrostatic/hydrophobic interactions governing dimer formation and hydrogen bonding governing tetramer formation. Association/dissociation kinetic parameters, i.e., steady-state jump rates, have been derived from exchange-induced line-width increases and from 1H NMR (500 MHz) saturation-transfer and spin-lattice (Tl) relaxation experiments. From dissociation jump rates and equilibrium constants, association rate constants were estimated. For dimer and tetramer equilibria at 30 degrees C, unimolecular dissociation rate constants are 35 +/- 10 s-1 for dimer dissociation and 6 +/- 2 s-1 for tetramer dissociation.(ABSTRACT TRUNCATED AT 250 WORDS)     

Helix-capping interaction in λ cro protein: A free energy simulation analysis

The stability mutant Tyr-26 → Asp was studied in the Cro protein from bacteriophage λ using free energy molecular dynamics simulations. The mutant was calculated to be more stable than the wild type by 3.0 ± 1.7 kcal/mol/monomer, in reasonable agreement with experiment (1.4 kcal/mol/monomer). Moreover, the aspartic acid in the mutant was found to form a capping interation with the amino terminus of the third α-helix of Cro. The simulations were analyzed to understand better the source of the stability of this helix-capping interaction and to examine the results in light of previous explanations of stabilizing helix caps-namely, a model of local unsatisfied hydrogen bonds at the helix termini and the helix macro dipole model. Analysis of the simulations shows that the stabilizing effect of this charged helical cap is due both to favorable hydrogen bonds with backbone NH groups at the helix terminus and to favorable electrostatic interactions (but not hydrogen bonds) with their carbonyls (effectively the next row of local dipoles in the helix). However, electrostatic interactions are weak or negligible with backbone dipolar groups in the helix further away from the terminus. Moreover, the importance of other local electrostatic interactions with polar side chains near the helix terminus, which are neglected in most treatments of this effect, are shown to be important. Thus, the results support a model that is intermediate between the two previous explanations: both unsatisfied hydrogen bonds at the helix terminus and other, local preoriented dipolar groups stabilize the helix cap. These findings suggest that similar interactions with preoriented dipolar groups may be important for cooperativity in other charge–dipole interactions and may be employed to advantage for molecular design. © 1994 Wiley-Liss, Inc.     

Computational study on C−H…π interactions of acetylene with benzene, 1,3,5-trifluorobenzene and coronene

Meta-hybrid density functional theory calculations using M06-2X/6-31+G(d,p) and M06-2X/6-311+G(d,p) levels of theory have been performed to understand the strength of C?H^…π interactions of two possible types for benzene-acetylene, 1,3,5-trifluorobenzene-acetylene and coronene-acetylene complexes. Our study reveals that the C?H^...π interaction complex where acetylene located above to the center of benzene ring (classical T-shaped) is the lowest energy structure. This structure is twice more stable than the configuration characterized by H atom of benzene interacting with the π-cloud of acetylene. The binding energy of 2.91 kcal/mol calculated at the M06-2X/6-311+G(d,p) level for the lowest energy configuration (1A) is in very good agreement with the experimental binding energy of 2.7?±?0.2 kcal/mol for benzene-acetylene complex. Interestingly, the C?H^...π interaction of acetylene above to the center of the aromatic ring is not the lowest energy configuration for 1,3,5-trifluorobenzene-acetylene and coronene-acetylene complexes. The lowest energy configuration (2A) for the former complex possesses both C?H^...π interaction and C?H^...F hydrogen bond, while the lowest energy structure for the coronene-acetylene complex involves both π-π and C?H^...π interactions. C?H stretching vibrational frequencies and the frequency shifts are reported and analyzed for all of the configurations. We observed red-shift of the vibrational frequency for the stretching mode of the C-H bond that interacts with the π-cloud. Acetylene in the lowest-energy structures of the complexes exhibits significant red-shift of the C?H stretching frequency and change in intensity of the corresponding vibrational frequency, compared to bare acetylene. We have examined the molecular electrostatic potential on the surfaces of benzene, 1,3,5-trifluorobenzene, coronene and acetylene to explain the binding strengths of various complexes studied here.     

Temperature dependence of kinetic interactions between progesterone and the uterine cytoplasmic receptor

The temperature dependence of the rates of dissociation and association for progesterone-receptor interactions was measured over the temperature range of 0–20°C. The dissociation process is biphasic indicating that either two forms of receptor are present or that the binding of progesterone to the receptor is a concatenated reaction.The enthalpy of activation for the dissociation of progesterone from the receptor is about 26–28 kcal/mol and the entropic energy of activation is about ?5 kcal/mol. The enthalpy of activation for the association of these molecules is about 3 kcal/mol and the entropic energy of activation is about 6 kcal/mol. These data are consistent with a model of progesterone binding to the receptor that includes hydrogen bonds between each of the two ketone groups and hydrogen donors on the receptor protein and involves van der Waals' interactions, due to the close proximity of the receptor binding site to a large fraction of the progesterone surface.     

Equilibrium folding of the core histones: the H3-H4 tetramer is less stable than the H2A-H2B dimer

To compare the stability of structurally related dimers and to aid in understanding the thermodynamics of nucleosome assembly, the equilibrium stabilities of the recombinant wild-type H3-H4 tetramer and H2A-H2B dimer have been determined by guanidinium-induced denaturation, using fluorescence and circular dichroism spectroscopies. The unfolding of the tetramer and dimer are highly reversible. The unfolding of the H2A-H2B dimer is a two-state process, with no detected equilibrium intermediates. The H3-H4 tetramer is unstable at moderate ionic strengths (mu approximately 0.2 M). TMAO (trimethylamine-N-oxide) was used to stabilize the tetramer; the stability of the H2A-H2B dimer was determined under the same solvent conditions. The equilibrium unfolding of H3-H4 was best described by a three-state mechanism, with well-folded H3-H4 dimers as a populated intermediate. When compared to H2A-H2B, the H3-H3 tetramer interface and the H3-H4 histone fold are strikingly less stable. The free energy of unfolding, in the absence of denaturant, for the H3-H4 and H2A-H2B dimers are 12.4 and 21.0 kcal mol(-)(1), respectively, in 1 M TMAO. It is postulated that the difference in stability between the histone dimers, which contain the same fold, is the result of unfavorable tertiary interactions, most likely the partial to complete burial of three salt bridges and burial of a charged hydrogen bond. Given the conservation of these buried interactions in histones from yeast to mammals, it is speculated that the H3-H4 tetramer has evolved to be unstable, and this instability may relate to its role in nucleosome dynamics.     

Thermodynamic and base-pairing studies of matched and mismatched DNA dodecamer duplexes containing cis-syn, (6-4) and Dewar photoproducts of TT.

Cis-syn dimers, (6-4) products and their Dewar valence isomers are the major photoproducts of DNA and have different mutagenic properties and rates of repair. To begin to understand the physical basis for these differences, the thermal stability and base pairing properties of the corresponding photoproducts of the TT site in d(GAGTATTATGAG) were investigated. The (6-4) and Dewar products destabilize the duplex form by approximately 6 kcal/mol of free energy at 37 degreesC relative to the parent, whereas a cis-syn dimer only destabilizes the duplex form by 1.5 kcal/mol. Duplexes with G opposite the 3'-T of the (6-4) and Dewar products are more stable than those with A by approximately 0.4 kcal/mol, whereas the cis-syn dimer prefers A over G by 0.7 kcal/mol. Proton NMR suggests that wobble base pairing takes place between the 3'-T of the cis-syn dimer and an opposed G, whereas there is no evidence of significant H-bonding between these two bases in the (6-4) product. The thermodynamic and H-bonding data for the (6-4) product are consistent with a 4 nt interior loop structure which may facilitate flipping of the photoproduct in and out of the helix.     

Membrane-Mediated Protein-Protein Interactions and Connection to Elastic Models: A Coarse-Grained Simulation Analysis of Gramicidin A Association

To further foster the connection between particle based and continuum mechanics models for membrane mediated biological processes, we carried out coarse-grained (CG) simulations of gramicidin A (gA) dimer association and analyzed the results based on the combination of potential of mean force (PMF) and stress field calculations. Similar to previous studies, we observe that the association of gA dimers depends critically on the degree of hydrophobic mismatch, with the estimated binding free energy of >10 kcal/mol in a distearoylphosphatidylcholine bilayer. Qualitative trends in the computed PMF can be understood based on the stress field distributions near a single gA dimer and between a pair of gA dimers. For example, the small PMF barrier, which is ∼1 kcal/mol independent of lipid type, can be captured nearly quantitatively by considering membrane deformation energy associated with the region confined by two gA dimers. However, the PMF well depth is reproduced poorly by a simple continuum model that only considers membrane deformation energy beyond the annular lipids. Analysis of lipid orientation, configuration entropy, and stress distribution suggests that the annular lipids make a significant contribution to the association of two gA dimers. These results highlight the importance of explicitly considering contributions from annular lipids when constructing approximate models to study processes that involve a significant reorganization of lipids near proteins, such as protein-protein association and protein insertion into biomembranes. Finally, large-scale CG simulations indicate that multiple gA dimers also form clusters, although the preferred topology depends on the protein concentration. Even at high protein concentrations, every gA dimer requires contact to lipid hydrocarbons to some degree, and at most three to four proteins are in contact with each gA dimer; this observation highlights another aspect of the importance of interactions between proteins and annular lipids.     

Investigation and optimisation of the drying of reduced glutathione by steered molecular dynamics simulation

The drying of reduced glutathione from a series of aqueous–ethanol binary solutions at 300 K (below human body temperature) and 330 K (above human body temperature) was investigated in detail by steered molecular simulation and an umbrella sampling method with the Gromacs software package and Gromos96(53a6) united atomic force field. The results show that electrostatic interactions between glutathione and solvent represent the main resistance to drying. When the aqueous solution was gradually changed to pure ethanol, the energy of electrostatic interaction between glutathione and solvent molecules increased by 445.088 kJ/mol, and the drying potential of mean force (PMF) free energy also fell by 253.040 kJ/mol. However, an increase in temperature from 300 to 330 K in the aqueous solution only results in an increase of 23.013 kJ/mol in electrostatic interaction energy and a decrease of 34.956 kJ/mol in drying PMF free energy. Furthermore, we show that hydrogen bonding is the major form of electrostatic interaction involved, and directly affects the drying of glutathione. Therefore, choosing water-miscible solvents that minimise hydrogen-bond formation with glutathione will enhance its drying rate, and this is likely to be more efficient than increasing the temperature of the process. Thus, a power-saving technology can be used to produce the high bioactivity medicines.     

Competition between halogen, dihalogen and hydrogen bonds in bromo- and iodomethanol dimers

O-H…X and O-H…O H-bonds as well as C-X…X dihalogen and C-X…O halogen bonds have been investigated in halomethanol dimers (bromomethanol dimer, iodomethanol dimer, difluorobromomethanol…bromomethanol complex and difluoroiodomethanol…iodomethanol complex). Structures of all complexes were optimized at the counterpoise-corrected MP2/cc-pVTZ level and single-point energies were calculated at the CCSD(T)/aug-cc-pVTZ level. Energy decomposition for the bromomethanol dimer complex was performed using the DFT-SAPT method based on the aug-cc-pVTZ basis set. OH…O and OH…X H-bonds are systematically the strongest in all complexes investigated, with the former being the strongest bond. Halogen and dihalogen bonds, being of comparable strength, are weaker than both H-bonds but are still significant. The strongest bonds were found in the difluoroiodomethanol…iodomethanol complex, where the O-H…O H-bond exceeds 7 kcal mol^-1, and the halogen and dihalogen bonds exceed 2.5 and 2.3 kcal mol^-1, respectively. Electrostatic energy is dominant for H-bonded structures, in halogen bonded structures electrostatic and dispersion energies are comparable, and, finally, for dihalogen structures the dispersion energy is clearly dominant.

Dissecting the Forces that Dominate Dimerization of the Nucleotide Binding Domains of ABCB1

P-glycoprotein, also known as multidrug resistance protein 1 or ABCB1, can export a wide range of chemically unrelated compounds, including chemotherapeutic drugs. ABCB1 consists of two transmembrane domains that form the substrate binding and translocation domain, and of two cytoplasmic nucleotide binding domains (NBDs) that energize substrate transport by ATP binding and hydrolysis. ATP binding triggers dimerization of the NBDs, which switches the transporter from an inward facing to an outward facing transmembrane domain conformation. We performed MD simulations to study the dynamic behavior of the NBD dimer in the presence or absence of nucleotides. In the apo configuration, the NBDs were overall attractive to each other as shown in the potential of mean force profile, but the energy well was shallow and broad. In contrast, a sharp and deep energy minimum (～?42 kJ/mol) was found in the presence of ATP, leading to a well-defined conformation. Motif interaction network analyses revealed that ATP stabilizes the NBD dimer by serving as the central hub for interdomain connections. Simulations showed that forces promoting dimerization are multilayered, dominated by electrostatic interactions between the nucleotide and conserved amino acids of the signature sequence and the Walker A motif. In addition, direct and water-bridged hydrogen bonds between NBDs provided conformation-defining interactions. Importantly, we characterized a largely unrecognized but essential contribution from hydrophobic interactions between the adenine moiety of the nucleotides and a hydrophobic surface of the X-loop to the stabilization of the nucleotide-bound NBD dimer. These hydrophobic interactions lead to a sharp energy minimum, thereby conformationally restricting the nucleotide-bound state.     

Conformational stability of dimeric proteins: quantitative studies by equilibrium denaturation.

The conformational stability of dimeric globular proteins can be measured by equilibrium denaturation studies in solvents such as guanidine hydrochloride or urea. Many dimeric proteins denature with a 2-state equilibrium transition, whereas others have stable intermediates in the process. For those proteins showing a single transition of native dimer to denatured monomer, the conformational stabilities, delta Gu (H2O), range from 10 to 27 kcal/mol, which is significantly greater than the conformational stability found for monomeric proteins. The relative contribution of quaternary interactions to the overall stability of the dimer can be estimated by comparing delta Gu (H2O) from equilibrium denaturation studies to the free energy associated with simple dissociation in the absence of denaturant. In many cases the large stabilization energy of dimers is primarily due to the intersubunit interactions and thus gives a rationale for the formation of oligomers. The magnitude of the conformational stability is related to the size of the polypeptide in the subunit and depends upon the type of structure in the subunit interface. The practical use, interpretation, and utility of estimation of conformational stability of dimers by equilibrium denaturation methods are discussed.