Prem P. Chapagain

Intricate functions of matrix metalloproteinases in physiological and pathological conditions

Matrix metalloproteinases (MMPs) are a diverse group of proteolytic enzymes and play an important role in the degradation and remodeling of the extracellular matrix (ECM). In normal physiological conditions, MMPs are usually minimally expressed. Despite their low expression, MMPs have been implicated in many cellular processes ranging from embryological development to apoptosis. The activity of MMPs is controlled at three different stages: (1) transcription; (2) zymogen activation; and (3) inhibition of active forms by tissue inhibitor metalloproteinases (TIMPs). They can collectively degrade any component of ECM and basement membrane, and their excessive activity has been linked to numerous pathologies mainly including, but not limited to, tumor invasion and metastasis. The lack of information about several MMPs and the steady stream of new discoveries suggest that there is much more to be studied in this field. In particular, there is a need for controlling their expression in disease states. Various studies over the past 30 years have found that each MMP has a specific mode of activation, action, and inhibition. Drugs specifically targeting individual MMPs could revolutionize the treatment of a great number of health conditions and tremendously reduce their burden. In this review article, we have summarized the recent advances in understanding the role of MMPs in physiological and pathological conditions. J. Cell. Physiol. 231: 2599–2621, 2016.

Self-organization in protein folding and the hydrophobic interaction.

Self-organization is a critical aspect of living systems. During the folding of protein molecules, the hydrophobic interaction plays an important role in the collapse of the peptide chain to a compact shape. As the hydrophobic core tightens and excludes water, not only does the number of hydrophobic side chain contacts increase, but stabilization is further enhanced by an increase in strength of each hydrophobic interaction between side chains in the core. Thus, the self-organization of the protein folding process augments itself by enhancing the stability of the core against large-scale motions that would unfold the protein. Through calculations and computer simulations on a model four-helix bundle protein, we show how the strengthening of the hydrophobic interaction is crucial for stabilizing the core long enough for completion of the folding process and quantitatively manifests self-organizing dynamical behavior.

Fluorescent protein barrel fluctuations and oxygen diffusion pathways in mCherry

Fluorescent proteins (FPs) are valuable tools as biochemical markers for studying cellular processes. Red fluorescent proteins (RFPs) are highly desirable for in vivo applications because they absorb and emit light in the red region of the spectrum where cellular autofluorescence is low. The naturally occurring fluorescent proteins with emission peaks in this region of the spectrum occur in dimeric or tetrameric forms. The development of mutant monomeric variants of RFPs has resulted in several novel FPs known as mFruits. Though oxygen is required for maturation of the chromophore, it is known that photobleaching of FPs is oxygen sensitive, and oxygen-free conditions result in improved photostabilities. Therefore, understanding oxygen diffusion pathways in FPs is important for both photostabilites and maturation of the chromophores. In this paper, we use molecular dynamics calculations to investigate the protein barrel fluctuations in mCherry, which is one of the most useful monomeric mFruit variant. We employ implicit ligand sampling to determine oxygen pathways from the bulk solvent into the mCherry chromophore in the interior of the protein. We also show that these pathways can be blocked or altered and barrel fluctuations can be reduced by strategic amino acid substitutions.

Exploring Structural and Optical Properties of Fluorescent Proteins by Squeezing: Modeling High-Pressure Effects on the mStrawberry and mCherry Red Fluorescent Proteins

Molecular dynamics calculations of pressure effects on mStrawberry and mCherry fluorescent proteins are reported. The simulations reveal that mStrawberry has much floppier structure at atmospheric pressure, as evidenced by larger backbone fluctuations and the coexistence of two conformers that differ by Ser146 orientation. Consequently, pressure increase has a larger effect on mStrawberry, making its structure more rigid and reducing the population of one of the conformers. The most significant effect of pressure increase is in the hydrogen-bonding network between the chromophore and the nearby residues. The quantum-mechanics/molecular mechanics calculations of excitation energies in mStrawberry explain the observed blue shift and identify Lys70 as the residue that has the most pronounced effect on the spectra. The results suggest that pressure increase causes an initial increase of fluorescence yield only for relatively floppy fluorescent proteins, whereas the fluorescent proteins that have more rigid structures have quantum yields close to their maximum. The results suggest that a low quantum yield in fluorescent proteins is dynamic in nature and depends on the range of thermal motions of the chromophore and fluctuations in the H-bonding network rather than on their average structure.

