Dynamic Cross-correlation Analysis Research Articles

The binding pattern of ferric iron and iron-binding protein in Botrytis cinerea

Iron-binding protein (Ibp) has protective effect on pathogen exposed to H2O2 in defense response of plants. Ibp in Botrytis cinerea (BcIbp) is related to its virulence. Bcibp mutation lead to virulence deficiencies in B. cinerea. BcIbp is involved in the Fe3+ homeostasis regulation. Recognition the binding site and binding pattern of ferric iron and iron-binding protein in B. cinerea are vital to understand its function. In this study, molecular dynamics (MD) simulations, gaussian accelerated molecular dynamics (GaMD) simulations, dynamic cross correlation analysis and quantum chemical energy calculation were used to explore binding pattern of ferric iron. MD results showed that the C-terminal region had little effect on the stability of residues in the Fe3+-binding pocket. Energy calculations suggested the most likely coordination pattern for ferric iron in iron-binding protein. These results will help to understand the binding of ferric iron to iron-binding protein and provide new ideas for regulating the virulence of B. cinerea.

Discovery of small molecule entry inhibitors targeting the linoleic acid binding pocket of SARS-CoV-2 spike protein

Spike glycoprotein has a significant role in the entry of SARS-CoV-2 to host cells, which makes it a potential drug target. Continued accumulation of non-synonymous mutations in the receptor binding domain of spike protein poses great challenges in identifying antiviral drugs targeting this protein. This study aims to identify potential entry inhibitors of SARS-CoV-2 using virtual screening and molecular dynamics (MD) simulations from three distinct chemical libraries including Pandemic Response Box, Drugbank and DrugCentral, comprising 6971 small molecules. The molecules were screened against a binding pocket identified in the receptor-binding domain (RBD) region of the spike protein which is known as the linoleic acid binding pocket, a highly conserved motif among several SARS-CoV-2 variants. Through virtual screening and binding free energy calculations, we identified four top-scoring compounds, MMV1579787 ([2-Oxo-2-[2-(3-phenoxyphenyl)ethylamino]ethyl]phosphonic acid), Tretinoin, MMV1633963 ((2E,4E)-5-[3-(3,5-dichlorophenoxy)phenyl]penta-2,4-dienoic acid) and Polydatin, which were previously reported to have antibacterial, antifungal or antiviral properties. These molecules showed stable binding on MD simulations over 100 ns and maintained stable interactions with TYR365, PHE338, PHE342, PHE377, TYR369, PHE374 and LEU368 of the spike protein RBD that are found to be conserved among SARS-CoV-2 variants. Our findings were further validated with free energy landscape, principal component analysis and dynamic cross-correlation analysis. Our in silico analysis of binding mode and MD simulation analyses suggest that the identified compounds may impede viral entrance by interacting with the linoleic acid binding site of the spike protein of SARS-CoV-2 regardless of its variants, and they thus demand for further in vitro and in vivo research. Communicated by Ramaswamy H. Sarma

The effect of phosphorylation on the conformational dynamics and allostery of the association of death-associated protein kinase with calmodulin

Protein phosphorylation plays an important role in the signal transduction and is capable of regulation of cell activity. The death-associated protein kinase 1 (DAPK1), as a Ser/Thr kinase, interacts with calmodulin (CaM) to regulate apoptotic and autophagic signaling. Autophosphorylation of DAPK1 at Ser308 located at the autoregulatory domain (ARD) blocks CaM binding and inhibits kinase catalytic activity. However, the mechanism underlying the influence of Ser308 phosphorylation (pS308) on the DAPK1 activity remains unclear. Here, we performed multiple, microsecond length molecular dynamics (MD) simulations, the molecular mechanics generalized Born/surface area (MM-GBSA) binding free energy calculations, principal component analysis, and dynamic cross-correlation analysis to unravel the conformational dynamics and allostery of the DAPK1 − CaM interaction triggered by the pS308 at the ARD. MD simulations showed that pS308 affected the conformational stability of the DAPK1 − CaM complex. Further energetic and structural exploration revealed that pS308 weakened the association of the phosphorylated DAPK1 to CaM, which lowered the susceptibility of DAPK1 to be activated by CaM. This result can provide mechanistic insights into the molecular underpinning through which the DAPK1 kinase activity is modulated by the auto-phosphorylation. Communicated by Ramaswamy H. Sarma

