Seung-gu Kang

Computational design of a protein crystal

Protein crystals have catalytic and materials applications and are central to efforts in structural biology and therapeutic development. Designing predetermined crystal structures can be subtle given the complexity of proteins and the noncovalent interactions that govern crystallization. De novo protein design provides an approach to engineer highly complex nanoscale molecular structures, and often the positions of atoms can be programmed with sub-Å precision. Herein, a computational approach is presented for the design of proteins that self-assemble in three dimensions to yield macroscopic crystals. A three-helix coiled-coil protein is designed de novo to form a polar, layered, three-dimensional crystal having the P6 space group, which has a “honeycomb-like” structure and hexameric channels that span the crystal. The approach involves: (i) creating an ensemble of crystalline structures consistent with the targeted symmetry; (ii) characterizing this ensemble to identify “designable” structures from minima in the sequence-structure energy landscape and designing sequences for these structures; (iii) experimentally characterizing candidate proteins. A 2.1 Å resolution X-ray crystal structure of one such designed protein exhibits sub-Å agreement [backbone root mean square deviation (rmsd)] with the computational model of the crystal. This approach to crystal design has potential applications to the de novo design of nanostructured materials and to the modification of natural proteins to facilitate X-ray crystallographic analysis.

How force unfolding differs from chemical denaturation

Significance Proteins can be unfolded either chemically by denaturing agents (such as urea) or by mechanical force (such as optical tweezers). If these two methods produce quite different unfolded structures, the course of the subsequent protein refolding when these agents are removed will be quite different. Our computer simulations show that the respective unfolded states are indeed very different, being fully extended under force, with no contacts, and partially extended in urea with many nonnative contacts. Understanding these differences is crucial because such configurations serve as starting points for different single-molecule experimental studies of protein folding. These differences can affect both the folding mechanism and the kinetics observed, leading to very different interpretations. Single-molecule force spectroscopies are remarkable tools for studying protein folding and unfolding, but force unfolding explores protein configurations that are potentially very different from the ones traditionally explored in chemical or thermal denaturation. Understanding these differences is crucial because such configurations serve as starting points of folding studies, and thus can affect both the folding mechanism and the kinetics. Here we provide a detailed comparison of both chemically induced and force-induced unfolded state ensembles of ubiquitin based on extensive, all-atom simulations of the protein either extended by force or denatured by urea. As expected, the respective unfolded states are very different on a macromolecular scale, being fully extended under force with no contacts and partially extended in urea with many nonnative contacts. The amount of residual secondary structure also differs: A significant population of α-helices is found in chemically denatured configurations but such helices are absent under force, except at the lowest applied force of 30 pN where short helices form transiently. We see that typical-size helices are unstable above this force, and β-sheets cannot form. More surprisingly, we observe striking differences in the backbone dihedral angle distributions for the protein unfolded under force and the one unfolded by denaturant. A simple model based on the dialanine peptide is shown to not only provide an explanation for these striking differences but also illustrates how the force dependence of the protein dihedral angle distributions give rise to the worm-like chain behavior of the chain upon force.

Non-destructive Inhibition of Metallofullerenol Gd@C 82 (OH) 22 on WW domain: Implication on Signal Transduction Pathway

Endohedral metallofullerenol Gd@C82(OH)22 has recently been shown to effectively inhibit tumor growth; however, its potential adverse bioeffects remain to be understood before its wider applications. Here, we present our study on the interaction between Gd@C82(OH)22 and WW domain, a representative protein domain involved in signaling and regulatory pathway, using all-atom explicit solvent molecular dynamics simulations. We find that Gd@C82(OH)22 has an intrinsic binding preference to the binding groove, particularly the key signature residues Y28 and W39. In its binding competition with the native ligand PRM, Gd@C82(OH)22 is shown to easily win the competition over PRM in occupying the active site, implying that Gd@C82(OH)22 can impose a potential inhibitory effect on the WW domain. Further analyses with binding free energy landscapes reveal that Gd@C82(OH)22 can not only directly block the binding site of the WW domain, but also effectively distract the PRM from its native binding pocket.

Structural and electronic properties of uranium-encapsulated Au 14 cage

