Massimiliano Bonomi

Crystal structure of a eukaryotic phosphate transporter

Phosphate is crucial for structural and metabolic needs, including nucleotide and lipid synthesis, signalling and chemical energy storage. Proton-coupled transporters of the major facilitator superfamily (MFS) are essential for phosphate uptake in plants and fungi, and also have a function in sensing external phosphate levels as transceptors. Here we report the 2.9u2009Å structure of a fungal (Piriformospora indica) high-affinity phosphate transporter, PiPT, in an inward-facing occluded state, with bound phosphate visible in the membrane-buried binding site. The structure indicates both proton and phosphate exit pathways and suggests a modified asymmetrical ‘rocker-switch’ mechanism of phosphate transport. PiPT is related to several human transporter families, most notably the organic cation and anion transporters of the solute carrier family (SLC22), which are implicated in cancer-drug resistance. We modelled representative cation and anion SLC22 transporters based on the PiPT structure to surmise the structural basis for substrate binding and charge selectivity in this important family. The PiPT structure demonstrates and expands on principles of substrate transport by the MFS transporters and illuminates principles of phosphate uptake in particular.

Molecular basis of cyclooxygenase enzymes (COXs) selective inhibition

The widely used nonsteroidal anti-inflammatory drugs block the cyclooxygenase enzymes (COXs) and are clinically used for the treatment of inflammation, pain, and cancers. A selective inhibition of the different isoforms, particularly COX-2, is desirable, and consequently a deeper understanding of the molecular basis of selective inhibition is of great demand. Using an advanced computational technique we have simulated the full dissociation process of a highly potent and selective inhibitor, SC-558, in both COX-1 and COX-2. We have found a previously unreported alternative binding mode in COX-2 explaining the time-dependent inhibition exhibited by this class of inhibitors and consequently their long residence time inside this isoform. Our metadynamics-based approach allows us to illuminate the highly dynamical character of the ligand/protein recognition process, thus explaining a wealth of experimental data and paving the way to an innovative strategy for designing new COX inhibitors with tuned selectivity.

The Unfolded Ensemble and Folding Mechanism of the C-Terminal GB1 β-Hairpin

In this work, we shed new light on a much-studied case of beta-hairpin folding by means of advanced molecular dynamics simulations. A fully atomistic description of the protein and the solvent molecule is used, together with metadynamics, to accelerate the sampling and estimate free-energy landscapes. This is achieved using the path collective variables approach, which provides an adaptive description of the mechanism under study. We discover that the folding mechanism is a multiscale process where the turn region conformation leads to two different energy pathways that are connected by elongated structures. The former displays a stable 2:4 native-like structure in which an optimal hydrophobic packing and hydrogen bond pattern leads to 8 kcal/mol of stabilization. The latter shows a less-structured 3:5 beta-sheet, where hydrogen bonds and hydrophobic packing provide only 2.5 kcal/mol of stability. This perspective is fully consistent with experimental evidence that shows this to be a prototypical two-state folder, while it redefines the nature of the unfolded state.

Accuracy of Current All-Atom Force-Fields in Modeling Protein Disordered States

Molecular Dynamics (MD) plays a fundamental role in characterizing protein disordered states that are emerging as crucial actors in many biological processes. Here we assess the accuracy of three current force-fields in modeling disordered peptides by combining enhanced-sampling MD simulations with NMR data. These force-fields generate significantly different conformational ensembles, and AMBER03w [ Best and Mittal J. Phys. Chem. B 2010 , 114 , 14916 - 14923 ] provides the best agreement with experiments, which is further improved by adding the ILDN corrections [ Lindorff-Larsen et al. Proteins 2010 , 78 , 1950 - 1958 ].

Efficient Simulation of Explicitly Solvated Proteins in the Well-Tempered Ensemble

Herein, we report significant reduction in the cost of combined parallel tempering and metadynamics simulations (PTMetaD). The efficiency boost is achieved using the recently proposed well-tempered ensemble (WTE) algorithm. We studied the convergence of PTMetaD-WTE conformational sampling and free energy reconstruction of an explicitly solvated 20-residue tryptophan-cage protein (trp-cage). A set of PTMetaD-WTE simulations was compared to a corresponding standard PTMetaD simulation. The properties of PTMetaD-WTE and the convergence of the calculations were compared. The roles of the number of replicas, total simulation time, and adjustable WTE parameter γ were studied.

Structural basis for alternating access of a eukaryotic calcium/proton exchanger

