Long‐range molecular dynamics show that inactive forms of Protein Kinase A are more dynamic than active forms

Many protein kinases are characterized by at least two structural forms corresponding to the highest level of activity (active) and low or no activity, (inactive). Further, protein dynamics is an important consideration in understanding the molecular and mechanistic basis of enzyme function. In this work, we use protein kinase A (PKA) as the model system and perform microsecond range molecular dynamics (MD) simulations on six variants which differ from one another in terms of active and inactive form, with or without bound ligands, C‐terminal tail and phosphorylation at the activation loop. We find that the root mean square fluctuations in the MD simulations are generally higher for the inactive forms than the active forms. This difference is statistically significant. The higher dynamics of inactive states has significant contributions from ATP binding loop, catalytic loop, and αG helix. Simulations with and without C‐terminal tail show this differential dynamics as well, with lower dynamics both in the active and inactive forms if C‐terminal tail is present. Similarly, the dynamics associated with the inactive form is higher irrespective of the phosphorylation status of Thr 197. A relatively stable stature of active kinases may be better suited for binding of substrates and detachment of the product. Also, phosphoryl group transfer from ATP to the phosphosite on the substrate requires precise transient coordination of chemical entities from three different molecules, which may be facilitated by the higher stability of the active state.     

Bending of the BST‐2 coiled‐coil during viral budding

BST‐2/tetherin is a human extracellular transmembrane protein that serves as a host defense factor against HIV‐1 and other viruses by inhibiting viral spreading. Structurally, BST‐2 is a homo‐dimeric coiled‐coil that is connected to the host cell membrane by N and C terminal transmembrane anchors. The C‐terminal membrane anchor of BST‐2 is inserted into the budding virus while the N‐terminal membrane anchor remains in the host cell membrane creating a viral tether. The structural mechanism of viral budding and tethering as mediated by BST‐2 is not clear. To more fully describe the mechanism of viral tethering, we created a model of BST‐2 embedded in a membrane and used steered molecular dynamics to simulate the transition from the host cell membrane associated form to the cell‐virus membrane bridging form. We observed that BST‐2 did not transition as a rigid structure, but instead bent at positions with a reduced interface between the helices of the coiled‐coil. The simulations for the human BST‐2 were then compared with simulations on the mouse homolog, which has no apparent weak spots. We observed that the mouse homolog spread the bending across the ectodomain, rather than breaking at discrete points as observed with the human homolog. These simulations support previous biochemical and cellular work suggesting some flexibility in the coiled‐coil is necessary for viral tethering, while also highlighting how subtle changes in protein sequence can influence the dynamics and stability of proteins with overall similar structure.     

Evidence for a matriptase-prostasin proteolytic cascade regulating terminal epidermal differentiation

Recent gene ablation studies in mice have shown that matriptase, a type II transmembrane serine protease, and prostasin, a glycosylphosphatidylinositol-anchored membrane serine protease, are both required for processing of the epidermis-specific polyprotein, profilaggrin, stratum corneum formation, and acquisition of epidermal barrier function. Here we present evidence that matriptase acts upstream of prostasin in a zymogen activation cascade that regulates terminal epidermal differentiation and is required for prostasin zymogen activation. Enzymatic gene trapping of matriptase combined with prostasin immunohistochemistry revealed that matriptase was co-localized with prostasin in transitional layer cells of the epidermis and that the developmental onset of expression of the two membrane proteases was coordinated and correlated with acquisition of epidermal barrier function. Purified soluble matriptase efficiently converted soluble prostasin zymogen to an active two-chain form that formed SDS-stable complexes with the serpin protease nexin-1. Whereas two forms of prostasin with molecular weights corresponding to the prostasin zymogen and active prostasin were present in wild type epidermis, prostasin was exclusively found in the zymogen form in matriptase-deficient epidermis. These data suggest that matriptase, an autoactivating protease, acts upstream from prostasin to initiate a zymogen cascade that is essential for epidermal differentiation.     