Exploring the diffusion of molecular oxygen in the red fluorescent protein mCherry using explicit oxygen molecular dynamics simulations.

The development of fluorescent proteins (FPs) has revolutionized cell biology research. The monomeric variants of red fluorescent proteins (RFPs), known as mFruits, have been especially valuable for tagging and tracking cellular processes in vivo. Determining oxygen diffusion pathways in FPs can be important for improving photostability and for understanding maturation of the chromophore. We use molecular dynamics (MD) calculations to investigate the diffusion of molecular oxygen in one of the most useful monomeric RFPs, mCherry. We describe a pathway that allows oxygen molecules to enter from the solvent and travel through the protein barrel to the chromophore. We calculate the free-energy of an oxygen molecule at points along the path. The pathway contains several oxygen hosting pockets, which are identified by the amino acid residues that form the pocket. We also investigate an RFP variant known to be significantly less photostable than mCherry and find much easier oxygen access in this variant. The results provide a better understanding of the mechanism of molecular oxygen access into the fully folded mCherry protein barrel and provide insight into the photobleaching process in these proteins.

Hydrogen Bond Flexibility Correlates with Stokes Shift in mPlum Variants

Fluorescent proteins have revolutionized molecular biology research and provide a means of tracking subcellular processes with extraordinary spatial and temporal precision. Species with emission beyond 650 nm offer the potential for deeper tissue penetration and lengthened imaging times; however, the origin of their extended Stokes shift is not fully understood. We employed spectrally resolved transient grating spectroscopy and molecular dynamics simulations to investigate the relationship between the flexibility of the chromophore environment and Stokes shift in mPlum. We examined excited state solvation dynamics in a panel of strategic point mutants of residues E16 and I65 proposed to participate in a hydrogen-bonding interaction thought responsible for its red-shifted emission. We observed two characteristic relaxation constants of a few picoseconds and tens of picoseconds that were assigned to survival times of direct and water-mediated hydrogen bonds at the 16-65 position. Moreover, variants of the largest Stokes shift (mPlum, I65V) exhibited significant decay on both time scales, indicating the bathochromic shift correlates with a facile switching between a direct and water-mediated hydrogen bond. This dynamic model underscores the role of environmental flexibility in the mechanism of excited state solvation and provides a template for engineering next-generation red fluorescent proteins.

The Ebola virus protein VP40 hexamer enhances the clustering of PI(4,5)P2 lipids in the plasma membrane

The Ebola virus is a lipid-enveloped virus that obtains its lipid coat from the plasma membrane of the host cell it infects during the budding process. The Ebola virus protein VP40 localizes to the inner leaflet of the plasma membrane and forms the viral matrix, which provides the major structure for the Ebola virus particles. VP40 is initially a dimer that rearranges to a hexameric structure that mediates budding. VP40 hexamers and larger filaments have been shown to be stabilized by PI(4,5)P2 in the plasma membrane inner leaflet. Reduction in the plasma membrane levels of PI(4,5)P2 significantly reduce formation of VP40 oligomers and virus-like particles. We investigated the lipid-protein interactions in VP40 hexamers at the plasma membrane. We quantified lipid-lipid self-clustering by calculating the fractional interaction matrix and found that the VP40 hexamer significantly enhances the PI(4,5)P2 clustering. The radial pair distribution functions suggest a strong interaction between PI(4,5)P2 and the VP40 hexamer. The cationic Lys side chains are found to mediate the PIP2 clustering around the protein, with cholesterol filling the space between the interacting PIP2 molecules. These computational studies support recent experimental data and provide new insights into the mechanisms by which VP40 assembles at the plasma membrane inner leaflet, alters membrane curvature, and forms new virus-like particles.