An in-silico approach to identify bioactive phytochemicals from Houttuynia cordata Thunb. As potential inhibitors of human glutathione reductase

Cellular infections are central to the etiology of various diseases, notably cancer and malaria. Counteracting cellular oxidative stress via the inhibition of glutathione reductase (GR) has emerged as a promising therapeutic strategy. Houttuynia cordata, a medicinal plant known for its potent antioxidant properties, has been the focus of our investigation. In this study, we conducted comprehensive in silico analyses involving the phytochemical constituents of H. cordata to identify potential natural GR inhibitors. Our methodological approach encompassed multiple in silico techniques, including molecular docking, molecular dynamics simulations, MMPBSA analysis, and dynamic cross-correlation analysis. Out of 13 docked phytochemicals, Quercetin, Quercitrin, and Sesamin emerged as particularly noteworthy due to their exceptional binding affinities for GR. Notably, our investigation demonstrated that Quercetin and Sesamin exhibited promising outcomes compared to the well-established pharmaceutical agent N-acetylcysteine (NAC). Molecular dynamics analyses provided insights into the ability of these phytochemicals to induce structural compaction and stabilization of the GR protein, as evidenced by changes in radius of gyration and solvent-accessible surface area. Moreover, MMPBSA analysis highlighted the crucial roles of specific residues, namely Gly27, Gly28, Ser51, His52, and Val61, in mediating essential interactions with these phytochemicals. Furthermore, an assessment of Absorption, Distribution, Metabolism, Excretion, and Toxicity (ADME-Tox) profiles underscored the favourable drug-like attributes of these phytochemicals. Thus, the current findings underscore the immense potential of Houttuynia cordata phytochemicals as potent antioxidants with the capacity to combat a spectrum of maladies, including malaria and cancer. This study not only unveils novel therapeutic avenues but also underscores the distinctive outcomes and paramount significance of harnessing H. cordata phytochemicals for their efficacious antioxidant properties. Communicated by Ramaswamy H. Sarma

Molecular docking and molecular dynamics simulations reveal the clinical resistance of the gatekeeper mutation V564F of FGFR2 against Infigratinib

Fibroblast growth factor receptor 2 (FGFR2), as a transmembrane receptor tyrosine kinase, is implicated in a plethora of human cancers, including intrahepatic cholangiocarcinomas, breast cancers, and non-small cell lung cancer. The clinically relevant V564F gatekeeper mutation conferred resistance to current FGFR2 drug − Infigratinib. In this study, the protein − ligand interactions between FGFR2 kinase domain (wild-type and V564F) and Infigratinib were compared through an integrated computational method. The multiple, large-scale molecular dynamics (MD) simulations, together with dynamic cross-correlation analysis and binding free energy calculations suggested that the resistant mutation may not trigger the conformational changes of the FGFR2 kinase domain. The simulation results also indicated that the driving force to decrease the binding affinity of Infigratinib to the FGFR2 V564F variant derived from the difference in the protein − ligand hydrogen bonding interactions. Moreover, the per-residue free energy decomposition analysis revealed that the reduced contributions from several residues in the ATP-binding site of FGFR2, especially Glu565 and Ala567 located at the kinase hinge domain, were the key residues responsible for the loss of binding affinity of Infigratinib to the V564F variant. This study may offer useful information for the design of novel selective kinase inhibitors targeting FGFR2.

Assessing the impact of the Russia–Ukraine war on energy prices: A dynamic cross-correlation analysis

Assessing the impact of the Russia–Ukraine war on energy prices: A dynamic cross-correlation analysis

Modeling the role of charged residues in thermophilic proteins by rotamer and dynamic cross correlation analysis.

