Wonmuk Hwang

Design of Nanostructured Biological Materials Through Self-Assembly of Peptides and Proteins

Several self-assembling peptide and protein systems that form nanotubes, helical ribbons and fibrous scaffolds have recently emerged as biological materials. Peptides and proteins have also been selected to bind metals, semiconductors and ions, inspiring the design of new materials for a wide range of applications in nano-biotechnology.

Force generation in kinesin hinges on cover-neck bundle formation.

In kinesin motors, a fundamental question concerns the mechanism by which ATP binding generates the force required for walking. Analysis of available structures combined with molecular dynamics simulations demonstrates that the conformational change of the neck linker involves the nine-residue-long N-terminal region, the cover strand, as an element that is essential for force generation. Upon ATP binding, it forms a beta sheet with the neck linker, the cover-neck bundle, which induces the forward motion of the neck linker, followed by a latch-type binding to the motor head. The estimated stall force and anisotropic response to external loads calculated from the model agree with force-clamp measurements. The proposed mechanism for force generation by the cover-neck bundle formation appears to apply to several kinesin families. It also elucidates the design principle of kinesin as the smallest known processive motor.

Computational analysis of viscoelastic properties of crosslinked actin networks

Mechanical force plays an important role in the physiology of eukaryotic cells whose dominant structural constituent is the actin cytoskeleton composed mainly of actin and actin crosslinking proteins (ACPs). Thus, knowledge of rheological properties of actin networks is crucial for understanding the mechanics and processes of cells. We used Brownian dynamics simulations to study the viscoelasticity of crosslinked actin networks. Two methods were employed, bulk rheology and segment-tracking rheology, where the former measures the stress in response to an applied shear strain, and the latter analyzes thermal fluctuations of individual actin segments of the network. It was demonstrated that the storage shear modulus (G′) increases more by the addition of ACPs that form orthogonal crosslinks than by those that form parallel bundles. In networks with orthogonal crosslinks, as crosslink density increases, the power law exponent of G′ as a function of the oscillation frequency decreases from 0.75, which reflects the transverse thermal motion of actin filaments, to near zero at low frequency. Under increasing prestrain, the network becomes more elastic, and three regimes of behavior are observed, each dominated by different mechanisms: bending of actin filaments, bending of ACPs, and at the highest prestrain tested (55%), stretching of actin filaments and ACPs. In the last case, only a small portion of actin filaments connected via highly stressed ACPs support the strain. We thus introduce the concept of a ‘supportive framework,’ as a subset of the full network, which is responsible for high elasticity. Notably, entropic effects due to thermal fluctuations appear to be important only at relatively low prestrains and when the average crosslinking distance is comparable to or greater than the persistence length of the filament. Taken together, our results suggest that viscoelasticity of the actin network is attributable to different mechanisms depending on the amount of prestrain.

Kinesin's cover-neck bundle folds forward to generate force

Each step of the kinesin motor involves a force-generating molecular rearrangement. Although significant progress has been made in elucidating the broad features of the kinesin mechanochemical cycle, molecular details of the force generation mechanism remain a mystery. Recent molecular dynamics simulations have suggested a mechanism in which the forward drive is produced when the N-terminal cover strand forms a β-sheet with the neck linker to yield the cover-neck bundle. We tested this proposal by comparing optical trapping motility measurements of cover strand mutants with the wild-type. Motility data, as well as kinetic analyses, revealed impairment of the force-generating capacity accompanied by a greater load dependence in the mechanochemical cycle. In particular, a mutant with the cover strand deleted functioned only marginally, despite the fact that the cover strand, the N-terminal “dangling end,” unlike the neck linker and nucleotide-binding pocket, is not involved with any previously considered energy transduction pathway. Furthermore, a constant assisting load, likely in lieu of a power stroke, was shown to rescue forward motility in the cover strand deletion mutant. Our results support a stepping mechanism driven by dynamic cover-neck bundle formation. They also suggest a strategy to generate motors with altered mechanical characteristics by targeting the force-generating element.

