Search

Coronavirus disease 2019 (COVID-19) is an ongoing pandemic. The virus that causes the disease, severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), predominantly infects the respiratory tract, which may lead to pneumonia and death in severe cases. Many marine compounds have been found to have immense medicinal value and have gained approval from the Food and Drug Administration (FDA), and some are being tested in clinical trials. In the current investigation, we redirected a number of marine compounds toward SARS-CoV-2 by targeting the main protease (Mpro, PDB ID: 6Y2F), subjecting them to several advanced computational techniques using co-crystallised ligand as the reference compound. The results of the binding affinity studies showed that two compounds, eribulin mesylate (eri) and soblidotin (sob), displayed higher docking scores than did the reference compound. When these compounds were assessed using molecular dynamics simulation, it was evident that they demonstrated stable binding at the binding pocket of the target protein. The systems demonstrated stable root mean square deviation and radius of gyration values, while occupying the binding pocket during the simulation run. Furthermore, the essential dynamics and free energy landscape exploration revealed that the protein had navigated through a minimal energy basin and demonstrated favourable conformation while binding to the proposed inhibitors. Collectively, our findings suggest that two marine compounds, namely eri and sob, show potential as SARS-CoV-2 main protease inhibitors., Graphical abstract Image 1

Human triokinase/flavin mononucleotide (FMN) cyclase (hTKFC) catalyzes the adenosine triphosphate (ATP)-dependent phosphorylation of D-glyceraldehyde and dihydroxyacetone (DHA), and the cyclizing splitting of flavin adenine dinucleotide (FAD). hTKFC structural models are dimers of identical subunits, each with two domains, K and L, with an L2-K1-K2-L1 arrangement. Two active sites lie between L2-K1 and K2-L1, where triose binds K and ATP binds L, although the resulting ATP-to-triose distance is too large (&asymp, 14 &Aring, ) for phosphoryl transfer. A 75-ns trajectory of molecular dynamics shows considerable, but transient, ATP-to-DHA approximations in the L2-K1 site (4.83 &Aring, or 4.16 &Aring, ). To confirm the trend towards site closure, and its relationship to kinase activity, apo-hTKFC, hTKFC:2DHA:2ATP and hTKFC:2FAD models were submitted to normal mode analysis. The trajectory of hTKFC:2DHA:2ATP was extended up to 160 ns, and 120-ns trajectories of apo-hTKFC and hTKFC:2FAD were simulated. The three systems were comparatively analyzed for equal lengths (120 ns) following the principles of essential dynamics, and by estimating site closure by distance measurements. The full trajectory of hTKFC:2DHA:2ATP was searched for in-line orientations and short distances of DHA hydroxymethyl oxygens to ATP &gamma, phosphorus. Full site closure was reached only in hTKFC:2DHA:2ATP, where conformations compatible with an associative phosphoryl transfer occurred in L2-K1 for significant trajectory time fractions.

The Epstein–Barr virus (EBV) is associated with a variety of malignancies and chronic active EBV infection is a severe systemic disease associated with high rates of mortality and morbidity. In this paper, the dynamics of the interaction of EBV-expressed latent membrane protein 1 (LMP1) with cellular signaling intermediate tumor necrosis factor receptor (TNFR)-associated factor 3 (TRAF3) is simulated using standard classical molecular dynamics (MD) protocols. For comparison, the dynamics of the interaction of TRAF3 with CD40, a TNFR mimicked by LMP1 to effect EBV infection is also calculated under similar conditions. Essential dynamics (ED) analysis is carried out to identify important degrees of vibrational freedom that relate to protein function and virus infection. Both the MD simulation and ED analysis reveal novel interactions that help explain the structural decoy action of LMP1 over CD40. These interactions involve the consensus sequence PXQXTXX shared by CD40 and LMP1. In LMP1, we have found novel interaction of Asp 209 with TRAF3 and the interaction is crucial although the adjacent Asp 210 was suggested to be essential by the X-ray analysis. In CD40, it is found that the hairpin formation is not indispensable for the interaction with TRAF3.

The main purpose of this work is to analyse, by means of molecular dynamics (MD) simulations both structural and mechanical-dynamical differences between Hepcidin-20 and Hepcidin-25 in both oxidized and reduced states in aqueous solution. Results indicate that the presence of disulfide bonds is essential, in both peptides, for maintaining their β-hairpin motif. As a matter of fact, the lack of this intra-peptide covalent interactions produces an almost immediate deviation from the oxidized, plausibly active, structure in both the systems. Interestingly, reduced Hepcidin-25 turns out to be characterized by a highly fluctuating structure which is found to rapidly span a large number of configurations at equilibrium. On the other hand, loss of disulfide bonds in the shorter peptide, results in a more compact and relatively rigid double-turn structure. Comparison of mechanical–dynamical properties and sidechains–sidechains interactions in oxidized Hepcidin-20 and Hepcidin-25 strongly suggest also the key role of N-terminus in the aggregation tendency of Hepcidin-25. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 917–926, 2010.

