Novel poly(3-hydroxybutyrate) nanocomposites containing WS2 inorganic nanotubes with improved thermal, mechanical and tribological properties

Poly(3-hydroxybutyrate) (PHB) nanocomposites containing environmentally-friendly tungsten disulphide inorganic nanotubes (INTeWS2) have been successfully prepared by a simple solution blending method. The dynamic and isothermal crystallization studies by differential scanning calorimetry (DSC) demonstrated that the INTeWS2 exhibits much more prominent nucleation activity on the crystallization of PHB than specific nucleating agents or other nanoscale fillers. Both crystallization rate and crystallinity significantly increase in the nanocomposites compared to neat PHB. These changes occur without modifying the crystalline structure of PHB in the nanocomposites, as shown by wide-angle X-ray diffraction (WAXS) and infrared/Raman spectroscopy. Other parameters such as the Avrami exponent, the equilibrium melting temperature, global rate constant and the fold surface free energy of PHB chains in the nanocomposites were obtained from the calorimetric data in order to determine the influence of the INTeWS2 filler. The addition of INTeWS2 remarkably influences the energetics and kinetics of nucleation and growth of PHB, reducing the fold surface free energy by up to 20%. Furthermore, these nanocomposites also show an improvement in both tribological and mechanical (hardness and modulus) properties with respect to pure PHB evidenced by friction and nanoindentation tests, which is of important potential interest for industrial and medical applications.

Brownian dynamics simulations of protein folding: Access to milliseconds time scale and beyond

Protein folding occurs on a time scale ranging from milliseconds to minutes for a majority of proteins. Computer simulation of protein folding, from a random configuration to the native structure, is nontrivial owing to the large disparity between the simulation and folding time scales. As an effort to overcome this limitation, simple models with idealized protein subdomains, e.g., the diffusion–collision model of Karplus and Weaver, have gained some popularity. We present here new results for the folding of a four-helix bundle within the framework of the diffusion–collision model. Even with such simplifying assumptions, a direct application of standard Brownian dynamics methods would consume 10,000 processor-years on current supercomputers. We circumvent this difficulty by invoking a special Brownian dynamics simulation. The method features the calculation of the mean passage time of an event from the flux overpopulation method and the sampling of events that lead to productive collisions even if their probability is extremely small (because of large free-energy barriers that separate them from the higher probability events). Using these developments, we demonstrate that a coarse-grained model of the four-helix bundle can be simulated in several days on current supercomputers. Furthermore, such simulations yield folding times that are in the range of time scales observed in experiments.

In disperse solution, “osmotic stress” is a restricted case of preferential interactions

In the practice of “osmotic stress,” the effect of excluded cosolvents on a biochemical equilibrium is interpreted as the number of water molecules participating in the reaction. This action is attributed to lowering of solvent water activity by the cosolvent. This concept of osmotic stress in disperse solution is erroneous: (i) A cosolvent cannot be both excluded and inert, i.e., noninteracting, because exclusion requires a positive free energy change; (ii) a decrease in water activity alone by addition of solute cannot affect an equilibrium when the reacting surface is in contact with the solvent; and (iii) osmotic stress in disperse solution is a restricted case of preferential interactions; the reaction is driven by the free energy of cosolvent exclusion, and the derived number of water molecules is solely a measure of the mutual perturbations of the chemical potentials of the cosolvent and the protein.

Antigen binding forces of individually addressed single-chain Fv antibody molecules

Antibody single-chain Fv fragment (scFv) molecules that are specific for fluorescein have been engineered with a C-terminal cysteine for a directed immobilization on a flat gold surface. Individual scFv molecules can be identified by atomic force microscopy. For selected molecules the antigen binding forces are then determined by using a tip modified with covalently immobilized antigen. An scFv mutant of 12% lower free energy for ligand binding exhibits a statistically significant 20% lower binding force. This strategy of covalent immobilization and measuring well separated single molecules allows the characterization of ligand binding forces in molecular repertoires at the single molecule level and will provide a deeper insight into biorecognition processes.