Supramolecular structure of helical ribbons self-assembled from a β-sheet peptide

We have investigated the supramolecular structure of helical ribbons formed during self-assembly of a β-sheet peptide using computer simulation. We tested a wide range of molecular packing geometries consistent with the experimental dimensions to identify the most stable structure, and then systematically changed the helical geometry to investigate its energy landscape. The effect of pH was incorporated by scaling the amount of charge on the side chains based on the electrostatic double layer theory. Our results suggest that these left-handed helical ribbons are comprised of a double β-sheet and that the experimentally measured dimensions correspond to a local energy minimum. Side chain interactions are found to be critical in determining the stability and curvature of the helix. Our approach has general applicability to the study of self-assembled nanostructures from β-sheet peptides where high resolution data are not yet available.

Site-directed nanoparticle labeling of cytochrome c

Although nanoparticle-protein conjugates have been synthesized for numerous applications, bioconjugation remains a challenge, often resulting in denaturation or loss of protein function. This is partly because the protein–nanoparticle interface is poorly understood, which impedes the use of nanoparticles in nanomedicine. Although the effects of nanoparticle ligand and material on protein structure have been explored, the choice of the labeling site on the protein has not yet been systematically studied. To address this issue, we label cytochrome c site-specifically with a negatively charged Au nanoparticle via a covalent thiol–Au bond. The attachment site is controlled by cysteine mutations of surface residues. The effect of labeling on protein structure is probed by circular dichroism. Protein unfolding is the most severe when the nanoparticle is attached to the N- and C-terminal foldon, the core motif of cytochrome c. Also, when the nanoparticle is attached in the vicinity of charged residues, the amount of structural damage is greater because of salt-dependent electrostatic interactions with charged ligand bis(p-sulfonatophenyl) phenylphosphine on the nanoparticle. Molecular dynamics simulations also elucidate local to global structural perturbation depending on labeling site. These results suggest that the labeling site must be considered as one of the main design criteria for nanoparticle–protein conjugates.

Thermodynamic Selection of Steric Zipper Patterns in the Amyloid Cross-β Spine

At the core of amyloid fibrils is the cross-β spine, a long tape of β-sheets formed by the constituent proteins. Recent high-resolution x-ray studies show that the unit of this filamentous structure is a β-sheet bilayer with side chains within the bilayer forming a tightly interdigitating “steric zipper” interface. However, for a given peptide, different bilayer patterns are possible, and no quantitative explanation exists regarding which pattern is selected or under what condition there can be more than one pattern observed, exhibiting molecular polymorphism. We address the structural selection mechanism by performing molecular dynamics simulations to calculate the free energy of incorporating a peptide monomer into a β-sheet bilayer. We test filaments formed by several types of peptides including GNNQQNY, NNQQ, VEALYL, KLVFFAE and STVIIE, and find that the patterns with the lowest binding free energy correspond to available atomistic structures with high accuracy. Molecular polymorphism, as exhibited by NNQQ, is likely because there are more than one most stable structures whose binding free energies differ by less than the thermal energy. Detailed analysis of individual energy terms reveals that these short peptides are not strained nor do they lose much conformational entropy upon incorporating into a β-sheet bilayer. The selection of a bilayer pattern is determined mainly by the van der Waals and hydrophobic forces as a quantitative measure of shape complementarity among side chains between the β-sheets. The requirement for self-complementary steric zipper formation supports that amyloid fibrils form more easily among similar or same sequences, and it also makes parallel β-sheets generally preferred over anti-parallel ones. But the presence of charged side chains appears to kinetically drive anti-parallel β-sheets to form at early stages of assembly, after which the bilayer formation is likely driven by energetics.

Region‐specific role of water in collagen unwinding and assembly

