A Comparison of Bonded and Nonbonded Zinc(II) Force Fields with NMR Data.

Classical molecular dynamics (MD) simulations are widely used to inspect the behavior of zinc(II)-proteins at the atomic level, hence the need to properly model the zinc(II) ion and the interaction with its ligands. Different approaches have been developed to represent zinc(II) sites, with the bonded and nonbonded models being the most used. In the present work, we tested the well-known zinc AMBER force field (ZAFF) and a recently developed nonbonded force field (NBFF) to assess how accurately they reproduce the dynamic behavior of zinc(II)-proteins. For this, we selected as benchmark six zinc-fingers. This superfamily is extremely heterogenous in terms of architecture, binding mode, function, and reactivity. From repeated MD simulations, we computed the order parameter (S2) of all backbone N-H bond vectors in each system. These data were superimposed to heteronuclear Overhauser effect measurements taken by NMR spectroscopy. This provides a quantitative estimate of the accuracy of the FFs in reproducing protein dynamics, leveraging the information about the protein backbone mobility contained in the NMR data. The correlation between the MD-computed S2 and the experimental data indicated that both tested FFs reproduce well the dynamic behavior of zinc(II)-proteins, with comparable accuracy. Thus, along with ZAFF, NBFF represents a useful tool to simulate metalloproteins with the advantage of being extensible to diverse systems such as those bearing dinuclear metal sites.

Learning to Identify Physiological and Adventitious Metal-Binding Sites in the Three-Dimensional Structures of Proteins by Following the Hints of a Deep Neural Network.

Thirty-eight percent of protein structures in the Protein Data Bank contain at least one metal ion. However, not all these metal sites are biologically relevant. Cations present as impurities during sample preparation or in the crystallization buffer can cause the formation of protein-metal complexes that do not exist in vivo. We implemented a deep learning approach to build a classifier able to distinguish between physiological and adventitious zinc-binding sites in the 3D structures of metalloproteins. We trained the classifier using manually annotated sites extracted from the MetalPDB database. Using a 10-fold cross validation procedure, the classifier achieved an accuracy of about 90%. The same neural classifier could predict the physiological relevance of non-heme mononuclear iron sites with an accuracy of nearly 80%, suggesting that the rules learned on zinc sites have general relevance. By quantifying the relative importance of the features describing the input zinc sites from the network perspective and by analyzing the characteristics of the MetalPDB datasets, we inferred some common principles. Physiological sites present a low solvent accessibility of the aminoacids forming coordination bonds with the metal ion (the metal ligands), a relatively large number of residues in the metal environment (≥20), and a distinct pattern of conservation of Cys and His residues in the site. Adventitious sites, on the other hand, tend to have a low number of donor atoms from the polypeptide chain (often one or two). These observations support the evaluation of the physiological relevance of novel metal-binding sites in protein structures.

Insights into the Dynamics of the Human Zinc Transporter ZnT8 by MD Simulations.

ZnT8 is a human zinc(II) transporter expressed at the membrane of secretory granules where it contributes to insulin storage importing zinc ions from the cytosol. In the human population, the two most common ZnT8 variants carry an arginine (R325) or a tryptophan (W325) in position 325. The former variant has the most efficient kinetics in zinc transport and has been correlated to a higher risk of developing insulin resistance. On the contrary, the W325 variant is less active and protects against type-2-diabetes. Here, we used molecular dynamics (MD) simulations to investigate the main differences between the R325 and W325 variants in the interaction with zinc(II) ions. Our simulations suggested that the position of the metal ion within the transport site was not the same for the two variants, underlying a different rearrangement of the transmembrane (TM) helices in the channel. The W325 variant featured a peculiar zinc environment not detected in the experimental structures. With respect to conformational dynamics, we observed that the R325 variant was significantly more flexible than W325, with the main role played by the transmembrane domain (TMD) and the C-terminal domain (CTD). This dynamics affected the packing of the TM helices and thus the channel accessibility from the cytosol. The dimer interface that keeps the two TM channels in contact became looser in both variants upon zinc binding to the transport site, suggesting that this may be an important step toward the switch from the inward- to the outward-facing state of the protein.

Automated Determination of Nuclear Magnetic Resonance Chemical Shift Perturbations in Ligand Screening Experiments: The PICASSO Web Server.