In this work we report the study of a peptide, the Contryphan Vn produced by Conus ventricosus, a vermivorous cone snail living in the temperate Mediterranean sea. This cyclic peptide of nine residues is a ring closed by a Cys-Cys (Cys: cysteine) disulfide bond containing two proline (Pro) residues and two tryptophans (Trp), one of them being a D-Trp. We present a statistical mechanical characterization of the peptide, simulated in water for about 200 ns with classical molecular dynamics (MD). In recent years there has been a growing interest in the study of the mechanics and dynamics of biological molecules, and in particular for proteins and peptides, about the relationship between collective motions and the active conformations which exert the biological function. To this aim we used the essential dynamics analysis on the MD trajectory and extracted, from the total fluctuations of the molecule, the dominant dynamical modes responsible of the principal conformational transitions. The Contryphan Vn small size allowed us to investigate in details the all-atoms dynamics and the corresponding thermodynamics in conformational space defined by the most significant intramolecular motions. © 2004 Wiley Periodicals, Inc. Biopoly- mers 74: 448 - 456, 2004

Essential Dynamics (ED) is a powerful tool for analyzing molecular dynamics (MD) simulations and it is widely adopted for conformational analysis of large molecular systems such as, for example, proteins and nucleic acids. In this study, we extend the use of ED to the study of clusters of arbitrary size constituted by weakly interacting particles, for example, atomic clusters and supramolecular systems. The key feature of the method we present is the identification of the relevant atomic-molecular clusters to be analyzed by ED for extracting the information of interest. The application of this computational approach allows a straightforward and unbiased conformational study of the local microstructures in liquids, as emerged from semiclassical MD simulations. The good performance of the method is demonstrated by calculating typical observables of liquid water, that is, NMR, NEXAFS O1s, and IR spectra, known to be rather sensitive both to the presence and to the conformational features of hydrogen-bonded clusters. � 2014 Wiley Periodicals, Inc.

The structural and dynamical behavior of the 41-56 β-hairpin from the protein G B1 domain (GB1) has been studied at different temperatures using molecular dynamics (MD) simulations in an aqueous environment. The purpose of these simulations is to establish the stability of this hairpin in view of its possible role as a nucleation site for protein folding. The conformation of the peptide in the crystallographic structure of the protein GB1 (native conformation) was lost in all simulations. The new equilibrium conformations are stable for several nanoseconds at 300 K (>10 ns), 350 K (>6.5 ns), and even at 450 K (up to 2.5 ns). The new structures have very similar hairpin-like conformations with properties in agreement with available experimental nuclear Overhauser effect (NOE) data. The stability of the structure in the hydrophobic core region during the simulations is consistent with the experimental data and provides further evidence for the role played by hydrophobic interactions in hairpin structures. Essential dynamics analysis shows that the dynamics of the peptide at different temperatures spans basically the same essential subspace. The main equilibrium motions in this subspace involve large fluctuations of the residues in the turn and ends regions. Of the six interchain hydrogen bonds, the inner four remain stable during the simulations. The space spanned by the first two eigenvectors, as sampled at 450 K, includes almost all of the 47 different hairpin structures found in the database. Finally, analysis of the hydration of the 300 K average conformations shows that the hydration sites observed in the native conformation are still well hydrated in the equilibrium MD ensemble.

The implementation of cutinase from Fusarium solani pisi as a fat-stain removing ingredient in laundry washing is hampered by its unfolding in the presence of anionic surfactants. In this work we present molecular dynamics (MD) computer simulations on cutinase and analysis procedures to distinguish the movements related to its functional behavior (e.g., substrate binding) from those related to the unfolding of the enzyme. Two kinds of MD-simulations were performed: a simulation mimicking the thermal motion at room temperature, and several simulations in which unfolding is induced either by high temperature or by using a modified water-protein interaction potential. Essential dynamics analyses (A. Amadei et al., Proteins 17:412–425, 1993) on the simulations identify distinct regions in the molecular structure of cutinase in which the motions occur for function and initial unfolding. The unfolding in various simulations starts in a similar way, suggesting that mutations in the regions involved might stabilize the enzyme without affecting its functionality. Proteins 33:253–264, 1998. © 1998 Wiley-Liss, Inc.