The structural properties of the uranium-encapsulated nano-cage U@Au14 are predicted using density functional theory. The presence of the uranium atom makes the Au14 structure more stable than the empty Au14-cage, with a triplet ground electronic state for U@Au14. Analysis of the electronic structure shows that the two frontier single-occupied molecular orbital electrons of U@Au14 mainly originate from the 5f shell of the U atom after charge transfer. Meanwhile, the bonding orbitals and charge population indicate that the designed U@Au14 nano-cage structure is stabilized by ionocovalent interactions. The current findings provide theoretical basis for future syntheses and further study of actinide doped gold nanoclusters, which might subsequently facilitate applications of such structure in radio-labeling, nanodrug carrier and other biomedical applications.

Salts drive controllable multilayered upright assembly of amyloid-like peptides at mica/water interface

Surface-assisted self-assembly of amyloid-like peptides has received considerable interest in both amyloidosis research and nanotechnology in recent years. Despite extensive studies, some controlling factors, such as salts, are still not well understood, even though it is known that some salts can promote peptide self-assemblies through the so-called “salting-out” effect. However, they are usually noncontrollable, disordered, amorphous aggregates. Here, we show via a combined experimental and theoretical approach that a conserved consensus peptide NH2-VGGAVVAGV-CONH2 (GAV-9) (from representative amyloidogenic proteins) can self-assemble into highly ordered, multilayered nanofilaments, with surprising all-upright conformations, under high-salt concentrations. Our atomic force microscopy images also demonstrate that the vertical stacking of multiple layers is highly controllable by tuning the ionic strength, such as from 0 mM (monolayer) to 100 mM (mainly double layer), and to 250 mM MgCl2 (double, triple, quadruple, and quintuple layers). Our atomistic molecular dynamics simulations then reveal that these individual layers have very different internal nanostructures, with parallel β-sheets in the first monolayer but antiparallel β-sheets in the subsequent upper layers due to their different microenvironment. Further studies show that the growth of multilayered, all-upright nanostructures is a common phenomenon for GAV-9 at the mica/water interface, under a variety of salt types and a wide range of salt concentrations.

Large scale molecular simulations of nanotoxicity

The widespread use of nanomaterials in biomedical applications has been accompanied by an increasing interest in understanding their interactions with tissues, cells, and biomolecules, and in particular, on how they might affect the integrity of cell membranes and proteins. In this mini‐review, we present a summary of some of the recent studies on this important subject, especially from the point of view of large scale molecular simulations. The carbon‐based nanomaterials and noble metal nanoparticles are the main focus, with additional discussions on quantum dots and other nanoparticles as well. The driving forces for adsorption of fullerenes, carbon nanotubes, and graphene nanosheets onto proteins or cell membranes are found to be mainly hydrophobic interactions and the so‐called π–π stacking (between aromatic rings), while for the noble metal nanoparticles the long‐range electrostatic interactions play a bigger role. More interestingly, there are also growing evidences showing that nanotoxicity can have implications in de novo design of nanomedicine. For example, the endohedral metallofullerenol Gd@C82(OH)22 is shown to inhibit tumor growth and metastasis by inhibiting enzyme MMP‐9, and graphene is illustrated to disrupt bacteria cell membranes by insertion/cutting as well as destructive extraction of lipid molecules. These recent findings have provided a better understanding of nanotoxicity at the molecular level and also suggested therapeutic potential by using the cytotoxicity of nanoparticles against cancer or bacteria cells. WIREs Syst Biol Med 2014, 6:265–279. doi: 10.1002/wsbm.1271

Interaction of Amyloid Inhibitor Proteins with Amyloid Beta Peptides: Insight from Molecular Dynamics Simulations

Knowledge of the detailed mechanism by which proteins such as human αB- crystallin and human lysozyme inhibit amyloid beta (Aβ) peptide aggregation is crucial for designing treatment for Alzheimers disease. Thus, unconstrained, atomistic molecular dynamics simulations in explicit solvent have been performed to characterize the Aβ17–42 assembly in presence of the αB-crystallin core domain and of lysozyme. Simulations reveal that both inhibitor proteins compete with inter-peptide interaction by binding to the peptides during the early stage of aggregation, which is consistent with their inhibitory action reported in experiments. However, the Aβ binding dynamics appear different for each inhibitor. The binding between crystallin and the peptide monomer, dominated by electrostatics, is relatively weak and transient due to the heterogeneous amino acid distribution of the inhibitor surface. The crystallin-bound Aβ oligomers are relatively long-lived, as they form more extensive contact surface with the inhibitor protein. In contrast, a high local density of arginines from lysozyme allows strong binding with Aβ peptide monomers, resulting in stable complexes. Our findings not only illustrate, in atomic detail, how the amyloid inhibitory mechanism of human αB-crystallin, a natural chaperone, is different from that of human lysozyme, but also may aid de novo design of amyloid inhibitors.