Nuclear magnetic resonance (NMR) is an effective, commonly used experimental approach to screen small organic molecules against a protein target. A very popular method consists of monitoring the changes of the NMR chemical shifts of the protein nuclei upon addition of the small molecule to the free protein. Multidimensional NMR experiments allow the interacting residues to be mapped along the protein sequence. A significant amount of human effort goes into manually tracking the chemical shift variations, especially when many signals exhibit chemical shift changes and when many ligands are tested. Some computational approaches to automate the procedure are available, but none of them as a web server. Furthermore, some methods require the adoption of a fairly specific experimental setup, such as recording a series of spectra at increasing small molecule:protein ratios. In this work, we developed a tool requesting a minimal amount of experimental data from the user, implemented it as an open-source program, and made it available as a web application. Our tool compares two spectra, one of the free protein and one of the small molecule:protein mixture, based on the corresponding peak lists. The performance of the tool in terms of correct identification of the protein-binding regions has been evaluated on different protein targets, using experimental data from interaction studies already available in the literature. For a total of 16 systems, our tool achieved between 79% and 100% correct assignments, properly identifying the protein regions involved in the interaction.

Structural insights into the molecular function of human [2Fe-2S] BOLA1-GRX5 and [2Fe-2S] BOLA3-GRX5 complexes.

Members of the monothiol glutaredoxin family and members of the BolA-like protein family have recently emerged as specific interacting partners involved in iron-sulfur protein maturation and redox regulation pathways. It is known that human mitochondrial BOLA1 and BOLA3 form [2Fe-2S] cluster-bridged dimeric heterocomplexes with the monothiol glutaredoxin GRX5. The structure and cluster coordination of the two [2Fe-2S] heterocomplexes as well as their molecular function are, however, not defined yet. Experimentally-driven structural models of the two [2Fe-2S] cluster-bridged dimeric heterocomplexes, the relative stability of the two complexes and the redox properties of the [2Fe-2S] cluster bound to these complexes are here presented on the basis of UV/vis, CD, EPR and NMR spectroscopies and computational protein-protein docking. While the BOLA1-GRX5 complex coordinates a reduced, Rieske-type [2Fe-2S]1+ cluster, an oxidized, ferredoxin-like [2Fe-2S]2+ cluster is present in the BOLA3-GRX5 complex. The [2Fe-2S] BOLA1-GRX5 complex is preferentially formed over the [2Fe-2S] BOLA3-GRX5 complex, as a result of a higher cluster binding affinity. All these observed differences provide the first indications discriminating the molecular function of the two [2Fe-2S] heterocomplexes.

Investigation of the Iron(II) Release Mechanism of Human H-Ferritin as a Function of pH.

We investigated the kinetics of the release of iron(II) ions from the internal cavity of human H-ferritin as a function of pH. Extensive molecular dynamics simulations of the entire 24-mer ferritin provided atomic-level information on the release mechanism. Double protonation of His residues at pH 4 facilitates the removal of the iron ligands within the C3 channel through the formation of salt bridges, resulting in a significantly lower release energy barrier than pH 9.

A protocol for the refinement of NMR structures using simultaneously pseudocontact shift restraints from multiple lanthanide ions.

The binding of paramagnetic metal ions to proteins produces a number of different effects on the NMR spectra of the system. In particular, when the magnetic susceptibility of the metal ion is anisotropic, pseudocontact shifts (PCSs) arise and can be easily measured. They constitute very useful restraints for the solution structure determination of metal-binding proteins. In this context, there has been great interest in the use of lanthanide(III) ions to induce PCSs in diamagnetic proteins, e.g. through the replacement native calcium(II) ions. By preparing multiple samples in each of which a different ion of the lanthanide series is introduced, it is possible to obtain multiple independent PCS datasets that can be used synergistically to generate protein structure ensembles (typically called bundles). For typical NMR-based determination of protein structure, it is necessary to perform an energetic refinement of such initial bundles to obtain final structures whose geometric quality is suitable for deposition in the PDB. This can be conveniently done by using restrained molecular dynamics simulations (rMD) in explicit solvent. However, there are no available protocols for rMD using multiple PCS datasets as part of the restraints. In this work, we extended the PCS module of the AMBER MD package to handle multiple datasets and tuned a previously developed protocol for NMR structure refinement to achieve consistent convergence with PCS restraints. Test calculations with real experimental data show that this new implementation delivers the expected improvement of protein geometry, resulting in final structures that are of suitable quality for deposition. Furthermore, we observe that also initial structures generated only with traditional restraints can be successfully refined using traditional and PCS restraints simultaneously.

G-triplex structure and formation propensity.