Discerning the determinants of protein thermostability is very important both from the theoretical and applied perspective. Different lines of evidence seem to indicate that a dynamical network of salt bridges/charged residues plays a fundamental role in the thermostability of enzymes. In this work, we applied measures of dynamic variance, like the Gini coefficients, Kullback-Leibler (KL) divergence and dynamic cross correlation (DCC) coefficients to compare the behavior of 3 pairs of homologous proteins from the thermophilic bacterium Thermus thermophilus and mesophilic Escherichia coli. Molecular dynamic (MD) simulations of these proteins were performed at 303K and 363K. In the characterization of their side chain rotamer distributions, the corresponding Gini coefficients and KL-divergence both revealed significant correlations with temperature. Similarly, a DCC analysis revealed a higher trend to de-correlate the movement of charged residues at higher temperatures in the thermophilic proteins, when compared with their mesophilic homologues. These results highlight the importance of dynamic electrostatic network interactions for the thermostability of enzymes.

Computational investigation on the effect of the lysine 2-hydroxyisobutyrylation on argininosuccinate synthetase 1 conformational dynamics in Botrytis cinerea.

Lysine 2-hydroxyisobutyrylation (Khib) is a newly discovered post-translational modification in recent years, which has been identified in several species and is associated with diverse cellular functions. Botrytis cinerea, as a broad host pathogen, is very destructive and causes serious losses to agricultural economy. Argininosuccinate synthetase (ASS, citrulline-aspartate ligase) is the rate-limiting enzyme in the catalytic arginine synthesis pathway. Arginine deficiency can affect the growth of Botrytis cinerea. The Khib site Lys120 was found in functional domain of argininosuccinate synthetase 1 from Botrytis cinerea (Bcass1), which is located in conserved loop. It is worth exploring how K120hib affects the conformation of Bcass1. In this study, molecular dynamics (MD) simulations, binding free energy calculation, principal component analysis (PCA), and dynamic cross-correlation analysis were used to explore the influence of K120hib on the conformation of Bcass1. The increase of root-mean-square fluctuation (RMSF) value of related residues and PCA results suggests that K120hib increases the flexibility of some regions of Bcass1. Moreover, K120hib weakens the binding free energy between Bcass1 and the two substrates. These results will help to understand the effects of K120hib on Bcass1 and provide new ideas for regulating the pathogenicity of Botrytis cinerea.

Can Second Coordination Sphere and Long-Range Interactions Modulate Hydrogen Atom Transfer in a Non-Heme Fe(II)-Dependent Histone Demethylase?

Fe(II)-dependent oxygenases employ hydrogen atom transfer (HAT) to produce a myriad of products. Understanding how such enzymes use dynamic processes beyond the immediate vicinity of the active site to control the selectivity and efficiency of HAT is important for metalloenzyme engineering; however, obtaining such knowledge by experiments is challenging. This study develops a computational framework for identifying second coordination sphere (SCS) and especially long-range (LR) residues relevant for catalysis through dynamic cross-correlation analysis (DCCA) using the human histone demethylase PHF8 (KDM7B) as a model oxygenase. Furthermore, the study explores the mechanistic pathways of influence of the SCS and LR residues on the HAT reaction. To demonstrate the plausibility of the approach, we investigated the effect of a PHF8 F279S clinical mutation associated with X-linked intellectual disability, which has been experimentally shown to ablate PHF8-catalyzed demethylation. In agreement, the molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) studies showed a change in the H31–14K9me2 substrate orientation and an increased HAT barrier. We systematically analyzed the pathways by which the identified SCS and LR residues may influence HAT by exploring changes in H3K9me2 substrate orientation, interdomain correlated motions, HAT transition state stabilization, reaction energetics, electron transfer mechanism, and alterations in the intrinsic electric field of PHF8. Importantly, SCS and LR variations decrease key motions of α9−α12 of the JmjC domain toward the Fe(IV)-center that are associated with tighter binding of the H31–14K9me2 substrate. SCS and LR residues alter the intrinsic electric field of the enzyme along the reaction coordinate and change the individual energetic contributions of residues toward TS stabilization. The overall results suggest that DCCA can indeed identify non-active-site residues relevant for catalysis. The substitutions of such dynamically correlated residues might be used as a tool to tune HAT in non-heme Fe(II)- and 2OG-dependent enzymes.

Insights into the Allosteric Effect of SENP1 Q597A Mutation on the Hydrolytic Reaction of SUMO1 via an Integrated Computational Study.