Eukaryotic Ca2+ regulation involves sequestration into intracellular organelles, and expeditious Ca2+ release into the cytosol is a hallmark of key signalling transduction pathways. Bulk removal of Ca2+ after such signalling events is accomplished by members of the Ca2+:cation (CaCA) superfamily. The CaCA superfamily includes the Na+/Ca2+ (NCX) and Ca2+/H+ (CAX) antiporters, and in mammals the NCX and related proteins constitute families SLC8 and SLC24, and are responsible for the re-establishment of Ca2+ resting potential in muscle cells, neuronal signalling and Ca2+ reabsorption in the kidney. The CAX family members maintain cytosolic Ca2+ homeostasis in plants and fungi during steep rises in intracellular Ca2+ due to environmental changes, or following signal transduction caused by events such as hyperosmotic shock, hormone response and response to mating pheromones. The cytosol-facing conformations within the CaCA superfamily are unknown, and the transport mechanism remains speculative. Here we determine a crystal structure of the Saccharomyces cerevisiae vacuolar Ca2+/H+ exchanger (Vcx1) at 2.3u2009Å resolution in a cytosol-facing, substrate-bound conformation. Vcx1 is the first structure, to our knowledge, within the CAX family, and it describes the key cytosol-facing conformation of the CaCA superfamily, providing the structural basis for a novel alternating access mechanism by which the CaCA superfamily performs high-throughput Ca2+ transport across membranes.

Metainference: A Bayesian inference method for heterogeneous systems

Researchers present a Bayesian inference method for heterogeneous systems that integrates prior information with noisy experimental data. Modeling a complex system is almost invariably a challenging task. The incorporation of experimental observations can be used to improve the quality of a model and thus to obtain better predictions about the behavior of the corresponding system. This approach, however, is affected by a variety of different errors, especially when a system simultaneously populates an ensemble of different states and experimental data are measured as averages over such states. To address this problem, we present a Bayesian inference method, called “metainference,” that is able to deal with errors in experimental measurements and with experimental measurements averaged over multiple states. To achieve this goal, metainference models a finite sample of the distribution of models using a replica approach, in the spirit of the replica-averaging modeling based on the maximum entropy principle. To illustrate the method, we present its application to a heterogeneous model system and to the determination of an ensemble of structures corresponding to the thermal fluctuations of a protein molecule. Metainference thus provides an approach to modeling complex systems with heterogeneous components and interconverting between different states by taking into account all possible sources of errors.

Principles of protein structural ensemble determination.

The biological functions of protein molecules are intimately dependent on their conformational dynamics. This aspect is particularly evident for disordered proteins, which constitute perhaps one-third of the human proteome. Therefore, structural ensembles often offer more useful representations of proteins than individual conformations. Here, we describe how the well-established principles of protein structure determination should be extended to the case of protein structural ensembles determination. These principles concern primarily how to deal with conformationally heterogeneous states, and with experimental measurements that are averaged over such states and affected by a variety of errors. We first review the growing literature of recent methods that combine experimental and computational information to model structural ensembles, highlighting their similarities and differences. We then address some conceptual problems in the determination of structural ensembles and define future goals towards the establishment of objective criteria for the comparison, validation, visualization and dissemination of such ensembles.

Ligand Discovery for the Alanine-Serine-Cysteine Transporter (ASCT2, SLC1A5) from Homology Modeling and Virtual Screening

The Alanine-Serine-Cysteine transporter ASCT2 (SLC1A5) is a membrane protein that transports neutral amino acids into cells in exchange for outward movement of intracellular amino acids. ASCT2 is highly expressed in peripheral tissues such as the lung and intestines where it contributes to the homeostasis of intracellular concentrations of neutral amino acids. ASCT2 also plays an important role in the development of a variety of cancers such as melanoma by transporting amino acid nutrients such as glutamine into the proliferating tumors. Therefore, ASCT2 is a key drug target with potentially great pharmacological importance. Here, we identify seven ASCT2 ligands by computational modeling and experimental testing. In particular, we construct homology models based on crystallographic structures of the aspartate transporter GltPh in two different conformations. Optimization of the models’ binding sites for protein-ligand complementarity reveals new putative pockets that can be targeted via structure-based drug design. Virtual screening of drugs, metabolites, fragments-like, and lead-like molecules from the ZINC database, followed by experimental testing of 14 top hits with functional measurements using electrophysiological methods reveals seven ligands, including five activators and two inhibitors. For example, aminooxetane-3-carboxylate is a more efficient activator than any other known ASCT2 natural or unnatural substrate. Furthermore, two of the hits inhibited ASCT2 mediated glutamine uptake and proliferation of a melanoma cancer cell line. Our results improve our understanding of how substrate specificity is determined in amino acid transporters, as well as provide novel scaffolds for developing chemical tools targeting ASCT2, an emerging therapeutic target for cancer and neurological disorders.

Efficient Sampling of High-Dimensional Free-Energy Landscapes with Parallel Bias Metadynamics.

Metadynamics accelerates sampling of molecular dynamics while reconstructing thermodynamic properties of selected descriptors of the system. Its main practical difficulty originates from the compromise between keeping the number of descriptors small for efficiently exploring their multidimensional free-energy landscape and biasing all of the slow motions of a process. Here we illustrate on a model system and on the tryptophan-cage miniprotein parallel bias metadynamics, a method that overcomes this issue by simultaneously applying multiple low-dimensional bias potentials.