Moon Ki Kim

Efficient Generation of Feasible Pathways for Protein Conformational Transitions

We develop a computationally efficient method to simulate the transition of a protein between two conformations. Our method is based on a coarse-grained elastic network model in which distances between spatially proximal amino acids are interpolated between the values specified by the two end conformations. The computational speed of this method depends strongly on the choice of cutoff distance used to define interactions as measured by the density of entries of the constant linking/contact matrix. To circumvent this problem we introduce the concept of using a cutoff based on a maximum number of nearest neighbors. This generates linking matrices that are both sparse and uniform, hence allowing for efficient computations that are independent of the arbitrariness of cutoff distance choices. Simulation results demonstrate that the method developed here reliably generates feasible intermediate conformations, because our method observes steric constraints and produces monotonic changes in virtual bond and torsion angles. Applications are readily made to large proteins, and we demonstrate our method on lactate dehydrogenase, citrate synthase, and lactoferrin. We also illustrate how this framework can be used to complement experimental techniques that partially observe protein motions.

Elastic models of conformational transitions in macromolecules.

We develop a computationally efficient and physically realistic method to simulate the transition of a macromolecule between two conformations. Our method is based on a coarse-grained elastic network model in which contact interactions between spatially proximal parts of the macromolecule are modelled with Gaussian/harmonic potentials. To delimit the interactions in such models, we introduce a cutoff to the permitted number of nearest neighbors. This generates stiffness (Hessian) matrices that are both sparse and quite uniform, hence, allowing for efficient computations. Several toy models are tested using our method to mimic simple classes of macromolecular motions such as stretching, hinge bending, shear, compression, ligand binding and nucleic acid structural transitions. Simulation results demonstrate that the method developed here reliably generates sequences of feasible intermediate conformations of macromolecules, since our method observes steric constraints and produces monotonic changes to virtual bond angles and torsion angles. A final application is made to the opening process of the protein lactoferrin.

Efficient Transfer of Large-Area Graphene Films onto Rigid Substrates by Hot Pressing

Graphene films grown on metal substrates by chemical vapor deposition (CVD) method have to be safely transferred onto desired substrates for further applications. Recently, a roll-to-roll (R2R) method has been developed for large-area transfer, which is particularly efficient for flexible target substrates. However, in the case of rigid substrates such as glass or wafers, the roll-based method is found to induce considerable mechanical damages on graphene films during the transfer process, resulting in the degradation of electrical property. Here we introduce an improved dry transfer technique based on a hot-pressing method that can minimize damage on graphene by neutralizing mechanical stress. Thus, we enhanced the transfer efficiency of the large-area graphene films on a substrate with arbitrary thickness and rigidity, evidenced by scanning electron microscope (SEM) and atomic force microscope (AFM) images, Raman spectra, and various electrical characterizations. We also performed a theoretical multiscale simulation from continuum to atomic level to compare the mechanical stresses caused by the R2R and the hot-pressing methods, which also supports our conclusion. Consequently, we believe that the proposed hot-pressing method will be immediately useful for display and solar cell applications that currently require rigid and large substrates.

Characterization of FBG sensor interrogation based on a FDML wavelength swept laser.

In this study, we develop an ultra-fast fiber Bragg grating sensor system that is based on the Fourier domain mode-locked (FDML) swept laser. A FDML wavelength swept laser has many advantages compared to the conventional wavelength swept laser source, such as high-speed interrogation, narrow spectral sensitivity, and high phase stability. The newly developed FDML wavelength swept laser shows a superior performance of a high scan rate of 31.3 kHz and a broad scan range of over 70 nm simultaneously. The performance of the grating sensor interrogating system using a FDML wavelength swept laser is characterized in both static and dynamic strain responses.

An elastic network model of HK97 capsid maturation

The structure of the capsid of bacteriophage HK97 has been solved at various stages of maturity by crystallography and cryo-electron microscopy, and has been reported previously in the literature. Typically the capsid assembles through polymerization and maturation processes. Maturation is composed of proteolytic cleavages to the precursor capsid (called Prohead II), expansion triggered by DNA packaging (in which the largest conformational changes of the capsid appear), and covalent cross-links of neighboring subunits to create the mature capsid called Head II. We apply a coarse-grained elastic network interpolation (ENI) to generate a feasible pathway for conformational change from Prohead II to Head II. The icosahedral symmetry of the capsid structure offers a significant computational advantage because it is not necessary to consider the whole capsid structure but only an asymmetric unit consisting of one hexamer plus an additional subunit from an adjacent pentamer. We also analyze normal modes of the capsid structure using an elastic network model which is also subject to symmetry constraints. Using our model, we can visualize the smooth evolution of capsid expansion and revisit in more detail several interesting geometric changes recognized in early experimental works such as rigid body motion of two compact domains (A and P) with two refolding extensions (N-arm and E-loop) and track the approach of the two particular residues associated with isopeptide bonds that make hexagonal cross-links in Head II. The feasibility of the predicted pathway is also supported by the results of our normal mode analysis.

A comparison between elastic network interpolation and MD simulation of 16S ribosomal RNA.