Structural and dynamical properties of KTS‐disintegrins: A comparison between Obtustatin and Lebestatin

Obtustatin and Lebestatin are lysine‐threonine‐serine (KTS)‐disintegrins, which are a family of low molecular weight polypeptides present in many viperidae venoms and are potent and specific inhibitors of collagen‐binding integrins. The integrin binding loop, harboring the ²¹KTS²³ motif, and the C‐terminal tail are known to be responsible for the selective binding to the α1β1 integrin. Despite a very high sequence homology (only two mutations are present in Lebestatin relative to Obtustatin, namely R24L and S38L), Lebestatin exhibits a higher inhibitory effect than Obtustatin on cell adhesion and cell migration to collagens I and IV. Here we show, by means of molecular dynamics simulations of the two polypeptides in aqueous solution, that Lebestatin possesses a higher flexibility of the C‐terminal tail and a greater solvent accessibility of the integrin binding loop than Obtustatin. It may be hypothesized that these properties may contribute to the higher binding‐affinity of Lebestatin to its biological partner. © 2012 Wiley Periodicals, Inc.     

The influenza fusion peptide promotes lipid polar head intrusion through hydrogen bonding with phosphates and N‐terminal membrane insertion depth

Influenza infection requires fusion between the virus envelope and a host cell endosomal membrane. The influenza hemagglutinin fusion peptide (FP) is essential to viral membrane fusion. It was recently proposed that FPs would fuse membranes by increasing lipid tail protrusion, a membrane fusion transition state. The details of how FPs induce lipid tail protrusion, however, remain to be elucidated. To decipher the molecular mechanism by which FPs promote lipid tail protrusion, we performed molecular dynamics simulations of the wild‐type (WT) FP, fusogenic mutant F9A, and nonfusogenic mutant W14A in model bilayers. This article presents the peptide–lipid interaction responsible for lipid tail protrusion and a related lipid perturbation, polar head intrusion, where polar heads are sunk under the membrane surface. The backbone amides from the four N‐terminal peptide residues, deeply inserted in the membrane, promoted both perturbations through H bonding with lipid phosphates. Polar head intrusion correlated with peptides N‐terminal insertion depth and activity: the N‐termini of WT and F9A were inserted deeper into the membrane than nonfusogenic W14A. Based on these results, we propose that FP‐induced polar head intrusion would complement lipid tail protrusion in catalyzing membrane fusion by reducing repulsions between juxtaposed membranes headgroups. The presented model provides a framework for further research on membrane fusion and influenza antivirals. Proteins 2014; 82:2118–2127. © 2014 Wiley Periodicals, Inc.     

Computational delineation of the catalytic step of a high‐fidelity DNA polymerase

The Bacillus fragment, belonging to a class of high‐fidelity polymerases, demonstrates high processivity (adding ～115 bases per DNA binding event) and exceptional accuracy (1 error in 10⁶ nucleotide incorporations) during DNA replication. We present analysis of structural rearrangements and energetics just before and during the chemical step (phosphodiester bond formation) using a combination of classical molecular dynamics, mixed quantum mechanics molecular mechanics simulations, and free energy computations. We find that the reaction is associative, proceeding via the two‐metal‐ion mechanism, and requiring the proton on the terminal primer O3′ to transfer to the pyrophosphate tail of the incoming nucleotide before the formation of the pentacovalent transition state. Different protonation states for key active site residues direct the system to alternative pathways of catalysis and we estimate a free energy barrier of ～12 kcal/mol for the chemical step. We propose that the protonation of a highly conserved catalytic aspartic acid residue is essential for the high processivity demonstrated by the enzyme and suggest that global motions could be part of the reaction free energy landscape.     

Conformational flexibility of the complete catalytic domain of Cdc25B phosphatases