The occurrence of a G-triplex folding intermediate of thrombin binding aptamer (TBA) has been recently predicted by metadynamics calculations, and experimentally supported by Nuclear Magnetic Resonance (NMR), Circular Dichroism (CD) and Differential Scanning Calorimetry (DSC) data collected on a 3' end TBA-truncated 11-mer oligonucleotide (11-mer-3'-t-TBA). Here we present the solution structure of 11-mer-3'-t-TBA in the presence of potassium ions. This structure is the first experimental example of a G-triplex folding, where a network of Hoogsteen-like hydrogen bonds stabilizes six guanines to form two G:G:G triad planes. The G-triplex folding of 11-mer-3'-t-TBA is stabilized by the potassium ion and destabilized by increasing the temperature. The superimposition of the experimental structure with that predicted by metadynamics shows a great similarity, with only significant differences involving two loops. These new structural data show that 11-mer-3'-t-TBA assumes a G-triplex DNA conformation as its stable form, reinforcing the idea that G-triplex folding intermediates may occur in vivo in human guanine-rich sequences. NMR and CD screening of eight different constructs obtained by removing from one to four bases at either the 3' and the 5' ends show that only the 11-mer-3'-t-TBA yields a relatively stable G-triplex.

Molecular view of an electron transfer process essential for iron-sulfur protein biogenesis.

Biogenesis of iron-sulfur cluster proteins is a highly regulated process that requires complex protein machineries. In the cytosolic iron-sulfur protein assembly machinery, two human key proteins--NADPH-dependent diflavin oxidoreductase 1 (Ndor1) and anamorsin--form a stable complex in vivo that was proposed to provide electrons for assembling cytosolic iron-sulfur cluster proteins. The Ndor1-anamorsin interaction was also suggested to be implicated in the regulation of cell survival/death mechanisms. In the present work we unravel the molecular basis of recognition between Ndor1 and anamorsin and of the electron transfer process. This is based on the structural characterization of the two partner proteins, the investigation of the electron transfer process, and the identification of those protein regions involved in complex formation and those involved in electron transfer. We found that an unstructured region of anamorsin is essential for the formation of a specific and stable protein complex with Ndor1, whereas the C-terminal region of anamorsin, containing the [2Fe-2S] redox center, transiently interacts through complementary charged residues with the FMN-binding site region of Ndor1 to perform electron transfer. Our results propose a molecular model of the electron transfer process that is crucial for understanding the functional role of this interaction in human cells.

Long-range correlated dynamics in intrinsically disordered proteins.

Intrinsically disordered proteins (IDPs) are involved in a wide variety of physiological and pathological processes and are best described by ensembles of rapidly interconverting conformers. Using fast field cycling relaxation measurements we here show that the IDP α-synuclein as well as a variety of other IDPs undergoes slow reorientations at time scales comparable to folded proteins. The slow motions are not perturbed by mutations in α-synuclein, which are related to genetic forms of Parkinson's disease, and do not depend on secondary and tertiary structural propensities. Ensemble-based hydrodynamic calculations suggest that the time scale of the underlying correlated motion is largely determined by hydrodynamic coupling between locally rigid segments. Our study indicates that long-range correlated dynamics are an intrinsic property of IDPs and offers a general physical mechanism of correlated motions in highly flexible biomolecular systems.

Mechanistic insights into the superoxide-cytochrome c reaction by lysine surface scanning.

This study summarizes results which have been obtained by a mutational study of human cytochrome c. The protein can be used as a recognition element in analytical assays and biosensors for superoxide radicals since ferricytochrome c reacts with superoxide to form ferrocytochrome c and oxygen. Here lysine mutagenesis of the distal surface (i.e., of exposed residues around the Met80 axial ligand) of human cytochrome c was pursued to evaluate the effect of the surface charges on the reaction rate with the superoxide anion radical and on the redox properties of the heme protein. The latter feature is particularly relevant when the protein is immobilized on a negatively charged self-assembled monolayer on an electrode to be used as a biosensor. The observed effects of the mutations are rationalized through structural investigations based on solution NMR spectroscopy and computational analysis of the surface electrostatics. The results suggest the presence of a specific path that guides superoxide toward an efficient reaction site. Localized positive charges at the rim of the entry channel are effective in increasing the reaction rate, whereas diffused positive charges or charges far from this area are not effective or are even detrimental, resulting in a misguided approach of the anion to the protein surface.

Cytochrome c and superoxide: a reply.

