Optimal ratio between the average conformational entropy and the average energy of interaction between residues for fast protein folding

Based on available experimental data and using a theoretical model of protein folding, we demonstrate that there is an optimal ratio between the average conformational entropy and the average contact energy per residue for fast protein folding. A statistical analysis of the conformational entropy and the number of contacts per residue for 5829 protein domains from four main classes (α, β, α/β, α+β) shows that each class has its own characteristic average number of contacts per residue and average conformational entropy per residue. These class-specific characteristics determine the protein folding rates: α-proteins are the fastest to fold, β-proteins are the second fastest, α+β-proteins are the third, and α/β-proteins are the last to fold.     

Surface free energy and interaction of Staphylococcus epidermidis with biomaterials

The adhesion of twenty nine Staphylococcus epidermidis strains to teflon, polyethylene, polycarbonate and bovine pericardium was studied in vitro and examined in relation to the surface free energies of both bacteria and biomaterials. All S. epidermidis strains had similar surface free energies, close to that of water, and adhered better to the materials with analogous surface free energies. There was a significant correlation (Kendall's Tau B = 1000) of biomaterial's surface free energy with the number of adhering bacteria. This correlation is inverse (Kendall's Tau B = -1000) when surface hydrophobicity is considered instead of surface free energy. This indicates that in Staphylococcus epidermidis adherence to biomaterials is inversely correlated to the surface hydrophobicity of the last, being so just the opposite of that occurring with other bacteria.     

Optimal relationship between average conformational entropy and average energy of interactions between residues for fast protein folding

Based on the known experimental data and using the theoretical modeling of protein folding, we demonstrate that there exists an optimal relationship between the average conformational entropy and the average energy of contacts per residue, that is an entropy capacity, for fast protein folding. Statistical analysis of conformational entropy and the number of contacts per residue for 5829 protein structures from four general structural classes (all-alpha, all-beta, +/-/beta, alpha+beta) demonstrates that each class of proteins has its own class-specific average number of contacts and average conformational entropy per residue. These class-specific features determine the folding rates: a proteins are the fastest folding proteins, then follow beta and alpha+beta proteins, and finally alpha/beta proteins are the slowest ones.     

Heat capacity-independent determination of differential free energy of stability between structurally homologous proteins

Under the assumption of equivalent heat capacity values, the differential free energy of stability for a pair of proteins midway between their thermal unfolding transition temperatures is shown to be independent of DeltaC(p) up to its cubic term in DeltaT(m). For model calculations reflecting the nearly 30 degrees C difference in T(m) for the adenylate kinases from the arctic bacterium Bacillus globisporus and the thermophilic bacterium Geobacillus stearothermophilus, the resultant error in estimating DeltaDeltaG by the formula 0.5 [DeltaS(T(m1))(1)+DeltaS(T(m2)) (2)] DeltaT(m) is less than 1%. Combined with the analogous thermal unfolding data for the adenylate kinase from Escherichia coli, these three homologous proteins exhibit T(m) and DeltaS(T(m)) values consistent with differential entropy and enthalpy contributions of equal magnitude. When entropy-enthalpy compensation holds for the differential free energy of stability, the incremental changes in T(m) values are shown to be proportionate to the changes in free energy.     

Measurement of histone-DNA interaction free energy in nucleosomes

Nucleosome positioning DNA sequences are of increasing interest because of their proposed roles in gene regulation and other chromosome functions in vivo, and because they have revealed new insights into the sequence-dependent structures and mechanics of DNA itself. Here, we describe methods to quantify the relative affinities of histone-DNA interactions in nucleosomes, i.e., the nucleosome positioning power of differing DNA sequences. We review methods developed by others and then discuss in detail our own approach to measurement of histone-DNA interaction free energies. Compared to earlier methods, our dialysis-based approach reduces the possibility that non-equilibrium or irreproducible results could be obtained. It facilitates a direct comparison of free energies for many sequences at the same time and it allows analysis of DNAs having a wide range of relative affinities.     

High free energy of lipid/protein interaction in biological membranes

The free energy of lipid/protein interaction in biological membranes is still unknown although extensive partitioning and modelling studies have revealed many partial energetic increments. Multiple site binding kinetics are now applied to four well-studied functional membrane proteins, and mean free energy values (+/-S.D.) of -4.23+/-0.49 kcal/mol for single lipid binding sites and of -89.7+/-35.4 kcal/mol for complete lipid substitution are obtained. These high free energy values point to an important bioenergetic role of lipid/protein interaction in membrane functions.     

Exploring the interaction between human focal adhesion kinase and inhibitors: a molecular dynamic simulation and free energy calculations

Focal adhesion kinase is an important target for the treatment of many kinds of cancers. Inhibitors of FAK are proposed to be the anticancer agents for multiple tumors. The interaction characteristic between FAK and its inhibitors is crucial to develop new inhibitors. In the present article, we used Molecular Dynamic (MD) simulation method to explore the characteristic of interaction between FAK and three inhibitors (PHM16, TAE226, and ligand3). The MD simulation results together with MM-GB/SA calculations show that the combinations are enthalpy-driven process. Cys502 and Asp564 are both essential residues due to the hydrogen bond interactions with inhibitors, which was in good agreement with experimental data. Glu500 can form a non-classical hydrogen bond with each inhibitor. Arg426 can form electrostatic interactions with PHM16 and ligand3, while weaker with TAE226. The electronic static potential was employed, and we found that the ortho-position methoxy of TAE226 has a weaker negative charge than the meta-position one in PHM16 or ligand3. Ile428, Val436, Ala452, Val484, Leu501, Glu505, Glu506, Leu553, Gly563 Leu567, Ser568 are all crucial residues in hydrophobic interactions. The key residues in this work will be available for further inhibitor design of FAK and also give assistance to further research of cancer.     

Insight into the dynamic interaction between different flavonoids and bovine serum albumin using molecular dynamics simulations and free energy calculations

In this study, the binding of Bovine serum albumin (BSA) with three flavonoids, kaempferol-3-O-a-L-rhamnopyranosyl-(1–3)-a-L-rhamnopyranosyl-(1–6)-b-D-galacto- pyranoside (drug 1),kaempfol-7-O-rhamnosyl-3-O-rutinoside (drug 2)andkaempferide-7-O-(4”-O-acetylrhamnosyl)-3-O-ruti- noside (drug 3) is investigated by molecular docking, molecular dynamics (MD) simulation, and binding free energy calculation. The free energies are consistent with available experimental results and suggest that the binding site of BSA-drug1 is more stable than those of BSA-drug2 and BSA-drug3. The energy decomposition analysis is performed and reveals that the electrostatic interactions play an important role in the stabilization of the binding site of BSA-drug1 while the van der Waals interactions contribute largely to stabilization of the binding site of BSA-drug2 and BSA-drug3. The key residues stabilizing the binding sites of BSA-drug1, BSA-drug2 and BSA-drug3 are identified based on the residue decomposition analysis.     

Relationship between physiological energy expenditure and biomechanical patterns associated with erect standing

The deviation of the static weight-bearing patterns under the feet of a lower extremity handicapped person may be measured quantitatively in terms of the Static Weight-Bearing Index. This paper describes the correlation between this index and the physiological energy expenditure associated with the erect standing posture. The Static Weight-Bearing Index can be conveniently used for evaluation of the functional status of the human lower extremity system in stance in case of lower extremity disabled persons who maintain basic weight-bearing mode involving both the limbs in double support condition.     

The exchange of energy between the medium and the active site

The idea of energy transduction realized by the enzyme was the objective of many studies during the last decades. This paper tries to analyze the possible effect of the transferral of energy from the environment to the active site through the enzyme-protein and the protein structure necessary to permit this kind of activity.     

Distances between active site probes in glutamine synthetase from Escherichia coli: fluorescence energy transfer in free and in stacked dodecamers

Probes for fluorescence energy transfer measurements were introduced into active sites of dodecameric glutamine synthetase from Escherichia coli by substituting appropriate ATP analogues for ATP in the autoinactivation reaction of this enzyme with L-Met-(S)-sulfoximine and Mn2+ [Maurizi, M. R., & Ginsburg, A. (1986) Biochemistry (preceding paper in this issue)]. Two fluorescent donors, 8-mercapto-ATP alkylated with either 5-[[[(iodoacetyl)amino]ethyl]amino]naphthalene-1-sulfonic acid (AEDANS-ATP) or 1,N6-etheno-2-aza-ATP (aza-epsilon-ATP), and two acceptors, 6-mercaptopurine ribonucleotide triphosphate or 8-mercapto-ATP alkylated with the chromophore 4'-[[4-(dimethylamino)-phenyl]azo]-2-iodoacetanilide (6-Y- or 8-Y-ATP), were used. Fluorescence emissions of enzyme derivatives with 1 or 2 equiv of fluorescent donor per dodecamer and either an acceptor (Y) or ADP at the remaining active sites were compared at pH 7.0. The results, together with the known geometry of the enzyme, indicate that active site probes in the dodecamer are widely separated and that energy transfer occurs from a single donor to two or three acceptors on adjacent subunits. The calculated distance between equidistant active site probes on heterologously bonded subunits within the same hexagonal ring is 56-60 A. Probes on isologously bonded subunits are no closer than 60 A and may be as far apart as 78 A. Thus, active sites are away from the 6-fold axis of symmetry toward the outer edges of the dodecamer and are located greater than or equal to 30 A from the plane separating the hexagonal rings. During Zn2+-induced stacking of the same enzyme derivatives along the 6-fold axes of symmetry, additional quenches of fluorescent probes were dependent on the presence of acceptors on separate dodecamers. The Zn2+-induced face to face aggregation of dodecamers in the presence of 46 microM ZnCl2 and 9 mM MgCl2 at pH 7.0 had an Arrhenius activation energy of 22.3 +/- 0.2 kcal/mol and a second-order rate constant at 25 degrees C of approximately 10(5) M-1 s-1 at early stages. Time-dependent fluorescence quenches correlated well with the degree of linear polymer formation and reached maximum values of 47-70% quench when the average n-mer was six dodecamers. After correction for unquenched polymer ends, a fluorescent donor and an acceptor probe in layered dodecamers were estimated to be approximately 36 A apart--an average value if there is some twisting of single strands.(ABSTRACT TRUNCATED AT 400 WORDS)     

A measure of DNA sequence dissimilarity based on free energy of nearest-neighbor interaction

We develop a novel method of asserting the similarity between two biological sequences without the need for alignment. The proposed method uses free energy of nearest-neighbor interactions as a simple measure of dissimilarity. It is used to perform a search for similarities of a query sequence against three complex datasets. The sensitivity and selectivity are computed and evaluated and the performance of the proposed distance measure is compared. Real data analysis shows that is a very efficient, sensitive and high-selective algorithm in comparing large dataset of DNA sequences.     

Targeting of the protein interaction site between FAK and IGF-1R

The interaction of focal adhesion kinase (FAK) and insulin-like growth factor-1 receptor (IGF-1R) plays an important role in cancer cell survival. Targeting this interaction with small molecule drugs could be a novel strategy in cancer therapy. By a series of pull-down assays using GST-tagged FAK fragments and His-tagged IGF-1R intracellular fragments, we showed that the FAK-NT2 (a.a. 127-243) domain directly interacts with the N-terminal part of the IGF-1R intracellular domain. Overexpressed FAK-NT2 domain was also shown to co-localize with IGF-1R in pancreatic cells. Computational modeling was used to predict the binding configuration of these two domains and to screen for small molecules binding to the interaction site. This strategy successfully identified a lead compound that disrupts FAK/IGF-1R interaction.     

Lipid bilayer deformation and the free energy of interaction of a Kv channel gating-modifier toxin

A number of membrane proteins act via binding at the water/lipid bilayer interface. An important example of such proteins is provided by the gating-modifier toxins that act on voltage-gated potassium (Kv) channels. They are thought to partition to the headgroup region of lipid bilayers, and so provide a good system for probing the nature of interactions of a protein with the water/bilayer interface. We used coarse-grained molecular dynamics simulations to compute the one-dimensional potential of mean force (i.e., free energy) profile that governs the interaction between a Kv channel gating-modifier toxin (VSTx1) and model phospholipid bilayers. The reaction coordinate sampled corresponds to the position of the toxin along the bilayer normal. The course-grained representation of the protein and lipids enabled us to explore extended time periods, revealing aspects of toxin/bilayer dynamics and energetics that would be difficult to observe on the timescales currently afforded by atomistic molecular dynamics simulations. In particular, we show for this model system that the bilayer deforms as it interacts with the toxin, and that such deformations perturb the free energy profile. Bilayer deformation therefore adds an additional layer of complexity to be addressed in investigations of membrane/protein systems. In particular, one should allow for local deformations that may arise due to the spatial array of charged and hydrophobic elements of an interfacially located membrane protein.     

Endothelial cells are a potential site of interaction between estrogens and glucocorticoids

Estrogens and glucocorticoids have synergistic effects in the micro and macrovasculature of endothelial cells. Whereas in the former they both have pro-inflammatory effects, in the latter they both inhibit the expression of adhesion molecules. Since the molecular basis of this synergy has not been defined, we propose possible mechanisms of interaction. An understanding of the functional interaction between estrogens and glucocorticoids in the endothelial cells of the micro and macrovasculature may contribute to clarifying the role of the endothelial cells of different vascular beds during the inflammatory response and in chronic inflammation. These advances could contribute to the design of more effective therapeutic strategies for the prevention of diseases, such as atherosclerosis, systemic lupus erythematosus and rheumatoid arthritis, in which the inflammatory process plays an important pathogenic role.