Abstract In this paper a coarse-grained method called elastic network interpolation (ENI) is used to generate feasible transition pathways between two given conformations of the core central domain of 16S Ribosomal RNA (16S rRNA). The two given conformations are the extremes generated by a molecular dynamics (MD) simulation, which differ from each other by 10Å in root-mean-square deviation (RMSD). It takes only several hours to build an ENI pathway on a 1.5GHz Pentium with 512 MB memory, while the MD takes several weeks on high-performance multi-processor servers such as the SGI ORIGIN 2000/2100. It is shown that multiple ENI pathways capture the essential anharmonic motions of millions of timesteps in a particular MD simulation. A coarse-grained normal mode analysis (NMA) is performed on each intermediate ENI conformation, and the lowest 1% of the normal modes (representing about 40 degrees of freedom (DOF)) are used to parameterize fluctuations. This combined ENI/NMA method captures all intermediate conformations in the MD run with 1.5Å RMSD on average. In addition, if we restrict attention to the time interval of the MD run between the two extreme conformations, the RMSD between the closest ENI/NMA pathway and the MD results is about 1Å. These results may serve as a paradigm for reducedDOF dynamic simulations of large biological macromolecules as well as a method for the reduced-parameter interpretation of massive amounts of MD data.

KOSMOS: a universal morph server for nucleic acids, proteins and their complexes

KOSMOS is the first online morph server to be able to address the structural dynamics of DNA/RNA, proteins and even their complexes, such as ribosomes. The key functions of KOSMOS are the harmonic and anharmonic analyses of macromolecules. In the harmonic analysis, normal mode analysis (NMA) based on an elastic network model (ENM) is performed, yielding vibrational modes and B-factor calculations, which provide insight into the potential biological functions of macromolecules based on their structural features. Anharmonic analysis involving elastic network interpolation (ENI) is used to generate plausible transition pathways between two given conformations by optimizing a topology-oriented cost function that guarantees a smooth transition without steric clashes. The quality of the computed pathways is evaluated based on their various facets, including topology, energy cost and compatibility with the NMA results. There are also two unique features of KOSMOS that distinguish it from other morph servers: (i) the versatility in the coarse-graining methods and (ii) the various connection rules in the ENM. The models enable us to analyze macromolecular dynamics with the maximum degrees of freedom by combining a variety of ENMs from full-atom to coarse-grained, backbone and hybrid models with one connection rule, such as distance-cutoff, number-cutoff or chemical-cutoff. KOSMOS is available at http://bioengineering.skku.ac.kr/kosmos.

A mass weighted chemical elastic network model elucidates closed form domain motions in proteins

An elastic network model (ENM), usually Cα coarse‐grained one, has been widely used to study protein dynamics as an alternative to classical molecular dynamics simulation. This simple approach dramatically saves the computational cost, but sometimes fails to describe a feasible conformational change due to unrealistically excessive spring connections. To overcome this limitation, we propose a mass‐weighted chemical elastic network model (MWCENM) in which the total mass of each residue is assumed to be concentrated on the representative alpha carbon atom and various stiffness values are precisely assigned according to the types of chemical interactions. We test MWCENM on several well‐known proteins of which both closed and open conformations are available as well as three α‐helix rich proteins. Their normal mode analysis reveals that MWCENM not only generates more plausible conformational changes, especially for closed forms of proteins, but also preserves protein secondary structures thus distinguishing MWCENM from traditional ENMs. In addition, MWCENM also reduces computational burden by using a more sparse stiffness matrix.

Spectrally-sampled OCT for sensitivity improvement from limited optical power

Although high optical illumination power is favored in optical coherence tomography (OCT) for better signal-to-noise ratio, optical power is often limited by a damaged threshold for biomedical living tissues and autocorrelation signals observed in tomograms. In order to improve signal sensitivity without increasing the optical illumination power, a spectrally sampled multi-wavelength light source is proposed for the OCT system. A fiber Sagnac comb filter was used to spectrally sample the output of a continuous spectral light source. Point spread function analysis shows that the spectrally sampled OCT has an almost 50% dynamic range improvement in comparison with a conventional continuous spectral light source OCT for the same average optical power of 6 mW.

Investigation of Ion Channel Activities of Gramicidin A in the Presence of Ionic Liquids Using Model Cell Membranes.

Ionic liquids (ILs) are considered to be green solvents because of their non-volatility. Although ILs are relatively safe in the atmospheric environment, they may be toxic in other environments. Our previous research showed that the cytotoxicity of ILs to biological organisms is attributable to interference with cell membranes by IL insertion. However, the effects of ILs on ion channels, which play important roles in cell homeostasis, have not been comprehensively studied to date. In this work, we studied the interactions between ILs and lipid bilayer membranes with gramicidin A ion channels. We used two methods, namely electrical and fluorescence measurements of ions that permeate the membrane. The lifetimes of channels were increased by all the ILs tested in this work via stabilizing the compressed structure of the lipid bilayer and the rate of ion flux through gA channels was decreased by changing the membrane surface charge. The former effect, which increased the rate of ion flux, was dominant at high salt concentrations, whereas the latter, which decreased the rate of ion flux, was dominant at low salt concentrations. The effects of ILs increased with increasing concentration and alkyl chain length. The experimental results were further studied using molecular dynamics simulations.