Thermodynamic analysis of self-assembly in coiled-coil biomaterials.

Coiled-coil protein structural motifs have proven amenable to the design of structurally well-defined biomaterials. Mesoscale structural properties can be fairly well predicted based on rules governing the chemical interactions between the helices that define this structural motif. We explore the role of the hydrophobic core residues on the self-assembly of a coiled-coil polymer through a mutational analysis coupled with a salting-out procedure. Because the resultant polymers remain in solution, a thermodynamic approach is applied to characterize the polymer assembly using conventional equations from polymer theory to extract nucleation and elongation parameters. The stabilities and lengths of the polymers are measured using circular dichroism spectropolarimetry, sizing methods including dynamic light scattering and analytical ultracentrifugation, and atomic force microscopy to assess mesoscale morphology. Upon mutating isoleucines at two core positions to serines, we find that polymer stability is decreased while the degree of polymerization is about the same. Differences in results from circular dichroism and dynamic light scattering experiments suggest the presence of a stable intermediate state, and a scheme is proposed for how this intermediate might relate to the monomer and polymer states.

Energetic differences at the subunit interfaces of normal human hemoglobins correlate with their developmental profile.

A previously unrecognized function of normal human hemoglobins occurring during protein assembly is described, i.e. self-regulation of subunit pairings and their durations arising from the variable strengths of their subunit interactions. Although many mutant human hemoglobins are known to have altered subunit interface strengths, those of the normal embryonic, fetal, and adult human hemoglobins have not been considered to differ significantly. However, in a comprehensive study of both types of subunit interfaces of seven of the eight normal oxy human hemoglobins, we found that the strengths, i.e., the free energies of the tetramer-dimer interfaces, contrary to previous reports, differ by 3 orders of magnitude and display an undulating profile similar to the transitions ("switches") of various globin subunit types over time. The dimer interface strengths are also variable and correlate linearly with their developmental profile. Embryonic hemoglobins are the weakest; fetal hemoglobin is of intermediate strength, and adult hemoglobins are the strongest. The pattern also correlates generally with their different O(2) affinities and responses to allosteric regulatory molecules. Acetylation of fetal hemoglobin weakens its unusually strong subunit interactions and occurs progressively as its level of expression diminishes and adult hemoglobin A formation begins; a causal relationship is suggested. The relative contributions of globin gene order and competition among subunits due to differences in their interface strengths were found to be complementary and establish a connection among genetics, thermodynamics, and development.

Quantifying the impact of simple DNA parameters on the cyclization J-factor for single-basepair-addition families.

We use Monte Carlo simulation to quantify the change in cyclization J-factor within a dramatically simplified model of DNA that involves parameters for uniform stiffnesses, intrinsic twist, and intrinsic bending (including nonplanar bending). Plots of J versus DNA length over multiple periods of helical repeat are fit to a simple functional form in order to project the behavior of J over a broad range of these model parameters. In some instances, this process allows us to find families of DNA molecules (within our model) with quite different material properties, but very similar plots of J versus length, so similar as to likely to be indistinguishable by experiments. This effect is seen both for the parameter-pair of bend angle and stiffness scaling, as well as for the parameter-trio of helical repeat, bend angle, and bend non-planarity.

Human embryonic, fetal, and adult hemoglobins have different subunit interface strengths. Correlation with lifespan in the red cell.

The different types of naturally occurring, normal human hemoglobins vary in their tetramer-dimer subunit interface strengths (stabilities) by three orders of magnitude in the liganded (CO or oxy) state. The presence of embryonic zeta-subunits leads to an average 20-fold weakening of tetramer-dimer interfaces compared to corresponding hemoglobins containing adult alpha-subunits. The dimer-monomer interfaces of these hemoglobins differ by at least 500-fold in their strengths; such interfaces are weak if they contain zeta-subunits and exchange with added beta-subunits in the form of beta(4) (HbH) significantly faster than do those with alpha-subunits. Subunit exchange occurs at the level of the dimer, although tetramer formation reciprocally influences the amount of dimer available for exchange. Competition between subunit types occurs so that pairs of weak embryonic hemoglobins can exchange subunits to form the stronger fetal and adult hemoglobins. The dimer strengths increase in the order Hb Portland-2 (zeta(2)beta(2)) < Hb Portland-1 (zeta(2)gamma(2)) approximately equal Hb Gower-1 (zeta(2)epsilon(2)) < Hb Gower-2 (alpha(2)epsilon(2)) < HbF(1) < HbF (alpha(2)gamma(2)) < HbA(2) (alpha(2)delta(2)), i.e., from embryonic to fetal to adult types, representing maturation from weaker to stronger monomer-monomer subunit contacts. This increasing order recapitulates the developmental order in which globins are expressed (embryonic --> fetal --> adult), suggesting that the intrinsic binding properties of the subunits themselves regarding the strengths of interfaces they form with competing subunits play an important role in the dynamics of protein assemblies and networks.

Sequence-Dependent Persistence Lengths of DNA.