DNA supercoiling, a critical signal regulating the basal expression of the lac operon in Escherichia coli

Escherichia coli lac repressor (LacI) is a paradigmatic transcriptional factor that controls the expression of lacZYA in the lac operon. This tetrameric protein specifically binds to the O1, O2 and O3 operators of the lac operon and forms a DNA loop to repress transcription from the adjacent lac promoter. In this article, we demonstrate that upon binding to the O1 and O2 operators at their native positions LacI constrains three (−) supercoils within the 401-bp DNA loop of the lac promoter and forms a topological barrier. The stability of LacI-mediated DNA topological barriers is directly proportional to its DNA binding affinity. However, we find that DNA supercoiling modulates the basal expression from the lac operon in E. coli. Our results are consistent with the hypothesis that LacI functions as a topological barrier to constrain free, unconstrained (−) supercoils within the 401-bp DNA loop of the lac promoter. These constrained (−) supercoils enhance LacI’s DNA-binding affinity and thereby the repression of the promoter. Thus, LacI binding is superhelically modulated to control the expression of lacZYA in the lac operon under varying growth conditions.

Structures of the kinetically trapped i-motif DNA intermediates

In the present work, the conformational dynamics and folding pathways of i-motif DNA were studied in solution and in the gas-phase as a function of the solution pH conditions using circular dichroism (CD), photoacoustic calorimetry analysis (PAC), trapped ion mobility spectrometry-mass spectrometry (TIMS-MS), and molecular dynamics (MD). Solution studies showed at thermodynamic equilibrium the existence of a two-state folding mechanism, whereas during the pH = 7.0 → 4.5 transition a fast and slow phase (ΔHfast + ΔHslow = 43 ± 7 kcal mol-1) with a volume change associated with the formation of hemiprotonated cytosine base pairs and concomitant collapse of the i-motif oligonucleotide into a compact conformation were observed. TIMS-MS experiments showed that gas-phase, kinetically trapped i-motif DNA intermediates produced by nanoESI are preserved, with relative abundances depending on the solution pH conditions. In particular, a folded i-motif DNA structure was observed in nanoESI-TIMS-MS for low charge states in both positive and negative ion mode (e.g., z = ±3 to ±5) at low pH conditions. As solution pH increases, the cytosine neutralization leads to the loss of cytosine-cytosine+ (C·CH+) base pairing in the CCC strands and in those conditions we observe partially unfolded i-motif DNA conformations in nanoESI-TIMS-MS for higher charge states (e.g., z = -6 to -9). Collisional induced activation prior to TIMS-MS showed the existence of multiple local free energy minima, associated with the i-motif DNA unfolding at z = -6 charge state. For the first time, candidate gas-phase structures are proposed based on mobility measurements of the i-motif DNA unfolding pathway. Moreover, the inspection of partially unfolded i-motif DNA structures (z = -7 and z = -8 charge states) showed that the presence of inner cations may or may not induce conformational changes in the gas-phase. For example, incorporation of ammonium adducts does not lead to major conformational changes while sodium adducts may lead to the formation of sodium mediated bonds between two negatively charged sides inducing the stabilization towards more compact structures in new local, free energy minima in the gas-phase.

Excluded volume entropic effects on protein unfolding times and intermediary stability

The dynamics of protein folding result from both enthalpic and entropic contributions to the free energy. In this paper we focus on entropic volume exclusion effects. We carry out computer simulations using a model that allows us to independently change the size or biochemical properties of amino acid residues. To determine the importance of excluded volume effects, we investigate the effects of changing the size of side chains on the unfolding dynamics of a model four-helix bundle protein. In addition, we also investigate the effects of changing the thickness of the chains backbone. This has relevance to the behavior of synthetic polymers where the size of the constituent units can be varied. We find that entropic excluded volume effects are crucially important for stabilizing the organized native state relative to the molten globule.