Cdc25B phosphatases are involved in cell cycle checkpoints and have become a possible target for developing new anticancer drugs. A more rational design of Cdc25B ligands would benefit from detailed knowledge of its tertiary structure. The conformational flexibility of the C‐terminal region of the Cdc25B catalytic domain has been debated recently and suggested to play an important structural role. Here, a combination of experimental NMR measurements and molecular dynamics simulations for the complete catalytic domain of the Cdc25B phosphatase is presented. The stability of the C‐terminal α‐helix is confirmed, but the last 20 residues in the complete catalytic domain are very flexible, partially occlude the active site and may establish transient contacts with the protein core. This flexibility in the C‐terminal tail may modulate the molecular recognition of natural substrates and competitive inhibitors by Cdc25B. Proteins 2016; 84:1567–1575. © 2016 Wiley Periodicals, Inc.     

Activation of the easter zymogen is regulated by five other genes to define dorsal-ventral polarity in the Drosophila embryo.

The product of the Drosophila easter gene, a member of the trypsin family of serine proteases, must be more active ventrally than dorsally to promote normal embryonic polarity. The majority of the easter protein in the embryo is present in the unprocessed zymogen form and appears to be evenly distributed in the extracellular space, indicating that the asymmetric activity of wild-type easter must arise post-translationally. A dominant mutant form of easter that does not require cleavage of the zymogen for activity (ea delta N) is active both dorsally and ventrally. The ea delta N mutant bypasses the requirement for five other maternal effect genes, indicating that these five genes exert their effects on dorsal-ventral patterning solely by controlling the activation of the easter zymogen. We propose that dorsal-ventral asymmetry is initiated by a ventrally-localized molecule in the vitelline membrane that nucleates an easter zymogen activation complex, leading to the production of ventrally active easter enzyme.     

Molecular insight into conformational transformation of human glucokinase: conventional and targeted molecular dynamics

Abstract

Glucokinase (GK) plays a key role in the regulation of hepatic glucose metabolism. An unusual mechanism of positive cooperativity of monomeric GK containing only a single binding site for glucose is very interesting and still unclear. The activation process of GK is associated with a large-scale conformational change from the inactive to the active state. Here, conventional and targeted molecular dynamics simulations were used to study the conformational dynamics of GK in the stable configurations and in the transition from active to inactive state. Three phases of the structural reorganization of GK were detected. The first step is a transformation of GK from the active state to the intermediate structure, where the cleft between the domains is open, but alpha helix 13 is still inside the small domain. From this point, there are two alternative paths. One path leads to the inactive state through the release of helix 13 from the inside of small domain to the outside. Other path goes back to the active state. Simulation results reveal the critical role of helix 13 in the transformation of GK from the open state to inactive one and the influence of the loop 2 on the protein transformation between the open and the closed active states. Principal component analysis and covariance matrix analysis were carried out to analyze the dynamics of protein. Importance of hydrogen bonds in the stability of the closed conformation is shown. Overall, our simulations provide new information about the dynamics of GK and its structural transformation.

Communicated by Ramaswamy H. Sarma     

Pepsinogen-like activation intermediate of plasmepsin II revealed by molecular dynamics analysis

Plasmepsins are pharmaceutically relevant aspartic proteases involved in haemoglobin degradation by the malaria causing parasites Plasmodium spp. They are translated as inactive proenzymes, with an elongated prosegment. On prosegment cleavage, plasmepsins undergo a series of hitherto unresolved conformational changes before becoming active. Here, the flexibility of plasmepsin and proplasmepsin and the activation process are investigated by multiple explicit water molecular dynamics simulations. The large N-terminal displacement and the interdomain shift from the proenzyme structure to active plasmepsin are promoted by essential dynamics sampling. An intermediate, stabilized by electrostatic interactions between the catalytic dyad and the N-terminus of mature plasmepsin, is observed along all activation trajectories. Notably, the stabilizing interactions in the activation intermediate of plasmepsin are similar to those in the X-ray structure of pepsinogen. In particular, the catalytic aspartates act as hydrogen bond acceptors for the N-terminal amino group and the Ser2 hydroxyl in plasmepsin, and the side chains of Lys36pro and Tyr9 in pepsinogen. The simulation results are used to suggest in vitro experiments to test the conformational transitions involved in the maturation of plasmepsin, and design small-molecule inhibitors.     