Structure-based conformational preferences of amino acids

Proteins can be very tolerant to amino acid substitution, even within their core. Understanding the factors responsible for this behavior is of critical importance for protein engineering and design. Mutations in proteins have been quantified in terms of the changes in stability they induce. For example, guest residues in specific secondary structures have been used as probes of conformational preferences of amino acids, yielding propensity scales. Predicting these amino acid propensities would be a good test of any new potential energy functions used to mimic protein stability. We have recently developed a protein design procedure that optimizes whole sequences for a given target conformation based on the knowledge of the template backbone and on a semiempirical potential energy function. This energy function is purely physical, including steric interactions based on a Lennard-Jones potential, electrostatics based on a Coulomb potential, and hydrophobicity in the form of an environment free energy based on accessible surface area and interatomic contact areas. Sequences designed by this procedure for 10 different proteins were analyzed to extract conformational preferences for amino acids. The resulting structure-based propensity scales show significant agreements with experimental propensity scale values, both for α-helices and β-sheets. These results indicate that amino acid conformational preferences are a natural consequence of the potential energy we use. This confirms the accuracy of our potential and indicates that such preferences should not be added as a design criterion.

Strained DNA is kinked by low concentrations of Zn2+

A novel atomic force microscope with a magnetically oscillated tip has provided unprecedented resolution of small DNA fragments spontaneously adsorbed to mica and imaged in situ in the presence of divalent ions. Kinks (localized bends of average angle 78°) were observed in axially strained minicircles consisting of tandemly repeated d(A)5 and d(GGGCC[C]) sequences. The frequency of kinks in identical minicircles increased 4-fold in the presence of 1 mM Zn2+ compared with 1 mM Mg2+. Kinking persisted in mixed Mg2+/Zn2+ electrolytes until the Zn2+ concentration dropped below 100 μM, indicating that this type of kinking may occur under physiological conditions. Kinking appears to replace intrinsic bending, and statistical analysis shows that kinks are not localized within any single sequence element. A surprisingly small free energy is associated with kink formation.

Semirational design of active tumor suppressor p53 DNA binding domain with enhanced stability

We have designed a p53 DNA binding domain that has virtually the same binding affinity for the gadd45 promoter as does wild-type protein but is considerably more stable. The design strategy was based on molecular evolution of the protein domain. Naturally occurring amino acid substitutions were identified by comparing the sequences of p53 homologues from 23 species, introducing them into wild-type human p53, and measuring the changes in stability. The most stable substitutions were combined in a multiple mutant. The advantage of this strategy is that, by substituting with naturally occurring residues, the function is likely to be unimpaired. All point mutants bind the consensus DNA sequence. The changes in stability ranged from +1.27 (less stable Q165K) to −1.49 (more stable N239Y) kcal mol−1, respectively. The changes in free energy of unfolding on mutation are additive. Of interest, the two most stable mutants (N239Y and N268D) have been known to act as suppressors and restored the activity of two of the most common tumorigenic mutants. Of the 20 single mutants, 10 are cancer-associated, though their frequency of occurrence is extremely low: A129D, Q165K, Q167E, and D148E are less stable and M133L, V203A and N239Y are more stable whereas the rest are neutral. The quadruple mutant (M133LV203AN239YN268D), which is stabilized by 2.65 kcal mol−1 and Tm raised by 5.6°C is of potential interest for trials in vivo.

A streptavidin mutant with altered ligand-binding specificity

