Search

Molecular modeling techniques have become indispensable in many fields of molecular sciences in which the details related to mechanisms and reactivity need to be studied at an atomistic level. This review article provides a collection of computational modeling works on a topic of enormous interest and urgent relevance: the properties of metalloenzymes involved in the degradation and valorization of natural biopolymers and synthetic plastics on the basis of both circular biofuel production and bioremediation strategies. In particular, we will focus on lytic polysaccharide monooxygenase, laccases, and various heme peroxidases involved in the processing of polysaccharides, lignins, rubbers, and some synthetic polymers. Special attention will be dedicated to the interaction between these enzymes and their substrate studied at different levels of theory, starting from classical molecular docking and molecular dynamics techniques up to techniques based on quantum chemistry. [ABSTRACT FROM AUTHOR]

Catalytic H2 oxidation has been dissected by means of DFT into the key steps common to the Fe2 unit of both the [FeFe]‐hydrogenase cofactor and selected biomimics. The aim was to elucidate the molecular details underlying the very different performances of the two systems. We found that the better enzyme performance is based on a single iron atom that is maintained electron‐poor, favoring H2 binding, although embedded within a highly electron‐rich cofactor, ensuring a facile oxidation of the Fe2–H2 adduct. This is due to 1) CN− coordinating to both iron atoms, due to their amphipathic Lewis acid/base properties, and 2) the 4Fe4S subunit further withdrawing electrons from the Fe2 core. Preserving a moderate electron deficiency at a single iron also helps the cofactor preserve hydride affinity, which favors H2 cleavage. Such valuable characteristics allow the biocatalyst to turnover close to equilibrium conditions. All previous biomimicry has shown, in contrast, the impossibility to properly balance the two apparently contrasting aforementioned requisites, although evident progress has been made by the H2‐ase community. Disclosure of the differences identified could inspire the design of novel biomimics, for instance, reconsidering the use of CN− in the catalyst architecture. Indeed, in the presence of bases normally employed in oxidative catalysis, undesired stable protonation at coordinated CN−, which affects the opposite process (proton reduction), could be overcome. Nature versus biomimicry: Dissecting H2 oxidation by the [FeFe]‐H2ase core in comparison with a series of Fe2S2 models has revealed that CN− ligands and the 4Fe4S subunit are the key to the higher catalytic performance of the enzyme compared with the synthetic compounds. Indeed, cyanides and cubane provide a single electron‐poor iron (which assists H2 binding) embedded within an electron‐rich environment, which in turn ensures a facile subsequent oxidation. [ABSTRACT FROM AUTHOR]

Among the different classes of enzymes involved in the ubiquitin pathway, E2 ubiquitin-conjugating enzymes occupy a central role in the ubiquitination cascade. Cdc34-like E2 enzymes are characterized by a 12–14 residue insertion in the proximity of the catalytic site, known as the acidic loop. Cdc34 ubiquitin-charging activity is regulated by CK2-dependent phosphorylation and the regulatory mechanism involves the acidic loop. Indeed, the phosphorylation stabilizes the loop in an open conformation that is competent for ubiquitin charging. Cdc34 is associated with a variety of diseases, such as hepatocellular carcinomas and prostatic adenocarcinomas. In light of its role, the discovery of potential inhibitory compounds would provide the mean to effectively modulate its activity. Here, we carried out a computational study based on molecular dynamics, virtual screening and docking to identify potential inhibitory compounds of Cdc34, modulating the acidic loop conformation. The molecules identified in this study have been designed to act as molecular hinges that can bind the acidic loop in its closed conformation, thus inhibiting the Cdc34-mediated ubiquitination cascade at the ubiquitin-charging step. In particular, we proposed a pharmacophore model featuring two amino groups in the central part of the model and two lateral aromatic chains, which respectively establish electrostatic interactions with the acidic loop (Asp 108 and Glu 109) and a hydrogen bond with Ser 139, which is one of the key residues for Cdc34 activity. [ABSTRACT FROM AUTHOR]