Conformational stability of the collagen triple helix affects its turnover and determines tissue homeostasis. Although it is known that the presence of imino acids (prolines or hydroxyprolines) confer stability to the molecule, little is known regarding the stability of the imino‐poor region lacking imino acids, which plays a key role in collagen cleavage. In particular, there have been continuing debates about the role of water in collagen stability. We addressed these issues using molecular dynamics simulations on 30‐residue long collagen triple helices, including a structure that has a biologically relevant 9‐residue imino‐poor region from type III collagen (PDB ID: 1BKV). A torsional map approach was used to characterize the conformational motion of the molecule that differ between imino‐rich and imino‐poor regions. At temperatures 300 K and above, unwinding initiates at a common cleavage site, the glycine‐isoleucine bond in the imino‐poor region. This provides a linkage between previous observations that unwinding of the imino‐poor region is a requirement for collagenase cleavage, and that isolated collagen molecules are unstable at body temperature. We found that unwinding of the imino‐poor region is controlled by dynamic water bridges between backbone atoms with average lifetimes on the order of a few picoseconds, as the degree of unwinding strongly correlated with the loss of water bridges, and unwinding could be either prevented or enhanced, respectively by enforcing or forbidding water bridge formation. While individual water bridges were short‐lived in the imino‐poor region, the hydration shell surrounding the entire molecule was stable even at 330 K. The diameter of the hydrated collagen including the first hydration shell was about 14 Å, in good agreement with the experimentally measured inter‐collagen distances. These results elucidate the general role of water in collagen turnover: water not only affects collagen cleavage by controlling its torsional motion, but it also forms a larger‐scale lubrication layer mediating collagen self‐assembly. Proteins 2008.

Molecular basis for optical clearing of collagenous tissues

Molecular interactions of optical clearing agents were investigated using a combination of molecular dynamics (MD) simulations and optical spectroscopy. For a series of sugar alcohols with low to high optical clearing potential, Raman spectroscopy and integrating sphere measurements were used to quantitatively characterize tissue water loss and reduction in light scattering following agent exposures. The rate of tissue water loss was found to correlate with agent optical clearing potential, but equivalent tissue optical clearing was measured in native and fixed tissue in vitro, given long-enough exposure times to the polyol series. MD simulations showed that the rate of tissue optical clearing correlated with the preferential formation of hydrogen bond bridges between agent and collagen. Hydrogen bond bridge formation disrupts the collagen hydration layer and facilitates replacement by a chemical agent to homogenize tissue refractive index. However, the reduction in tissue light scattering did not correlate with the agent index of refraction. Our results suggest that a necessary property of optical clearing agents is hyperosmolarity to tissue, but that the most effective agents with the highest rates of optical clearing are a subset with the highest collagen solubilities.

Epitaxially guided assembly of collagen layers on mica surfaces.

Ordered assembly of collagen molecules on flat substrates has potential for various applications and serves as a model system for studying the assembly process. While previous studies demonstrated self-assembly of collagen on muscovite mica into highly ordered layers, the mechanism by which different conditions affect the resulting morphology remains to be elucidated. Using atomic force microscopy, we follow the assembly of collagen on muscovite mica at a concentration lower than the critical fibrillogenesis concentration in bulk. Initially, individual collagen molecules adsorb to mica and subsequently nucleate into fibrils possessing the 67 nm D-periodic bands. Emergence of fibrils aligned in parallel despite large interfibril distances agrees with an alignment mechanism guided by the underlying mica. The epitaxial growth was further confirmed by the formation of novel triangular networks of collagen fibrils on phlogopite mica, whose surface lattice is known to have a hexagonal symmetry, whereas the more widely used muscovite does not. Comparing collagen assembly on the two types of mica at different potassium concentrations revealed that potassium binds to the negatively charged mica surface and neutralizes it, thereby reducing the binding affinity of collagen and enhancing surface diffusion. These results suggest that collagen assembly on mica follows the surface adsorption, diffusion, nucleation, and growth pathway, where the growth direction is determined at the nucleation step. Comparison with other molecules that assemble similarly on mica supports generality of the proposed assembly mechanism, the knowledge of which will be useful for controlling the resulting surface morphologies.