Model-free methods are introduced to determine quantities pertaining to protein domain motions from normal mode analyses and molecular dynamics simulations, For the normal mode analysis, the methods are based on the assumption that in low frequency modes, domain motions can be well approximated by modes of motion external to the domains, To analyze the molecular dynamics trajectory, a principal component analysis tailored specifically to analyze interdomain motions is applied, A method based on the curl of the atomic displacements is described, which yields a sharp discrimination of domains, and which defines a unique interdomain screw-axis, Hinge axes are defined and classified as twist or closure axes depending on their direction, The methods have been tested on lysozyme, A remarkable correspondence was found between the first normal mode axis and the first principal mode axis, with both axes passing within 3 Angstrom of the alpha-carbon atoms of residues 2, 39, and 56 of human lysozyme, and near the interdomain helix, The axes of the first modes are overwhelmingly closure axes, A lesser degree of correspondence is found for the second modes, but in both cases they are more twist axes than closure axes, Both analyses reveal that the interdomain connections allow only these two degrees of freedom, one more than provided by a pure mechanical hinge. (C) 1997 Wiley-Liss, Inc.

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 21KTS23 motif, and the C-terminal tail are known to be responsible for the selective binding to the a1 beta 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. (c) 2012 Wiley Periodicals, Inc.

The mechanism of nucleotide-regulated assembly and disassembly of the prokaryotic cell division protein FtsZ is not yet clearly understood. In this work, we attempt to characterize the functional motions in monomeric FtsZ through molecular dynamics simulations and essential dynamics (ED) analyses and correlate those motions to FtsZ assembly and disassembly. Results suggest that the nucleotide binding subdomain of FtsZ can switch between multitudes of curved conformations in all nucleotide states, but it prefers to be in an assembly competent less curved conformation in the GTP-bound state. Further, the GDP to GTP exchange invokes a subtle conformational change in the nucleotide binding pocket that tends to align the top portion of core helix H7 along the longitudinal axis of the protein. ED analyses suggest that the longitudinal movements of H7 and the adjacent H6-H7 region modulate the motions of C-domain elements coherently. These longitudinal movements of functionally relevant H7, H6-H7, T3, T7, and H10 regions are likely to facilitate the assembly of GTP-FtsZ into straight filament. On the other hand, the observed radial or random movements of FtsZ residues in the GDP state might not allow the monomers to assemble as efficiently as GTP-bound monomers and could produce curved filaments. Our results correlate very well with recent mutagenesis data that inferred FtsZ conformational flexibility and the involvement of the H6-H7 region in assembly. � 2013 American Chemical Society.

The Gal4p mediated transcriptional activation of GAL genes requires the interaction between Gal3p bound with ATP and galactose and Gal80p. Though numerous studies suggest that galactose and ATP activate Gal3p/Gal1p interaction with Gal80p, neither the mechanism of activation nor the interacting surface that binds to Gal80p is well understood. In this study we investigated the dynamics of Gal3p and Gal1p in the presence and absence of ligands ATP and galactose to understand the role played by dynamics in the function of these proteins through molecular dynamics simulation and protein–protein docking studies. We performed simulations totaling to 510 ns on both Gal1p and Gal3p proteins in the presence and absence of ligands ATP and galactose. We find that, while binding of ligands ATP and galactose to Gal3p/Gal1p do not affect the global conformation of proteins, some local conformational changes around upper-lip helix including insertion domain are observed. We observed that only in the presence of ATP and galactose, Gal3p displays opening and closing motion between the two domains. And because of this motion, a binding interface, which is largely hydrophobic, opens up on the surface of Gal3p and this surface can bind to Gal80p. From our simulation studies we infer probable docking sites for Gal80p on Gal3p/Gal1p, which were further ascertained by the docking of Gal80p on to ligand bound Gal1p and Gal3p proteins, and the residues at the interface between Gal3p and Gal80p are identified. Our results correlate quite well with the existing body of literature on functional and dynamical aspects of Gal1p and Gal3p proteins.