A Monte Carlo code applied to the cgDNA coarse-grain rigid-base model of B-form double-stranded DNA is used to predict a sequence-averaged persistence length of lF = 53.5 nm in the sense of Flory, and of lp = 160 bp or 53.5 nm in the sense of apparent tangent-tangent correlation decay. These estimates are slightly higher than the consensus experimental values of 150 bp or 50 nm, but we believe the agreement to be good given that the cgDNA model is itself parametrized from molecular dynamics simulations of short fragments of length 10-20 bp, with no explicit fit to persistence length. Our Monte Carlo simulations further predict that there can be substantial dependence of persistence lengths on the specific sequence [Formula: see text] of a fragment. We propose, and confirm the numerical accuracy of, a simple factorization that separates the part of the apparent tangent-tangent correlation decay [Formula: see text] attributable to intrinsic shape, from a part [Formula: see text] attributable purely to stiffness, i.e., a sequence-dependent version of what has been called sequence-averaged dynamic persistence length l̅d (=58.8 nm within the cgDNA model). For ensembles of both random and λ-phage fragments, the apparent persistence length [Formula: see text] has a standard deviation of 4 nm over sequence, whereas our dynamic persistence length [Formula: see text] has a standard deviation of only 1 nm. However, there are notable dynamic persistence length outliers, including poly(A) (exceptionally straight and stiff), poly(TA) (tightly coiled and exceptionally soft), and phased A-tract sequence motifs (exceptionally bent and stiff). The results of our numerical simulations agree reasonably well with both molecular dynamics simulation and diverse experimental data including minicircle cyclization rates and stereo cryo-electron microscopy images.

Diagnostic schemes for reducing epidemic size of African viral hemorrhagic fever outbreaks.

INTRODUCTION: Viral hemorrhagic fever (VHF) outbreaks, with high mortality rates, have often been amplified in African health institutions due to person-to-person transmission via infected body fluids.  By collating and analyzing epidemiological data from documented outbreaks, we observed that diagnostic delay contributes to epidemic size for Ebola and Marburg hemorrhagic fever outbreaks. METHODOLOGY: We used a susceptible-exposed-infectious-removed (SEIR) model and data from the 1995 outbreak in Kikwit, Democratic Republic of Congo, to simulate Ebola hemorrhagic fever epidemics. Our model allows us to describe the dynamics for hospital staff separately from that for the general population, and to implement health worker-specific interventions. RESULTS: The model illustrates that implementing World Health Organization/US Centers for Disease Control and Prevention guidelines of isolating patients who do not respond to antimalarial and antibacterial chemotherapy reduces total outbreak size, from a median of 236, by 90% or more. Routinely employing diagnostic testing in post-mortems of patients that died of refractory fevers reduces the median outbreak size by a further 60%. Even greater reductions in outbreak size were seen when all febrile patients were tested for endemic infections or when febrile health-care workers were tested.  The effect of testing strategies was not impaired by the 1-3 day delay that would occur if testing were performed by a reference laboratory. CONCLUSION: In addition to improving the quality of care for common causes of febrile infections, increased and strategic use of laboratory diagnostics for fever could reduce the chance of hospital amplification of VHFs in resource-limited African health systems.

Quantitative atomic force microscopy image analysis of unusual filaments formed by the Acanthamoeba castellanii myosin II rod domain.

We describe a quantitative analysis of Acanthamoeba castellanii myosin II rod domain images collected from atomic force microscope experiments. These images reveal that the rod domain forms a novel filament structure, most likely requiring unusual head-to-tail interactions. Similar filaments are seen also in negatively stained electron microscopy images. Truncated myosins from Acanthamoeba and other model organisms have been visualized before, revealing laterally associated bipolar minifilaments. In contrast, the filament structures that we observe are dominated by axial rather than lateral polymerization. The unusually small features in this structure (1-5 nm) required the development of quantitative and statistical techniques for filament image analysis. These techniques enhance the extraction of features that hitherto have been difficult to ascertain from more qualitative imaging approaches. The heights of the filaments are observed to have a bimodal distribution consistent with the diameters of a single rod domain and a pair of close-packed rod domains. Further quantitative analysis indicates that in-plane association is limited to at most a pair of rod domains. Taken together, this implies that the filaments contain no more than four rod domains laterally associated with one another, somewhat less than that seen in bipolar minifilaments. Analysis of images of the filaments decorated with an anti-FLAG antibody reveals head-to-tail association with mean distances between the antibodies of 75 +/- 15 nm. We consider a set of molecular models to help interpret possible structures of the filaments.

Link, twist, energy, and the stability of DNA minicircles.

We describe how the stability properties of DNA minicircles can be directly read from plots of various biologically intuitive quantities along families of equilibrium configurations. Our conclusions follow from extensions of the mathematical theory of distinguished bifurcation diagrams that are applied within the specific context of an elastic rod model of minicircles. Families of equilibria arise as a twisting angle alpha is varied. This angle is intimately related to the continuously varying linking number Lk for nicked DNA configurations that is defined as the sum of Twist and Writhe. We present several examples of such distinguished bifurcation diagrams involving plots of the energy E, linking number Lk, and a twist moment m3, along families of cyclized equilibria of both intrinsically straight and intrinsically curved DNA fragments.