Conformational dynamics of the EGFR kinase domain reveals structural features involved in activation

The epidermal growth factor receptor (EGFR) has been the focus of intensive studies because of its importance in cancer research. Thus, a broader understanding of the molecular mechanism of activation of the EGFR kinase will have profound significance for the development of novel therapeutics. Numerous crystal structures of EGFR kinase, including the structure of the activating‐kinase dimer, have provided snapshots of the specific pathway. Herein, we performed unrestrained‐, as well as targeted‐molecular dynamics simulations based on these data, to gain further insight into the conformational changes responsible for activation. Comparison of the monomer‐ versus activating‐EGFR‐dimer simulations indicates that the dimerization is stabilizing structural elements associated with the activated state and predicts new salt‐bridge interactions involving activation‐loop residues that may also be associated with that state. Targeted molecular dynamics simulations of the inactive‐to‐active EGFR transition, as well as the reverse pathway, confirm the formation of conserved structural features of functional importance for the activity or stabilization of either conformation. Interestingly, simulations of the L834R mutant, which is associated with cancer, suggest that the structural basis of the activation induced by that mutation might be the ability of the mutated R834 residue to consecutively form salt bridges with neighboring acidic residues and cause destabilization of a hydrophobic cluster in the inactive state. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

The N-terminal Arg Residue Is Essential for Autocatalytic Activation of a Lipopolysaccharide-responsive Protease Zymogen

Factor C, a serine protease zymogen involved in innate immune responses in horseshoe crabs, is known to be autocatalytically activated on the surface of bacterial lipopolysaccharides, but the molecular mechanism of this activation remains unknown. In this study, we show that wild-type factor C expressed in HEK293S cells exhibits a lipopolysaccharide-induced activity equivalent to that of native factor C. Analysis of the N-terminal addition, deletion, or substitution mutants shows that the N-terminal Arg residue and the distance between the N terminus and the tripartite of lipopolysaccharide-binding site are essential factors for autocatalytic activation, and that the positive charge of the N terminus may interact with an acidic amino acid(s) of the molecule to convert the zymogen into an active form. Chemical cross-linking experiments indicate that the N terminus is required to form a complex of the factor C molecules in a sufficiently close vicinity to be chemically cross-linked on the surface of lipopolysaccharides. We propose a molecular mechanism of the autocatalytic activation of the protease zymogen on lipopolysaccharides functioning as a platform to induce specific protein-protein interaction between the factor C molecules.     

Is there nascent structure in the intrinsically disordered region of troponin I?

In striated muscle, the binding of calcium to troponin C (TnC) results in the removal of the C‐terminal region of the inhibitory protein troponin I (TnI) from actin. While structural studies of the muscle system have been successful in determining the overall organization of most of the components involved in force generation at the atomic level, the structure and dynamics of the C‐terminal region of TnI remains controversial. This domain of TnI is highly flexible, and it has been proposed that this intrinsically disordered region (IDR) regulates contraction via a “fly‐casting” mechanism. Different structures have been presented for this region using different methodologies: a single α‐helix, a “mobile domain” containing a small β‐sheet, an unstructured region, and a two helix segment. To investigate whether this IDR has in fact any nascent structure, we have constructed a skeletal TnC‐TnI chimera that contains the N‐domain of TnC (1–90), a short linker (GGAGG), and the C‐terminal region of TnI (97–182) and have acquired ¹⁵N NMR relaxation data for this chimera. We compare the experimental relaxation parameters with those calculated from molecular dynamic simulations using four models based upon the structural studies. Our experimental results suggest that the C‐terminal region of TnI does not contain any defined secondary structure, supporting the “fly‐casting” mechanism. We interpret the presence of a “plateau” in the ¹⁵N NMR relaxation data as being an intrinsic property of IDRs. We also identified a more rigid adjacent region of TnI that has implications for muscle performance under ischemic conditions. Proteins 2011. © 2011 Wiley‐Liss, Inc.     