In his comments, W.H. Koppenol criticizes our article with respect to our conclusions and procedures. In this answer, we respond in detail to his objections, demonstrating that the approaches used are commonly accepted in the literature and that he makes a number of assumptions regarding our proposed mechanism that are not justified. Our study is thus a contribution to the ongoing investigation of the behavior of cytochrome c.

A Grid-enabled web portal for NMR structure refinement with AMBER.

MOTIVATION: The typical workflow for NMR structure determination involves collecting thousands of conformational restraints, calculating a bundle of 20-40 conformers in agreement with them and refining the energetics of these conformers. The structure calculation step employs simulated annealing based on molecular dynamics (MD) simulations with very simplified force fields. The value of refining the calculated conformers using restrained MD (rMD) simulations with state-of-art force fields is documented. This refinement however presents various subtleties, from the proper formatting of conformational restraints to the definition of suitable protocols. RESULTS: We describe a web interface to set up and run calculations with the AMBER package, which we called AMPS-NMR (AMBER-based Portal Server for NMR structures). The interface allows the refinement of NMR structures through rMD. Some predefined protocols are provided for this purpose, which can be personalized; it is also possible to create an entirely new protocol. AMPS-NMR can handle various restraint types. Standard rMD refinement in explicit water of the structures of three different proteins are shown as examples. AMPS-NMR additionally includes a workspace for the user to store different calculations. As an ancillary service, a web interface to AnteChamber is available, enabling the calculation of force field parameters for organic molecules such as ligands in protein-ligand adducts. AVAILABILITY AND IMPLEMENTATION: AMPS-NMR is embedded within the NMR services of the WeNMR project and is available at http://py-enmr.cerm.unifi.it/access/index/amps-nmr; its use requires registration with a digital certificate. CONTACT: ivanobertini@cerm.unifi.it SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Molecular Structure of Cu(II)-Bound Amyloid-ß Monomer Implicated in Inhibition of Peptide Self-Assembly in Alzheimer's Disease.

Metal ions, such as copper and zinc ions, have been shown to strongly modulate the self-assembly of the amyloid-ß (Aß) peptide into insoluble fibrils, and elevated concentrations of metal ions have been found in amyloid plaques of Alzheimer's patients. Among the physiological transition metal ions, Cu(II) ions play an outstanding role since they can trigger production of neurotoxic reactive oxygen species. In contrast, structural insights into Cu(II) coordination of Aß have been challenging due to the paramagnetic nature of Cu(II). Here, we employed specifically tailored paramagnetic NMR experiments to determine NMR structures of Cu(II) bound to monomeric Aß. We found that monomeric Aß binds Cu(II) in the N-terminus and combined with molecular dynamics simulations, we could identify two prevalent coordination modes of Cu(II). For these, we report here the NMR structures of the Cu(II)-bound Aß complex, exhibiting heavy backbone RMSD values of 1.9 and 2.1 Å, respectively. Further, applying aggregation kinetics assays, we identified the specific effect of Cu(II) binding on the Aß nucleation process. Our results show that Cu(II) efficiently retards Aß fibrillization by predominately reducing the rate of fibril-end elongation at substoichiometric ratios. A detailed kinetic analysis suggests that this specific effect results in enhanced Aß oligomer generation promoted by Cu(II). These results can quantitatively be understood by Cu(II) interaction with the Aß monomer, forming an aggregation inert complex. In fact, this mechanism is strikingly similar to other transition metal ions, suggesting a common mechanism of action of retarding Aß self-assembly, where the metal ion binding to monomeric Aß is a key determinant.

Structural Biology in the Clouds: The WeNMR-EOSC Ecosystem.

Structural biology aims at characterizing the structural and dynamic properties of biological macromolecules at atomic details. Gaining insight into three dimensional structures of biomolecules and their interactions is critical for understanding the vast majority of cellular processes, with direct applications in health and food sciences. Since 2010, the WeNMR project (www.wenmr.eu) has implemented numerous web-based services to facilitate the use of advanced computational tools by researchers in the field, using the high throughput computing infrastructure provided by EGI. These services have been further developed in subsequent initiatives under H2020 projects and are now operating as Thematic Services in the European Open Science Cloud portal (www.eosc-portal.eu), sending >12 millions of jobs and using around 4,000 CPU-years per year. Here we review 10 years of successful e-infrastructure solutions serving a large worldwide community of over 23,000 users to date, providing them with user-friendly, web-based solutions that run complex workflows in structural biology. The current set of active WeNMR portals are described, together with the complex backend machinery that allows distributed computing resources to be harvested efficiently.