The complex and specific pMHC interactions with diverse HIV-1 TCR clonotypes reveal a structural basis for alterations in CTL function.

Immune control of viral infections is modulated by diverse T cell receptor (TCR) clonotypes engaging peptide-MHC class I complexes on infected cells, but the relationship between TCR structure and antiviral function is unclear. Here we apply in silico molecular modeling with in vivo mutagenesis studies to investigate TCR-pMHC interactions from multiple CTL clonotypes specific for a well-defined HIV-1 epitope. Our molecular dynamics simulations of viral peptide-HLA-TCR complexes, based on two independent co-crystal structure templates, reveal that effective and ineffective clonotypes bind to the terminal portions of the peptide-MHC through similar salt bridges, but their hydrophobic side-chain packings can be very different, which accounts for the major part of the differences among these clonotypes. Non-specific hydrogen bonding to viral peptide also accommodates greater epitope variants. Furthermore, free energy perturbation calculations for point mutations on the viral peptide KK10 show excellent agreement with in vivo mutagenesis assays, with new predictions confirmed by additional experiments. These findings indicate a direct structural basis for heterogeneous CTL antiviral function.

Dual Inhibitory Pathways of Metallofullerenol Gd@C82(OH)22 on Matrix Metalloproteinase-2: Molecular insight into drug-like nanomedicine

Cancer metastasis is an important criterion to evaluate tumor malignancy. Matrix metalloproteinases (MMPs) play a crucial role in cancer proliferation and migration by virtue of their proteolytic functions in angiogenesis and extracelluar matrix (ECM) degradation, making them potential targets of anti-metastaic therapeutics. Recently we showed with both in vivo and in vitro experiments that metallofullerenol Gd@C82(OH)22 can effectively inhibit MMP-2 and MMP-9 with high antitumoral efficacy. Furthermore, our in silico study revealed that Gd@C82(OH)22 could indirectly inhibit the proteolysis of MMP-9 via allosteric modulation exclusively at the ligand specificity S1′ loop. Here, we expand our study toward another gelatinase, MMP-2, using molecular dynamics simulations. Despite the high structural similarity with 64.3% sequence identity, their responses to Gd@C82(OH)22 were quite different. Toward MMP-2, Gd@C82(OH)22 could block either the Zn2+-catalylitic site directly or the S1′ loop indirectly. Surface electrostatics uniquely determines the initial adsorption of Gd@C82(OH)22 on MMP-2, and then its further location of the most favorable binding site(s). These findings not only illustrated how the inhibitory mechanism of Gd@C82(OH)22 is distinguished between the two gelatinase MMPs with atomic details, but also shed light on the de novo design of anti-metastatic nanotherapeutics with enhanced target specificity.

Robust Denaturation of Villin Headpiece by MoS2 Nanosheet: Potential Molecular Origin of the Nanotoxicity

MoS2 nanosheet, a new two-dimensional transition metal dichalcogenides nanomaterial, has attracted significant attentions lately due to many potential promising biomedical applications. Meanwhile, there is also a growing concern on its biocompatibility, with little known on its interactions with various biomolecules such as proteins. In this study, we use all-atom molecular dynamics simulations to investigate the interaction of a MoS2 nanosheet with Villin Headpiece (HP35), a model protein widely used in protein folding studies. We find that MoS2 exhibits robust denaturing capability to HP35, with its secondary structures severely destroyed within hundreds of nanosecond simulations. Both aromatic and basic residues are critical for the protein anchoring onto MoS2 surface, which then triggers the successive protein unfolding process. The main driving force behind the adsorption process is the dispersion interaction between protein and MoS2 monolayer. Moreover, water molecules at the interface between some key hydrophobic residues (e.g. Trp-64) and MoS2 surface also help to accelerate the process driven by nanoscale drying, which provides a strong hydrophobic force. These findings might have shed new light on the potential nanotoxicity of MoS2 to proteins with atomic details, which should be helpful in guiding future biomedical applications of MoS2 with its nanotoxicity mitigated.