The biotin-binding site of streptavidin was modified to alter its ligand-binding specificity. In natural streptavidin, the side chains of N23 and S27 make two of the three hydrogen bonds with the ureido oxygen of biotin. These two residues were mutated to severely weaken biotin binding while attempting to maintain the affinity for two biotin analogs, 2-iminobiotin and diaminobiotin. Redesigning of the biotin-binding site used the difference in local electrostatic charge distribution between biotin and these biotin analogs. Free energy calculations predicted that the introduction of a negative charge at the position of S27 plus the mutation N23A should disrupt two of the three hydrogen bonds between natural streptavidin and the ureido oxygen of biotin. In contrast, the imino hydrogen of 2-iminobiotin should form a hydrogen bond with the side chain of an acidic amino acid at position 27. This should reduce the biotin-binding affinity by approximately eight orders of magnitude, while leaving the affinities for these biotin analogs virtually unaffected. In good agreement with these predictions, a streptavidin mutant with the N23A and S27D substitutions binds 2-iminobiotin with an affinity (Ka) of 1 × 106 M−1, two orders of magnitude higher than that for biotin (1 × 104 M−1). In contrast, the binding affinity of this streptavidin mutant for diaminobiotin (2.7 × 104 M−1) was lower than predicted (2.9 × 105 M−1), suggesting the position of the diaminobiotin in the biotin-binding site was not accurately determined by modeling.

Folding thermodynamics of a model three-helix-bundle protein

The calculated folding thermodynamics of a simple off-lattice three-helix-bundle protein model under equilibrium conditions shows the experimentally observed protein transitions: a collapse transition, a disordered-to-ordered globule transition, a globule to native-state transition, and the transition from the active native state to a frozen inactive state. The cooperativity and physical origin of the various transitions are explored with a single “optimization” parameter and characterized with the Lindemann criterion for liquid versus solid-state dynamics. Below the folding temperature, the model has a simple free energy surface with a single basin near the native state; the surface is similar to that calculated from a simulation of the same three-helix-bundle protein with an all-atom representation [Boczko, E. M. & Brooks III, C. L. (1995) Science 269, 393–396].

Peptide interfacial adsorption is kinetically limited by the thermodynamic stability of self association

We present a study of the adsorption of two peptides at the octane–water interface. The first peptide, Lac21, exists in mixed monomer–tetramer equilibrium in bulk solution with an appreciable monomer concentration. The second peptide, Lac28, exists as a tetramer in solution, with minimal exposed hydrophobic surface. A kinetic limitation to interfacial adsorption exists for Lac28 at moderate to high surface coverage that is not observed for Lac21. We estimate the potential energy barrier for Lac28 adsorption to be 42 kJ/mol and show that this is comparable to the expected free energy barrier for tetramer dissociation. This finding suggests that, at moderate to high surface coverage, adsorption is kinetically limited by the availability of interfacially active monomeric “domains” in the subinterfacial region. We also show how the commonly used empirical equation for protein adsorption dynamics can be used to estimate the potential energy barrier for adsorption. Such an approach is shown to be consistent with a formal description of diffusion–adsorption, provided a large potential energy barrier exists. This work demonstrates that the dynamics of interfacial adsorption depend on protein thermodynamic stability, and hence structure, in a quantifiable way.

Transition-state structure as a unifying basis in protein-folding mechanisms: Contact order, chain topology, stability, and the extended nucleus mechanism

I attempt to reconcile apparently conflicting factors and mechanisms that have been proposed to determine the rate constant for two-state folding of small proteins, on the basis of general features of the structures of transition states. Φ-Value analysis implies a transition state for folding that resembles an expanded and distorted native structure, which is built around an extended nucleus. The nucleus is composed predominantly of elements of partly or well-formed native secondary structure that are stabilized by local and long-range tertiary interactions. These long-range interactions give rise to connecting loops, frequently containing the native loops that are poorly structured. I derive an equation that relates differences in the contact order of a protein to changes in the length of linking loops, which, in turn, is directly related to the unfavorable free energy of the loops in the transition state. Kinetic data on loop extension mutants of CI2 and α-spectrin SH3 domain fit the equation qualitatively. The rate of folding depends primarily on the interactions that directly stabilize the nucleus, especially those in native-like secondary structure and those resulting from the entropy loss from the connecting loops, which vary with contact order. This partitioning of energy accounts for the success of some algorithms that predict folding rates, because they use these principles either explicitly or implicitly. The extended nucleus model thus unifies the observations of rate depending on both stability and topology.

Kinetic equilibrium of forces and molecular events in muscle contraction