Small ubiquitin-related modifier (SUMO)-specific protease 1 (SENP1) is a cysteine protease that catalyzes the cleavage of the C-terminus of SUMO1 for the processing of SUMO precursors and deSUMOylation of target proteins. SENP1 is considered to be a promising target for the treatment of hepatocellular carcinoma (HCC) and prostate cancer. SENP1 Gln597 is located at the unstructured loop connecting the helices α4 to α5. The Q597A mutation of SENP1 allosterically disrupts the hydrolytic reaction of SUMO1 through an unknown mechanism. Here, extensive multiple replicates of microsecond molecular dynamics (MD) simulations, coupled with principal component analysis, dynamic cross-correlation analysis, community network analysis, and binding free energy calculations, were performed to elucidate the detailed mechanism. Our MD simulations showed that the Q597A mutation induced marked dynamic conformational changes in SENP1, especially in the unstructured loop connecting the helices α4 to α5 which the mutation site occupies. Moreover, the Q597A mutation caused conformational changes to catalytic Cys603 and His533 at the active site, which might impair the catalytic activity of SENP1 in processing SUMO1. Moreover, binding free energy calculations revealed that the Q597A mutation had a minor effect on the binding affinity of SUMO1 to SENP1. Together, these results may broaden our understanding of the allosteric modulation of the SENP1−SUMO1 complex.

Identification and functional assays of single-nucleotide variants of opsins genes in melanocytic tumors.

Epidermal melanocytes sense solar light via the opsin-coupled signaling pathway which is involved in a range of biological functions, including regulating pigmentation, proliferation, apoptosis, and tumorigenesis. However, it remains unclear whether there are genetic variants within these opsins that affect opsin protein structure and function, and further melanocyte biological behaviors. Here, we examined single-nucleotide variants (SNVs) of five opsin (RGR, OPN1SW, OPN2, OPN4, and OPN5) genes in MM (malignant melanoma; n=76) and MN (melanocytic nevi; n=157), using next-generation sequencing. The effects of these pathogenic single-nucleotide variants (SNVs) on opsin structure and function were further investigated using molecular dynamics (MD) simulations, dynamic cross-correlation (DCC), and site-directed mutagenesis. In total, 107 SNV variants were identified. Of these variants, 14 nonsynonymous SNVs (nsSNVs) of opsin genes were detected, including three mutations in the RGR gene, three mutations in the OPN1SW gene, two mutations in the OPN2 gene, and six mutations in the OPN4 gene. The effect of these missense mutations on opsin function was then assessed using eight prediction tools to estimate the potential impact of an amino acid substitution. The impact of each nsSNV was investigated using MD simulations and DCC analysis. Furthermore, we performed in vitro fluorescence calcium imaging to assess the functional properties of nsSNV proteins using a site-directed mutagenesis method. Taken together, these results revealed that p.A103V (RGR), p.T167I (RGR), p.G141S (OPN1SW), p.R144C (OPN1SW), and p.S231F (OPN4) had more deleterious effects on protein structure and function among the 14 nsSNVs. Opsin gene alterations showed the low frequency of missense mutations in melanocytic tumors, and although rare, some mutations in these opsin genes disrupt the canonical function of opsin. Our findings provide new insight into the role of opsin variants in the loss of function.

Effect of Sulfation on the Conformational Dynamics of Dermatan Sulfate Glycosaminoglycan: A Gaussian Accelerated Molecular Dynamics Study.

Glycosaminoglycans (GAGs) are anionic biopolymers present on cell surfaces as a part of proteoglycans. The biological activities of GAGs depend on the sulfation pattern. In our study, we have considered three octadecasaccharide dermatan sulfate (DS) chains with increasing order of sulfation (dp6s, dp7s, and dp12s) to illuminate the role of sulfation on the GAG units and its chain conformation through 10 μs-long Gaussian accelerated molecular dynamics simulations. DS is composed of repeating disaccharide units of iduronic acid (IdoA) and N-acetylgalactosamine (N-GalNAc). Here, N-GalNAc is linked to IdoA via β(1-4), while IdoA is linked to N-GalNAc through α(1-3). With the increase in sulfation, the DS structure becomes more rigid and linear, as is evident from the distribution of root-mean-square deviations (RMSDs) and end-to-end distances. The tetrasaccharide linker region of the main chain shows a rigid conformation in terms of the glycosidic linkage. We have observed that upon sulfation (i.e., dp12s), the ring flip between two chair forms vanished for IdoA. The dynamic cross-correlation analysis reveals that the anticorrelation motions in dp12s are reduced significantly compared to dp6s or dp7s. An increase in sulfation generates relatively more stable hydrogen-bond networks, including water bridging with the neighboring monosaccharides. Despite the favorable linear structures of the GAG chains, our study also predicts few significant bendings related to the different puckering states, which may play a notable role in the function of the DS. The relation between the global conformation with the micro-level parameters such as puckering and water-mediated hydrogen bonds shapes the overall conformational space of GAGs. Overall, atomistic details of the DS chain provided in this study will help understand their functional and mechanical roles, besides developing new biomaterials.