Conformational space of flexible biological macromolecules from average data.

The concept of maximum occurrence (MO), i.e., the maximum percent of time that flexible proteins can spend in any given conformation, is introduced, and a rigorous method is developed to extensively sample the conformational space and to construct MO maps from experimental data. The method is tested in a case study, the flexible two-domain protein calmodulin (CaM), using SAXS and NMR data (i.e., pseudocontact shifts and self-orientation residual dipolar couplings arising from the presence of paramagnetic lanthanide ions), revealing that the "closed" and "fully extended" conformations trapped in the crystalline forms of CaM have MOs of only 5 and 15%, respectively. Compact conformations in general have small MOs, whereas some extended conformations have MO as high as 35%, strongly suggesting these conformations to be most abundant in solution. The method is universally applicable as it requires only standard SAXS data and specific NMR data on lanthanide derivatives of the protein (using native metal sites or lanthanide tagging). The computer program is publicly available using the grid computing infrastructure through the authors' Web portal.

A protocol to automatically calculate homo-oligomeric protein structures through the integration of evolutionary constraints and NMR ambiguous contacts.

Protein assemblies are involved in many important biological processes. Solid-state NMR (SSNMR) spectroscopy is a technique suitable for the structural characterization of samples with high molecular weight and thus can be applied to such assemblies. A significant bottleneck in terms of both effort and time required is the manual identification of unambiguous intermolecular contacts. This is particularly challenging for homo-oligomeric complexes, where simple uniform labeling may not be effective. We tackled this challenge by exploiting coevolution analysis to extract information on homo-oligomeric interfaces from NMR-derived ambiguous contacts. After removing the evolutionary couplings (ECs) that are already satisfied by the 3D structure of the monomer, the predicted ECs are matched with the automatically generated list of experimental contacts. This approach provides a selection of potential interface residues that is used directly in monomer-monomer docking calculations. We validated the protocol on tetrameric L-asparaginase II and dimeric Sod1.

An atomistic view of the YiiP structural changes upon zinc(II) binding.

BACKGROUND: YiiP is a bacterial zinc-for-proton antiporter belonging to the cation diffusion facilitator family. The zinc(II) ions are transported across the cell membrane, from the cytosol to the extracellular space. METHODS: We performed atomistic molecular dynamics simulations of the YiiP dimer with zinc(II) ions in solution to elucidate how the metal ions interact with the protein while moving from the cytosol to the transport site. RESULTS: We observed that of the two cavities of the dimer, only one was accessible from the cytosol during transport. Zinc(II) binding to D49 of the transport site triggered a rearrangement of the transmembrane domain that closed the accessible cavity. Finally, we analyzed the free-energy profiles of metal transit in the channel and observed the existence of a high barrier preventing release from the transport site. CONCLUSIONS: The observed dynamics is consistent with the dimer-dimer interface forming a stable scaffold against which the rest of the trans-membrane rearranges. GENERAL SIGNIFICANCE: Zinc(II) transporters are present in all kingdoms of life. The present study highlights structural features that might be of general relevance.

Conformational characterization of full-length X-chromosome-linked inhibitor of apoptosis protein (XIAP) through an integrated approach.

The X-chromosome-linked inhibitor of apoptosis protein (XIAP) is a multidomain protein whose main function is to block apoptosis by caspase inhibition. XIAP is also involved in other signalling pathways, including NF-κB activation and copper homeostasis. XIAP is overexpressed in tumours, potentiating cell survival and resistance to chemotherapeutics, and has therefore become an important target for the treatment of malignancy. Despite the fact that the structure of each single domain is known, the conformation of the full-length protein has never been determined. Here, the first structural model of the full-length XIAP dimer, determined by an integrated approach using nuclear magnetic resonance, small-angle X-ray scattering and electron paramagnetic resonance data, is presented. It is shown that XIAP adopts a compact and relatively rigid conformation, implying that the spatial arrangement of its domains must be taken into account when studying the interactions with its physiological partners and in developing effective inhibitors.

West-Life: A Virtual Research Environment for structural biology.

The West-Life project (https://about.west-life.eu/) is a Horizon 2020 project funded by the European Commission to provide data processing and data management services for the international community of structural biologists, and in particular to support integrative experimental approaches within the field of structural biology. It has developed enhancements to existing web services for structure solution and analysis, created new pipelines to link these services into more complex higher-level workflows, and added new data management facilities. Through this work it has striven to make the benefits of European e-Infrastructures more accessible to life-science researchers in general and structural biologists in particular.