Cyclin-dependent kinase inhibitors (CKIs) are key regulatory proteins of the eukaryotic cell cycle, which modulate cyclin-dependent kinase (Cdk) activity. CKIs perform their inhibitory effect by the formation of ternary complexes with a target kinase and its cognate cyclin. These regulators generally belong to the class of intrinsically disordered proteins (IDPs), which lack a well-defined and organized three-dimensional (3D) structure in their free state, undergoing folding upon binding to specific partners. Unbound IDPs are not merely random-coil structures, but can present intrinsically folded structural units (IFSUs) and collapsed conformations. These structural features can be relevant to protein function in vivo. The yeast CKI Sic1 is a 284-amino acid IDP that binds to Cdk1 in complex with the Clb5,6 cyclins, preventing phosphorylation of G1 substrates and, therefore, entrance to the S phase. Sic1 degradation, triggered by multiple phosphorylation events, promotes cell-cycle progression. Previous experimental studies pointed out a propensity of Sic1 and its isolated domains to populate both extended and compact conformations. The present contribution provides models for compact conformations of the Sic1 kinase-inhibitory domain (KID) by all-atom molecular dynamics (MD) simulations in explicit solvent and in the absence of interactors. The results are integrated by spectroscopic and spectrometric data. Helical IFSUs are identified, along with networks of intramolecular interactions. The results identify a group of putative hub residues and networks of electrostatic interactions, which are likely to be involved in the stabilization of the globular states. [ABSTRACT FROM AUTHOR]

Networks and clusters of intramolecular interactions, as well as their "communication" across the three-dimensional architecture have a prominent role in determining protein stability and function. Special attention has been dedicated to their role in thermal adaptation. In the present contribution, seven previously experimentally characterized mutants of a cold-adapted a-amylase, featuring mesophilic-like behavior, have been investigated by multiple molecular dynamics simulations, essential dynamics and analyses of correlated motions and electrostatic interactions. Our data elucidate the molecular mechanisms underlying the ability of single and multiple mutations to globally modulate dynamic properties of the cold-adapted a-amylase, including both local and complex unpredictable distal effects. Our investigation also shows, in agreement with the experimental data, that the conversion of the cold-adapted enzyme in a warm-adapted variant cannot be completely achieved by the introduction of few mutations, also providing the rationale behind these effects. Moreover, pivotal residues, which are likely to mediate the effects induced by the mutations, have been identified from our analyses, as well as a group of suitable candidates for protein engineering. In fact, a subset of residues here identified (as an isoleucine, or networks of mesophilic-like salt bridges in the proximity of the catalytic site) should be considered, in experimental studies, to get a more efficient modification of the features of the cold-adapted enzyme. [ABSTRACT FROM AUTHOR]

E2 ubiquitin-conjugating enzymes are crucial mediators of protein ubiquitination, which strongly influence the ultimate fate of the target substrates. Recently, it has been shown that the activity of several enzymes of the ubiquitination pathway is finely tuned by phosphorylation, an ubiquitous mechanism for cellular regulation, which modulates protein conformation. In this contribution, we provide the first rationale, at the molecular level, of the regulatory mechanism mediated by casein kinase 2 (CK2) phosphorylation of E2 Cdc34-like enzymes. In particular, we identify two co-evolving signature elements in one of the larger families of E2 enzymes: an acidic insertion in b4a2 loop in the proximity of the catalytic cysteine and two conserved key serine residues within the catalytic domain, which are phosphorylated by CK2. Our investigations, using yeast Cdc34 as a model, through 2.5 μs molecular dynamics simulations and biochemical assays, define these two elements as an important phosphorylation-controlled switch that modulate opening and closing of the catalytic cleft. The mechanism relies on electrostatic repulsions between a conserved serine phosphorylated by CK2 and the acidic residues of the b4a2 loop, promoting E2 ubiquitin charging activity. Our investigation identifies a new and unexpected pivotal role for the acidic loop, providing the first evidence that this loop is crucial not only for downstream events related to ubiquitin chain assembly, but is also mandatory for the modulation of an upstream crucial step of the ubiquitin pathway: the ubiquitin charging in the E2 catalytic cleft. [ABSTRACT FROM AUTHOR]