Structural analysis of the N‐terminal fragment of the antiangiogenic protein endostatin: A molecular dynamics study

Endostatin is a potent antiangiogenic protein derived from the noncollagenous domain 1 (NC1) of collagen XVIII. The mechanism by which endostatin exerts its antiangiogenic effect is still incompletely understood. It has been shown that the 27 amino acid N‐terminal fragment of murine endostatin has antitumor, antimigration, and antipermeability activities comparable to the full soluble protein. To understand how this peptide can exert such elaborate function, we performed structural analysis using molecular dynamics to evaluate the behavior of this fragment in aqueous environment. Here, we show that the N‐terminal peptide of murine endostatin is able to assume a well‐defined structure, folding into a zinc‐dependent β‐hairpin conformation. Analyzing the folding mechanism, we were able to understand why the N‐terminal peptide of human endostatin with the same length failed to acquire a stable conformation. Conversely, we were able to predict the successful folding of the R4Q mutant and of a shorter form of the human peptide with 25 residues. Finally, we show that the β‐hairpin conformation assumed by the zinc‐bound peptide of murine endostatin has a high structural similarity with fragments of another family of angiogenesis inhibitors: the integrin‐binding portion of the NC1 domain of collagen IV. Indeed, our docking simulations show that arresten, canstatin, and the endostatin peptide bind to the same spot of αVβ3 integrin, suggesting similar interactions via a common binding site on this receptor. Proteins 2011;. © 2011 Wiley‐Liss, Inc.     

Structure-function analysis of grass clip serine protease involved in Drosophila Toll pathway activation

Grass is a clip domain serine protease (SP) involved in a proteolytic cascade triggering the Toll pathway activation of Drosophila during an immune response. Epistasic studies position it downstream of the apical protease ModSP and upstream of the terminal protease Spaetzle-processing enzyme. Here, we report the crystal structure of Grass zymogen. We found that Grass displays a rather deep active site cleft comparable with that of proteases of coagulation and complement cascades. A key distinctive feature is the presence of an additional loop (75-loop) in the proximity of the activation site localized on a protruding loop. All biochemical attempts to hydrolyze the activation site of Grass failed, strongly suggesting restricted access to this region. The 75-loop is thus proposed to constitute an original mechanism to prevent spontaneous activation. A comparison of Grass with clip serine proteases of known function involved in analogous proteolytic cascades allowed us to define two groups, according to the presence of the 75-loop and the conformation of the clip domain. One group (devoid of the 75-loop) contains penultimate proteases whereas the other contains terminal proteases. Using this classification, Grass appears to be a terminal protease. This result is evaluated according to the genetic data documenting Grass function.     

A true autoactivating enzyme. Structural insight into mannose-binding lectin-associated serine protease-2 activations

Few reports have described in detail a true autoactivation process, where no extrinsic cleavage factors are required to initiate the autoactivation of a zymogen. Herein, we provide structural and mechanistic insight into the autoactivation of a multidomain serine protease: mannose-binding lectin-associated serine protease-2 (MASP-2), the first enzymatic component in the lectin pathway of complement activation. We characterized the proenzyme form of a MASP-2 catalytic fragment encompassing its C-terminal three domains and solved its crystal structure at 2.4 A resolution. Surprisingly, zymogen MASP-2 is capable of cleaving its natural substrate C4, with an efficiency about 10% that of active MASP-2. Comparison of the zymogen and active structures of MASP-2 reveals that, in addition to the activation domain, other loops of the serine protease domain undergo significant conformational changes. This additional flexibility could play a key role in the transition of zymogen MASP-2 into a proteolytically active form. Based on the three-dimensional structures of proenzyme and active MASP-2 catalytic fragments, we present model for the active zymogen MASP-2 complex and propose a mechanism for the autoactivation process.     