An insight into the inhibitory mechanism of phytochemicals and FDA-approved drugs on the ACE2-Spike complex of SARS-CoV-2 using computational methods.

The S-glycoprotein (Spike) of the SARS-CoV-2 forms a complex with the human transmembrane protein angiotensin-converting enzyme 2 (ACE2) during infection. It forms the first line of contact with the human cell. The FDA-approved drugs and phytochemicals from Indian medicinal plants were explored. Molecular docking and simulations of these molecules targeting the ACE2–Spike complex were performed. Rutin DAB10 and Swertiapuniside were obtained as the top-scored drugs as per the docking protocol. The MD simulations of ligand-free, Rutin DAB10-bound, and Swertiapuniside-bound ACE2–Spike complex revealed abrogation of the hydrogen bonding network between the two proteins. The principal component and dynamic cross-correlation analysis pointed out conformational changes in both the proteins unique to the ligand-bound systems. The interface residues, His34, and Lys353 from ACE2 and Arg403, and Tyr495 from the Spike protein formed significant strong interactions with the ligand molecules, inferring the inhibition of ACE2–Spike complex. Few novel interactions specific to Rutin-DAB10 and Swertiapuniside were also identified. The conformational flexibility of the drug-binding pocket was captured using the RMSD-based clustering of the ligand-free simulations. Ensemble docking was performed wherein the FDA-approved database and phytochemical dataset were docked on each of the cluster representatives of the ACE2–Spike. The phytochemicals identified belonged to Withania somnifera, Swertia chirayita, Tinospora cordifolia and Rutin DAB10, fulvestrant, elbasvir from FDA.Supplementary InformationThe online version contains supplementary material available at 10.1007/s11696-021-01680-1.

Conserved Conformational Hierarchy across Functionally Divergent Glycosyltransferases of the GT-B Structural Superfamily as Determined from Microsecond Molecular Dynamics.

It has long been understood that some proteins undergo conformational transitions en route to the Michaelis Complex to allow chemistry. Examination of crystal structures of glycosyltransferase enzymes in the GT-B structural class reveals that the presence of ligand in the active site triggers an open-to-closed conformation transition, necessary for their catalytic functions. Herein, we describe microsecond molecular dynamics simulations of two distantly related glycosyltransferases that are part of the GT-B structural superfamily, HepI and GtfA. Simulations were performed using the open and closed conformations of these unbound proteins, respectively, and we sought to identify the major dynamical modes and communication networks that interconnect the open and closed structures. We provide the first reported evidence within the scope of our simulation parameters that the interconversion between open and closed conformations is a hierarchical multistep process which can be a conserved feature of enzymes of the same structural superfamily. Each of these motions involves of a collection of smaller molecular reorientations distributed across both domains, highlighting the complexities of protein dynamic involved in the interconversion process. Additionally, dynamic cross-correlation analysis was employed to explore the potential effect of distal residues on the catalytic efficiency of HepI. Multiple distal nonionizable residues of the C-terminal domain exhibit motions anticorrelated to positively charged residues in the active site in the N-terminal domain involved in substrate binding. Mutations of these residues resulted in a reduction in negatively correlated motions and an altered enzymatic efficiency that is dominated by lower Km values with kcat effectively unchanged. The findings suggest that residues with opposing conformational motions involved in the opening and closing of the bidomain HepI protein can allosterically alter the population and conformation of the “closed” state, essential to the formation of the Michaelis complex. The stabilization effects of these mutations likely equally influence the energetics of both the ground state and the transition state of the catalytic reaction, leading to the unaltered kcat. Our study provides new insights into the role of conformational dynamics in glycosyltransferase’s function and new modality to modulate enzymatic efficiency.