In biomolecular systems, the mechanical transfer of free energy occurs with both high efficiency and high speed. It is shown here that such a transfer can be achieved only if the participating free-energy-storing elements exhibit opposing relationships between their content of free energy and the force they exert in the transfer direction. A kinetic equilibrium of forces (KEF) results, in which the transfer of free energy is mediated essentially by thermal molecular motion. On the basis of present evidence, KEF is used as a guiding principle in developing a mechanical model of the crossbridge cycle in muscle contraction. The model allows the basic features of molecular events to be visualized in terms of plausible structures. Real understanding of the process will require identification of the elements that perform the functions described here. Besides chemomechanical energy transduction, KEF may have a role in other biomolecular processes in which free energy is transferred mechanically over large distances.

A quantitative method for evaluating the stabilities of nucleic acids

We report a general method for screening, in solution, the impact of deviations from canonical Watson-Crick composition on the thermodynamic stability of nucleic acid duplexes. We demonstrate how fluorescence resonance energy transfer (FRET) can be used to detect directly free energy differences between an initially formed “reference” duplex (usually a Watson-Crick duplex) and a related “test” duplex containing a lesion/alteration of interest (e.g., a mismatch, a modified, a deleted, or a bulged base, etc.). In one application, one titrates into a solution containing a fluorescently labeled, FRET-active, reference duplex, an unlabeled, single-stranded nucleic acid (test strand), which may or may not compete successfully to form a new duplex. When a new duplex forms by strand displacement, it will not exhibit FRET. The resultant titration curve (normalized fluorescence intensity vs. logarithm of test strand concentration) yields a value for the difference in stability (free energy) between the newly formed, test strand-containing duplex and the initial reference duplex. The use of competitive equilibria in this assay allows the measurement of equilibrium association constants that far exceed the magnitudes accessible by conventional titrimetric techniques. Additionally, because of the sensitivity of fluorescence, the method requires several orders of magnitude less material than most other solution methods. We discuss the advantages of this method for detecting and characterizing any modification that alters duplex stability, including, but not limited to, mutagenic lesions. We underscore the wide range of accessible free energy values that can be defined by this method, the applicability of the method in probing for a myriad of nucleic acid variations, such as single nucleotide polymorphisms, and the potential of the method for high throughput screening.

Magnesium-dependent folding of self-splicing RNA: Exploring the link between cooperativity, thermodynamics, and kinetics

Folding of the Tetrahymena self-splicing RNA into its active conformation involves a set of discrete intermediate states. The Mg2+-dependent equilibrium transition from the intermediates to the native structure is more cooperative than the formation of the intermediates from the unfolded states. We show that the degree of cooperativity is linked to the free energy of each transition and that the rate of the slow transition from the intermediates to the native state decreases exponentially with increasing Mg2+ concentration. Monovalent salts, which stabilize the folded RNA nonspecifically, induce states that fold in less than 30 s after Mg2+ is added to the RNA. A simple model is proposed that predicts the folding kinetics from the Mg2+-dependent change in the relative stabilities of the intermediate and native states.

Stretching single-domain proteins: Phase diagram and kinetics of force-induced unfolding

Single-molecule force spectroscopy reveals unfolding of domains in titin on stretching. We provide a theoretical framework for these experiments by computing the phase diagrams for force-induced unfolding of single-domain proteins using lattice models. The results show that two-state folders (at zero force) unravel cooperatively, whereas stretching of non-two-state folders occurs through intermediates. The stretching rates of individual molecules show great variations reflecting the heterogeneity of force-induced unfolding pathways. The approach to the stretched state occurs in a stepwise “quantized” manner. Unfolding dynamics and forces required to stretch proteins depend sensitively on topology. The unfolding rates increase exponentially with force f till an optimum value, which is determined by the barrier to unfolding when f = 0. A mapping of these results to proteins shows qualitative agreement with force-induced unfolding of Ig-like domains in titin. We show that single-molecule force spectroscopy can be used to map the folding free energy landscape of proteins in the absence of denaturants.