The design of protein-metal complexes is rapidly advancing, with applications spanning catalysis, sensing, and bioremediation. We report a comprehensive investigation of METPsc1, a Miniaturized Electron Transfer Protein, in complex with cadmium. This study elucidates the impact of metal coordination on protein folding and structural dynamics across temperatures from 100 K to 300 K. Our findings reveal that METPsc1, composed of two similar halves stabilized by intramolecular hydrogen bonds, exhibits a unique "clothespin-like" recoil mechanism. This allows it to adapt to metal ions of varying radii, mirroring the flexibility observed in natural rubredoxins. High-resolution crystallography and molecular dynamics simulations unveil concerted backbone motions and subtle temperature-dependent shifts in side-chain conformations, particularly for residues involved in crystal packing. Notably, Cd S bond lengths increase with temperature, correlating with anisotropic motions of the sulfur atoms involved in second-shell hydrogen bonding. This suggests a dynamic role of protein matrix upon redox cycling. These insights into METPsc1 highlight its potential for catalysis and contribute to the designing of artificial metalloproteins with functional plasticity. The METPsc1 (Miniaturized Electron Transfer Protein single chain) designed protein adapts to different metal ion sizes through its unique backbone motion as shown by temperature-dependent experimental and computational analyses. [Display omitted] • METPsc1, a single-chain Miniaturized Electron Transfer Protein, adapts to metals of varying radii through a "clothespin-like" recoil mechanism. • Temperature-dependent X-ray diffraction reveals subtle conformational changes. • Molecular dynamics simulations reveal a glass transition around 190 K in aqueous solution. • METPsc1 demonstrates adaptability similar to natural rubredoxins. • Findings validate the design approach for creating functional protein-metal complexes. [ABSTRACT FROM AUTHOR]

Abstract: Adaptation to both high and low temperatures requires proteins with special properties. While organisms living at or close to the boiling point of water need to have proteins with increased stability, other properties are required at temperatures close to the freezing point of water. Indeed, it has been shown that enzymes adapted to cold environments are less resistant to heat with a concomitant increased activity as compared to their warm-active counter-parts. Several recent studies have pointed in the direction that electrostatic interactions play a central role in temperature adaptation, and in this study we investigate the role such interactions have in adaptation of elastase from Atlantic salmon and pig. Molecular dynamics (MD) simulations have been used to generate structural ensembles at 283 and 310K of the psychrophilic and mesophilic elastase, and a total of eight 12ns simulations have been carried out. Even though the two homologues have a highly similar three-dimensional structure, the location and number of charged amino acids are very different. Based on the simulated structures we find that very few salt-bridges are stable throughout the simulations, and provide little stabilization/destabilization of the proteins as judged by continuum electrostatic calculations. However, the mesophilic elastase is characterized by a greater number of salt-bridges as well as a putative salt-bridge network close to the catalytic site, indicating a higher rigidity of the components involved in the catalytic cycle. In addition, subtle differences are also found in the electrostatic potentials in the vicinity of the catalytic residues, which may explain the increased catalytic efficiency of the cold-adapted elastase. [Copyright &y& Elsevier]

Protein structure and dynamics are crucial for protein function. Thus, the study of conformational properties can be very informative for characterizing new proteins and to rationalize how residue substitutions at specific protein sites affect its dynamics, activity and thermal stability. Here, we investigate the structure and dynamics of the recently isolated cold-adapted acylaminoacyl peptidase from Sporosarcina psychrophila (SpAAP) by the integration of simulations, circular dichroism, mass spectrometry and other experimental data. Our study notes traits of cold-adaptation, such as lysine-to-arginine substitutions and a lack of disulphide bridges. Cold-adapted enzymes are generally characterized by a higher number of glycine residues with respect to their warm-adapted counterparts. Conversely, the SpAAP glycine content is lower than that in the warm-adapted variants. Nevertheless, glycine residues are strategically located in proximity to the functional sites in SpAAP, such as the active site and the linker between the two domains.. In particular, G457 reduces the steric hindrance around the nucleophile elbow. Our results suggest a local weakening of the intramolecular interactions in the cold-adapted enzyme. This study offers a basis for the experimental mutagenesis of SpAAP and related enzymes. The approaches employed in this study may also provide a more general framework to characterize new protein structures in the absence of X-ray or NMR data. [ABSTRACT FROM AUTHOR]

The DFTB and DFT methods are applied to the study ofdifferentforms of the [Cu(HGGG)(Py)] complex in water, with the aim of identifyingthe most stable isomer. The DFTB calculations were possible thanksto a careful parametrization of the atom–atom repulsive energyterms for Cu–H, Cu–C, Cu–N, and Cu–O.The speciation process is carried out by computing different DFTB-steeredmolecular dynamics (SMD) trajectories, each of which ends in a well-defineddifferent form. The last frame of each trajectory is subjected togeometry optimization at both DFTB and DFT levels, leading to a differentisomer. From the corresponding energy values, a rank of relative stabilityof the isomers can be established. The computational protocol developedhere is of general applicability to other metal–peptide systemsand represents a new powerful tool for the study of speciation ofmetal-containing systems in water solution, particularly useful whenthe full characterization of the compound cannot be carried out onthe basis of experimental results only. [ABSTRACT FROM AUTHOR]