Dynamic cross correlation analysis of Thermus thermophilus alkaline phosphatase and determinants of thermostability

BackgroundUnderstanding the determinants of protein thermostability is very important both from the theoretical and applied perspective. One emerging view in thermostable enzymes seems to indicate that a salt bridge/charged residue network plays a fundamental role in their thermostability. MethodsThe structure of alkaline phosphatase (AP) from Thermus thermophilus HB8 was solved by X-ray crystallography at 2.1 Å resolution. The obtained structure was further analyzed by molecular dynamics studies at different temperatures (303 K, 333 K and 363 K) and compared to homologous proteins from the cold-adapted organisms Shewanella sp. and Vibrio strain G15–21. To analyze differences in measures of dynamic variation, several data reduction techniques like principal component analysis (PCA), residue interaction network (RIN) analysis and rotamer analysis were used. Using hierarchical clustering, the obtained results were combined to determine residues showing high degree dynamical variations due to temperature jumps. Furthermore, dynamic cross correlation (DCC) analysis was carried out to characterize networks of charged residues. ResultsTop clustered residues showed a higher propensity for thermostabilizing mutations, indicating evolutionary pressure acting on thermophilic organisms. The description of rotamer distributions by Gini coefficients and Kullback–Leibler (KL) divergence both revealed significant correlations with temperature. DCC analysis revealed a significant trend to de-correlation of the movement of charged residues at higher temperatures. SignificanceThe de-correlation of charged residues detected in Thermus thermophilus AP, highlights the importance of dynamic electrostatic network interactions for the thermostability of this enzyme.

Investigating the mechanism of recognition and structural dynamics of nucleoprotein-RNA complex from Peste des petits ruminants virus via Gaussian accelerated molecular dynamics simulations

The nucleocapsid (N) protein from Peste des petits ruminants virus enwraps the nascent genomic RNA primarily to protect it from cellular enzymatic degradation. Here, we have combined molecular modeling and 1 μs long Gaussian accelerated molecular dynamics simulations to study the structural and conformational properties of the apo N-protein and three protein-RNA complexes, namely wild type (WT), MT1 (I46T/E297D/G342E), and MT2 (I46T/E297D/G342E/W333G). The root-mean-squared deviation (RMSD) analysis reveals that although MT2 deviates most significantly from the initial structure, the RNA binding region is more stable compared to WT or MT1. Further, the flexible nature of the N protein is revealed from the principal component analysis. Our study shows that the solvent accessible surface area of the binding region increases drastically for both mutant complexes compared to WT. The dynamic cross-correlation analysis suggests that the overall anticorrelated motion weakens after the RNA binding. MT2 displays a relatively larger positive correlation than WT or MT1. The binding free energy is estimated for all three complexes via the molecular mechanics/generalized Born surface (MM/GBSA) scheme, and it decreases in the order WT > MT1 > MT2. In all cases, the protein-RNA binding is mainly driven by the van der Waals interactions since the intermolecular electrostatic interaction is overcompensated by desolvation energy. All hotspot residues are identified from the per-residue decomposition of the total binding free energy. In MT2, RNA contributes most favorably compared to WT or MT1. We believe our study will help in understanding the mechanism of RNA recognition by N proteins. Communicated by Ramaswamy H. Sarma

New Insights into the Binding Mechanism of Co-regulator BUD31 to AR AF2 Site: Structural Determination and Analysis of the Mutation Effect.

Introduction Androgen Receptor (AR) plays a pivotal role in the development of male sex and contributes to prostate cancer growth. Different from other nuclear receptors that bind to the co-regulator LxxLL motif in coregulator peptide interaction, the AR Ligand Binding Domain (LBD) prefers to bind to the FxxLF motif. BUD31, a novel co-regulator with FxxLF motif, has been demonstrated to suppress wild-type and mutated AR-mediated prostate cancer growth.Methods To find out the interaction mechanisms of BUD31 with WT/T877A/W741L AR complex, molecular dynamics simulations were employed to study the complex BUD31 and WT/mutant ARs. The molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) results demonstrated that T877A and W741L point mutations can reduce the binding affinity between BUD31 and AR. The RMSF and dynamic cross-correlation analysis indicated that amino acid point mutations can affect the motions of loop residues in the AR structure.Results These results indicated that AR co-regulator binding site AF2 can serve as a target for drug discovery to solve the resistance problem.