NMR structures are typically deposited in databases such as the PDB in the form of an ensemble of structures. Generally, each of the models in such an ensemble satisfies the experimental data and is equally valid. No unique solution can be calculated because the experimental NMR data is insufficient, in part because it reflects the conformational variability and dynamical behavior of the molecule in solution. Even for relatively rigid molecules, the limited number of structures that are typically deposited cannot completely encompass the structural diversity allowed by the observed NMR data, but they can be chosen to try and maximize its representation. We describe here the adaptation and application of techniques more commonly used to examine large ensembles from molecular dynamics simulations, to the analysis of NMR ensembles. The approach, which is based on principal component analysis, we call COCO ("Complementary Coordinates"). The COCO approach analyses the distribution of an NMR ensemble in conformational space, and generates a new ensemble that fills "gaps" in the distribution. The method is very rapid, and analysis of a 25-member ensemble and generation of a new 25 member ensemble typically takes 1-2 min on a conventional workstation. Applied to the 545 structures in the RECOORD database, we find that COCO generates new ensembles that are as structurally diverse-both from each other and from the original ensemble-as are the structures within the original ensemble. The COCO approach does not explicitly take into account the NMR restraint data, yet in tests on selected structures from the RECOORD database, the COCO ensembles are frequently good matches to this data, and certainly are structures that can be rapidly refined against the restraints to yield high-quality, novel solutions. COCO should therefore be a useful aid in NMR structure refinement and in other situations where a richer representation of conformational variability is desired-for example in docking studies. COCO is freely accessible via the website www.ccpb.ac.uk/COCO.

The correlation between protein motions and function is a central problem in protein science. Several studies have demonstrated that ligand binding and protein dynamics are strongly correlated in intracellular lipid binding proteins (iLBPs), in which the high degree of flexibility, principally occurring at the level of helix-II, CD, and EF loops (the so-called portal area), is significantly reduced upon ligand binding. We have recently investigated by NMR the dynamic properties of a member of the iLBP family, chicken liver bile acid binding protein (cL-BABP), in its apo and holo form, as a complex with two bile salts molecules. Binding was found to be regulated by a dynamic process and a conformational rearrangement was associated with this event. We report here the results of molecular dynamics (MD) simulations performed on apo and holo cL-BABP with the aim of further characterizing the protein regions involved in motion propagation and of evaluating the main molecular interactions stabilizing bound ligands. Upon binding, the root mean square fluctuation values substantially decrease for CD and EF loops while increase for the helix-loop-helix region, thus indicating that the portal area is the region mostly affected by complex formation. These results nicely correlate with backbone dynamics data derived from NMR experiments. Essential dynamics analysis of the MD trajectories indicates that the major concerted motions involve the three contiguous structural elements of the portal area, which however are dynamically coupled in different ways whether in the presence or in the absence of the ligands. Motions of the EF loop and of the helical region are part of the essential space of both apo and holo-BABP and sample a much wider conformational space in the apo form. Together with NMR results, these data support the view that, in the apo protein, the flexible EF loop visits many conformational states including those typical of the holo state and that the ligand acts stabilizing one of these pre-existing conformations. The present results, in agreement with data reported for other iLBPs, sharpen our knowledge on the binding mechanism for this protein family.

The active site of thermolysin-like proteases (TLPs) is located at the bottom of a cleft between the N- and C-terminal domains. Crystallographic studies have shown that the active-site cleft is more closed in Ligand-binding TLPs than in Ligand-free TLPs. Accordingly, it has been proposed that TLPs undergo a hinge-bending motion during catalysis resulting in "closure" and "opening" of the active-site cleft. Two hinge regions have been proposed. One is located around a conserved glycine 78; the second involves residues 135 and 136. The importance of conserved glycine residues in these hinge regions was studied experimentally by analyzing the effects of Gly --> Ala mutations on catalytic activity. Eight such mutations were made in the TLP of Bacillus stearothermophilus (TLP-ste) and their effects on activity toward casein and various peptide substrates were determined. Only the Gly78Ala, Gly136Ala, and Gly135Ala + Gly136Ala mutants decreased catalytic activity significantly. These mutants displayed a reduction in k(cat)/K-m for 3-(2-furylacryloyl)-L-glycyl-L-leucine amide of 73%, 62%, and 96%, respectively. Comparisons of effects on k(cat)/K-m for various substrates with effects on the K-i for phosphoramidon suggested that the mutation at position 78 primarily had an effect on substrate binding, whereas the mutations at positions 135 and 136 primarily influence k(cat). The apparent importance of conserved glycine residues in proposed hinge-bending regions for TLP activity supports the idea that hinge-bending is an essential part of catalysis.