Abstract: Natriuretic peptides (NPs) are a family of structurally related hormone/paracrine factors (ANP, BNP and CNP), which mediate a broad array of physiological effects by interacting with specific guanylyl cyclase receptors (NPR) and have promising therapeutic and clinical applications. NPs are specific for different NPRs and share a common ring structure in which a disulfide bond between two conserved cysteine residues is formed. Residues within the cyclic loop are largely responsible for receptor selectivity. Structural features of free NPs in solution have not been investigated in details even if their characterization would be very useful in order to identify important aspects related to NPs function and receptor selectivity. In light of the above scenario, we carried out a 0.1μs molecular dynamics investigation of NPs with the aim of providing a high-resolution atomistic view of specific of their conformational ensemble in solution. Our results clearly indicate that NP receptor-bound conformations are not stable solution structure and that induced-fit mechanisms are involved in the formation of NP-NPR complexes. Moreover, in agreement with the current view on strictly relationship between protein dynamics and protein function and activity, it turns out that differences in activity and NPR specificity of CNP and ANP/BNP might be correlated to different amino acid composition of the cyclic loop, propensity to form β-sheet structures, flexibility patterns, dynamics properties and free conformations explored during the simulations. [Copyright &y& Elsevier]

Abstract: We have studied the accessibility of the structural calcium ion in the Burkholderia glumae lipase and the consequences of its removal on the protein conformation by different biophysical techniques (circular dichroism, fluorimetry, and mass spectrometry) and by molecular-dynamics simulations. We show that, in the native protein, calcium is not accessible unless specific flexible loops are displaced, for example, by a temperature increase. Such movements concern the whole calcium-binding pocket and particularly the environment of the coordinating aspartate residue 241. As a consequence of metal depletion the protein unfolds irreversibly and undergoes aggregation. The removal of the metal ion causes major structural transitions and leads to an increase in β-structure, in particular in protein regions that are largely unstructured in the native protein and encompass the calcium coordination residues. [Copyright &y& Elsevier]

Abstract: Several molecular dynamics (MD) simulations were used to sample conformations in the neighborhood of the native structure of holo-myoglobin (holo-Mb), collecting trajectories spanning 0.22μs at 300K. Principal component (PCA) and free-energy landscape (FEL) analyses, integrated by cluster analysis, which was performed considering the position and structures of the individual helices of the globin fold, were carried out. The coherence between the different structural clusters and the basins of the FEL, together with the convergence of parameters derived by PCA indicates that an accurate description of the Mb conformational space around the native state was achieved by multiple MD trajectories spanning at least 0.14μs. The integration of FEL, PCA, and structural clustering was shown to be a very useful approach to gain an overall view of the conformational landscape accessible to a protein and to identify representative protein substates. This method could be also used to investigate the conformational and dynamical properties of Mb apo-, mutant, or delete versions, in which greater conformational variability is expected and, therefore identification of representative substates from the simulations is relevant to disclose structure–function relationship. [Copyright &y& Elsevier]

Abstract: Molecular dynamics simulations of representative mesophilic and psycrophilic elastases have been carried out at different temperatures to explore the molecular basis of cold adaptation inside a specific enzymatic family. The molecular dynamics trajectories have been compared and analyzed in terms of secondary structure, molecular flexibility, intramolecular and protein–solvent interactions, unravelling molecular features relevant to rationalize the efficient catalytic activity of psychrophilic elastases at low temperature. The comparative molecular dynamics investigation reveals that modulation of the number of protein–solvent interactions is not the evolutionary strategy followed by the psycrophilic elastase to enhance catalytic activity at low temperature. In addition, flexibility and solvent accessibility of the residues forming the catalytic triad and the specificity pocket are comparable in the cold- and warm-adapted enzymes. Instead, loop regions with different amino acid composition in the two enzymes, and clustered around the active site or the specificity pocket, are characterized by enhanced flexibility in the cold-adapted enzyme. Remarkably, the psycrophilic elastase is characterized by reduced flexibility, when compared to the mesophilic counterpart, in some scattered regions distant from the functional sites, in agreement with hypothesis suggesting that local rigidity in regions far from functional sites can be beneficial for the catalytic activity of psychrophilic enzymes. [Copyright &y& Elsevier]