A beginner's guide to molecular dynamics simulations and the identification of cross-correlation networks for enzyme engineering.

The functional properties of proteins are decided not only by their relatively rigid overall structures, but even more importantly, by their dynamic properties. In a protein, some regions of structure exhibit highly correlated or anti-correlated motions with others, some are highly dynamic but uncorrelated, while other regions are relatively static. The residues with correlated or anti-correlated motions can form a so-called dynamic cross-correlation network, through which information can be transmitted. Such networks have been shown to be critical to allosteric transitions, and ligand binding, and have also been shown to be able to mediate epistatic interactions between mutations. As a result, they are likely to play a significant role in the development of new enzyme engineering strategies. In this chapter, protocols are provided for the assessment of dynamic cross-correlation networks, and for their application in protein engineering. Transketolase from E. coli is used as a model and the software GROMACS is applied for carrying out MD simulations to generate trajectories containing structural ensembles. The trajectory is then used for a dynamic cross correlation analysis using the R package, Bio3D. A matrix of all atom-wise cross-correlation coefficients is finally obtained, which can be displayed in a graphical representation termed a dynamical cross-correlation matrix.

Investigating Phosphorylation-Induced Conformational Changes in WNK1 Kinase by Molecular Dynamics Simulations.

The With-No-Lysine (WNK) kinase is considered to be a master regulator for various cation-chloride cotransporters involved in maintaining cell-volume and ion homeostasis. Here, we have investigated the phosphorylation-induced structural dynamics of the WNK1 kinase bound to an inhibitor via atomistic molecular dynamics simulations. Results from our simulations show that the phosphorylation at Ser382 could stabilize the otherwise flexible activation loop (A-loop). The intrahelix salt-bridge formed between Arg264 and Glu268 in the unphosphorylated system is disengaged after the phosphorylation, and Glu268 reorients itself and forms a stable salt-bridge with Arg348. The dynamic cross-correlation analysis shows that phosphorylation diminishes anticorrelated motions and increases correlated motions between different domains. Structural network analysis reveals that the phosphorylation causes structural rearrangements and shortens the communication path between the αC-helix and catalytic loop, making the binding pocket more suitable for accommodating the ligand. Overall, we have characterized the structural changes in the WNK kinase because of phosphorylation in the A-loop, which might help in designing rational drugs.

Conformation and Domain Movement Analysis of Human Matrix Metalloproteinase-2: Role of Associated Zn2+ and Ca2+ Ions.

Matrix metaloproteinase-2 (MMP-2) is an extracellular Zn2+ protease specific to type I and IV collagens. Its expression is associated with several inflammatory, degenerative, and malignant diseases. Conformational properties, domain movements, and interactions between MMP-2 and its associated metal ions were characterized using a 1.0 µs molecular dynamics simulation. Dihedral principle component analysis revealed ten families of conformations with the greatest degree of variability occurring in the link region connecting the catalytic and hemopexin domains. Dynamic cross-correlation analysis indicated domain movements corresponding to the opening and closing of the hemopexin domain in relation to the fibronectin and catalytic domains facilitated by the link region. Interaction energies were calculated using the molecular mechanics Poisson Boltzman surface area-interaction entropy (MMPBSA-IE) analysis method and revealed strong binding energies for the catalytic Zn2+ ion 1, Ca2+ ion 1, and Ca2+ ion 3 with significant conformational stability at the binding sites of Zn2+ ion 1 and Ca2+ ion 1. Ca2+ ion 2 diffuses freely away from its crystallographically defined binding site. Zn2+ ion 2 plays a minor role in conformational stability of the catalytic domain while Ca2+ ion 3 is strongly attracted to the highly electronegative sidechains of the Asp residues around the central β-sheet core of the hemopexin domain; however, the interacting residue sidechain carboxyl groups are outside of Ca2+ ion 3′s coordination sphere.