A method is presented that generates random protein structures that fulfil a set of upper and lower interatomic distance limits. These limits depend on distances mea- sured in experimental structures and the strength of the interatomic interaction. Struc- tural differences between generated struc- tures are similar to those obtained from experi- ment and from MD simulation. Although detailed aspects of dynamical mechanisms are not covered and the extent of variations are only estimated in a relative sense, applications to an IgG-binding domain, an SH3 binding domain, HPr, calmodulin, and lysozyme are presented which illustrate the use of the method as a fast and simple way to predict structural variability in proteins. The method may be used to support the design of mutants, when structural fluctuations for a large num- ber of mutants are to be screened. The results suggest that motional freedom in proteins is ruled largely by a set of simple geometric constraints. Proteins 29:240-251, 1997. r 1997 Wiley-Liss, Inc.

The crystal structure of calcium-calmodulin (CaM) reveals a protein with a typical dumbbell structure. Various spectroscopic studies have suggested that the central linker region of CaM, which is alpha-helical in the crystal structure, is flexible in solution. In particular, NMR studies have indicated the presence of a flexible backbone between residues Lys 77 and Asp 80. This flexibility is related directly to the function of the protein because it enables the N- and C-terminal domains of the protein to move toward each other and bind to the CaM-binding domain of a target protein. We have investigated the flexibility of the CaM central helix by a variety of computational techniques: molecular dynamics (MD) simulations, normal mode analysis (NMA), and essential dynamics (ED) analysis. Our MD results reproduce the experimentally determined location of the bend in a simulation of only the CaM central helix, indicating that the bending point is an intrinsic property of the alpha-helix, for which the remainder of the protein is not important. Interestingly, the modes found by the ED analysis of the MD trajectory are very similar to the lowest frequency modes from the NM analysis and to modes found by an ED analysis of different structures in a set of NMR structures. Electrostatic interactions involving residues Arg 74 and Asp 80 seem to be important for these bending motions and unfolding, which is in line with pH-dependent NMR and CD studies.

We address the question how well proteins can be modelled on the basis of NMR data, when these data are incorporated into the protein model using distance restraints in a molecular dynamics simulation. We found, using HPr as a model protein, that distance restraining freezes the essential motion of proteins, as defined by Amadei et al. [Amadei, A., Linssen, A.B.M. and Berendsen, H.J.C. (1993) Protein Struct. Funct. Genet., 17, 412-425]. We discuss how modelling protocols can be improved in order to solve this problem.

The effect of solvent hydrophobicity on activation of Candida rugosa lipase (CRL) was investigated by performing molecular dynamics simulations for four nano seconds (ns). The closed/inactive conformer of CRL (PDB code 1TRH) was solvated in three alkane–aqueous environments. The alkanes aggregated in a predominantly aqueous environment and by 1ns a stable spherical alkane–aqueous interface had formed. This led to the interfacial activation of CRL. On analyzing the simulated conformers with the closed conformer of CRL, the flap was found to have opened from a closed state by 7.7Å, 10.2Å, 13.1Å at hexane–aqueous, octane–aqueous, and decane–aqueous interfaces. Further, essential dynamics analysis revealed that major anharmonic fluctuations were confined to residues 64–81, the flap of CRL.

Background To understand the effect of the long intracellular loop 3 (ICL3) on the intrinsic dynamics of human β2-adrenergic receptor, molecular dynamics (MD) simulations were performed on two different models, both of which were based on the inactive crystal structure in complex with carazolol (after removal of carazolol and T4-lysozyme). In the so-called loop model, the ICL3 region that is missing in available crystal structures was modeled as an unstructured loop of 32-residues length, whereas in the clipped model, the two open ends were covalently bonded to each other. The latter model without ICL3 was taken as a reference, which has also been commonly used in recent computational studies. Each model was embedded into POPC bilayer membrane with explicit water and subjected to a 1 μs molecular dynamics (MD) simulation at 310 K. Results After around 600 ns, the loop model started a transition to a “very inactive” conformation, which is characterized by a further movement of the intracellular half of transmembrane helix 6 (TM6) towards the receptor core, and a close packing of ICL3 underneath the membrane completely blocking the G-protein’s binding site. Concurrently, the binding site at the extracellular part of the receptor expanded slightly with the Ser207-Asp113 distance increasing to 18 Å from 11 Å, which was further elaborated by docking studies. Conclusions The essential dynamics analysis indicated a strong coupling between the extracellular and intracellular parts of the intact receptor, implicating a functional relevance for allosteric regulation. In contrast, no such transition to the “very inactive” state, nor any structural correlation, was observed in the clipped model without ICL3. Furthermore, elastic network analysis using different conformers for the loop model indicated a consistent picture on the specific ICL3 conformational change being driven by global modes.