Terminal sialic acids on CD44 N‐glycans can block hyaluronan binding by forming competing intramolecular contacts with arginine sidechains

Specific sugar residues and their linkages form the basis of molecular recognition for interactions of glycoproteins with other biomolecules. Seemingly small changes, like the addition of a single monosaccharide in the covalently attached glycan component of glycoproteins, can greatly affect these interactions. For instance, the sialic acid capping of glycans affects protein‐ligand binding involved in cell–cell and cell–matrix interactions. CD44 is a single‐pass transmembrane glycoprotein whose binding with its carbohydrate ligand hyaluronan (HA), an extracellular matrix component, mediates processes such as leukocyte homing, cell adhesion, and tumor metastasis. This binding is highly regulated by glycosylation of the N‐terminal extracellular hyaluronan‐binding domain (HABD); specifically, sialic acid capped N‐glycans of HABD inhibit ligand binding. However, the molecular mechanism behind this sialic acid mediated regulation has remained unknown. Two of the five N‐glycosyation sites of HABD have been previously identified as having the greatest inhibitory effect on HA binding, but only if the glycans contain terminal sialic acid residues. These two sites, Asn25 and Asn120, were chosen for in silico glycosylation in this study. Here, from extensive standard molecular dynamics simulations and biased simulations, we propose a molecular mechanism for this behavior based on spontaneously‐formed charge‐paired hydrogen bonding interactions between the negatively‐charged sialic acid residues and positively‐charged Arg sidechains known to be critically important for binding to HA, which itself is negatively charged. Such intramolecular hydrogen bonds would preclude associations critical to hyaluronan binding. This observation suggests how CD44 and related glycoprotein binding is regulated by sialylation as cellular environments fluctuate. Proteins 2014; 82:3079–3089. © 2014 Wiley Periodicals, Inc.     

Molecular Basis of Enhanced Activity in Factor VIIa-Trypsin Variants Conveys Insights into Tissue Factor-mediated Allosteric Regulation of Factor VIIa Activity

The complex of coagulation factor VIIa (FVIIa), a trypsin-like serine protease, and membrane-bound tissue factor (TF) initiates blood coagulation upon vascular injury. Binding of TF to FVIIa promotes allosteric conformational changes in the FVIIa protease domain and improves its catalytic properties. Extensive studies have revealed two putative pathways for this allosteric communication. Here we provide further details of this allosteric communication by investigating FVIIa loop swap variants containing the 170 loop of trypsin that display TF-independent enhanced activity. Using x-ray crystallography, we show that the introduced 170 loop from trypsin directly interacts with the FVIIa active site, stabilizing segment 215–217 and activation loop 3, leading to enhanced activity. Molecular dynamics simulations and novel fluorescence quenching studies support that segment 215–217 conformation is pivotal to the enhanced activity of the FVIIa variants. We speculate that the allosteric regulation of FVIIa activity by TF binding follows a similar path in conjunction with protease domain N terminus insertion, suggesting a more complete molecular basis of TF-mediated allosteric enhancement of FVIIa activity.     

Modeling of substrate and inhibitory complexes of histidine-aspartic protease

A three-dimensional structure of histo-aspartic protease (HAP), a pepsin-like enzyme from the causative agent of malaria Plasmodium falciparum, is suggested on the basis of homologous modeling followed by equilibration by the method of molecular dynamics. The presence of a His residue in the catalytic site instead of an Asp residue, which is characteristic of pepsin-like enzymes, and replacement of some other conserved residues in the active site make it possible for the enzyme to function by the covalent mechanism inherent in serine proteases. The detailed structures of HAP complexes with pepstatin, a noncovalent inhibitor of aspartic proteases, and phenylmethylsulfonyl fluoride, a covalent inhibitor of serine proteases, as well as with a pentapeptide substrate are discussed.