Gating topology of the proton-coupled oligopeptide symporters.

Proton-coupled oligopeptide transporters belong to the major facilitator superfamily (MFS) of membrane transporters. Recent crystal structures suggest the MFS fold facilitates transport through rearrangement of their two six-helix bundles around a central ligand binding site; how this is achieved, however, is poorly understood. Using modeling, molecular dynamics, crystallography, functional assays, and site-directed spin labeling combined with double electron-electron resonance (DEER) spectroscopy, we present a detailed study of the transport dynamics of two bacterial oligopeptide transporters, PepTSo and PepTSt. Our results identify several salt bridges that stabilize outward-facing conformations and we show that, for all the current structures of MFS transporters, the first two helices of each of the four inverted-topology repeat units form half of either the periplasmic or cytoplasmic gate and that these function cooperatively in a scissor-like motion to control access to the peptide binding site during transport.

Constraint Network Analysis (CNA): a Python software package for efficiently linking biomacromolecular structure, flexibility, (thermo-)stability, and function.

For deriving maximal advantage from information on biomacromolecular flexibility and rigidity, results from rigidity analyses must be linked to biologically relevant characteristics of a structure. Here, we describe the Python-based software package Constraint Network Analysis (CNA) developed for this task. CNA functions as a front- and backend to the graph-based rigidity analysis software FIRST. CNA goes beyond the mere identification of flexible and rigid regions in a biomacromolecule in that it (I) provides a refined modeling of thermal unfolding simulations that also considers the temperature-dependence of hydrophobic tethers, (II) allows performing rigidity analyses on ensembles of network topologies, either generated from structural ensembles or by using the concept of fuzzy noncovalent constraints, and (III) computes a set of global and local indices for quantifying biomacromolecular stability. This leads to more robust results from rigidity analyses and extends the application domain of rigidity analyses in that phase transition points ("melting points") and unfolding nuclei ("structural weak spots") are determined automatically. Furthermore, CNA robustly handles small-molecule ligands in general. Such advancements are important for applying rigidity analysis to data-driven protein engineering and for estimating the influence of ligand molecules on biomacromolecular stability. CNA maintains the efficiency of FIRST such that the analysis of a single protein structure takes a few seconds for systems of several hundred residues on a single core. These features make CNA an interesting tool for linking biomacromolecular structure, flexibility, (thermo-)stability, and function. CNA is available from http://cpclab.uni-duesseldorf.de/software for nonprofit organizations.

Identification of molecular hinge points mediating alternating access in the vesicular monoamine transporter VMAT2.

Vesicular monoamine transporter 2 (VMAT2) catalyzes transport of monoamines into storage vesicles in a process that involves exchange of the charged monoamine with two protons. VMAT2 is a member of the DHA12 family of multidrug transporters that belongs to the major facilitator superfamily (MFS) of secondary transporters. Here we present a homology model of VMAT2, which has the standard MFS fold, that is, with two domains of six transmembrane helices each which are related by twofold pseudosymmetry and whose axis runs normal to the membrane and between the two halves. Demonstration of the essential role of a membrane-embedded glutamate and confirmation of the existence of a hydrogen bond probably involved in proton transport provide experimental evidence that validates some of the predictions inherent to the model. Moreover, we show the essential role of residues at two anchor points between the two bundles. These residues appear to function as molecular hinge points about which the two six transmembrane-helix bundles flex and straighten to open and close the pathways on either side of the membrane as required for transport. Polar residues that create a hydrogen bond cluster form one of the anchor points of VMAT2. The other results from hydrophobic interactions. Residues at the anchor points are strongly conserved in other MFS transporters in one way or another, suggesting that interactions at these locations will be critical in most, if not all, MFS transporters.

Global and local indices for characterizing biomolecular flexibility and rigidity.

Understanding flexibility and rigidity characteristics of biomolecules is a prerequisite for understanding biomolecular structural stability and function. Computational methods have been implemented that directly characterize biomolecular flexibility and rigidity by constraint network analysis. For deriving maximal advantage from these analyses, their results need to be linked to biologically relevant characteristics of a structure. Such links are provided by global and local measures ("indices") of biomolecular flexibility and rigidity. To date, more than 14 indices are available with sometimes overlapping or only vague definitions. We present concise definitions of these indices, analyze the relation between, and the scope and limitations of them, and compare their informative value. For this, we probe the structural stability of the calcium binding protein α-lactalbumin as a showcase, both in the "ground state" and after perturbing the system by changing the network topology. In addition, we introduce three indices for the first time that extend the application domain of flexibility and rigidity analyses. The results allow us to provide guidelines for future studies suggesting which of these indices could best be used for analyzing, understanding, and quantifying structural features that are important for biomolecular stability and function. Finally, we make suggestions for proper index notations in future studies to prevent the misinterpretation and to facilitate the comparison of results obtained from flexibility and rigidity analyses.

Thermostabilizing mutations preferentially occur at structural weak spots with a high mutation ratio.

We apply Constraint Network Analysis (CNA) to investigate the relationship between structural rigidity and thermostability of five citrate synthase (CS) structures over a temperature range from 37 °C to 100 °C. For the first time, we introduce an ensemble-based variant of CNA and model the temperature-dependence of hydrophobic interactions in the constraint network. A very good correlation between the predicted thermostabilities of CS and optimal growth temperatures of their source organisms (R²=0.88, p=0.017) is obtained, which validates that CNA is able to quantitatively discriminate between less and more thermostable proteins even within a series of orthologs. Structural weak spots on a less thermostable CS, predicted by CNA to be in the top 5% with respect to the frequency of occurrence over an ensemble, have a higher mutation ratio in a more thermostable CS than other sequence positions. Furthermore, highly ranked weak spots that are also highly conserved with respect to the amino acid type found at that sequence position are nevertheless found to be mutated in the more stable CS. As for mechanisms at an atomic level that lead to a reinforcement of weak spots in more stable CS, we observe that the thermophilic CS achieve a higher thermostability by better hydrogen bonding networks whereas hyperthermophilic CS incorporate more hydrophobic contacts to reach the same goal. Overall, these findings suggest that CNA can be applied as a pre-filter in data-driven protein engineering to focus on residues that are highly likely to improve thermostability upon mutation.

The alternating-access mechanism of MFS transporters arises from inverted-topology repeats.

Lactose permease (LacY) is the prototype of the major facilitator superfamily (MFS) of secondary transporters. Available structures of LacY reveal a state in which the substrate is exposed to the cytoplasm but is occluded from the periplasm. However, the alternating-access transport mechanism requires the existence of a periplasm-facing state. We recently showed that inverted-topology structural repeats provide the foundation for the mechanisms of two transporter families with folds distinct from the MFS. Here, we generated a structural model of LacY by swapping the conformations of inverted-topology repeats identified in its two domains. The model exhibits all required properties of an outward-facing conformation, i.e., closure of the binding site to the cytoplasm and exposure to the periplasm. Furthermore, the model agrees with double electron-electron resonance distance changes, accessibility to cysteine-modifying reagents, cysteine cross-linking data, and a recent structure of a distantly related transporter. Analysis of the intradomain differences between the two states suggests a role for conserved sequence motifs in occluding the central pathway through kinking of the pore-lining helices. In addition, predicted re-pairing of critical salt-bridging residues in the binding sites agrees remarkably well with previous proposals, allowing a description of the proton/sugar transport mechanism. More fundamentally, our model demonstrates that inverted-topology repeats provide the foundation for the alternating-access mechanisms of MFS transporters.

Protein rigidity and thermophilic adaptation.

We probe the hypothesis of corresponding states, according to which homologues from mesophilic and thermophilic organisms are in corresponding states of similar rigidity and flexibility at their respective optimal temperatures. For this, the local distribution of flexible and rigid regions in 19 pairs of homologous proteins from meso- and thermophilic organisms is analyzed and related to activity characteristics of the enzymes by constraint network analysis (CNA). Two pairs of enzymes are considered in more detail: 3-isopropylmalate dehydrogenase and thermolysin-like protease. By comparing microscopic stability features of homologues with the help of stability maps, introduced for the first time, we show that adaptive mutations in enzymes from thermophilic organisms maintain the balance between overall rigidity, important for thermostability, and local flexibility, important for activity, at the appropriate working temperature. Thermophilic adaptation in general leads to an increase of structural rigidity but conserves the distribution of functionally important flexible regions between homologues. This finding provides direct evidence for the hypothesis of corresponding states. CNA thereby implicitly captures and unifies many different mechanisms that contribute to increased thermostability and to activity at high temperatures. This allows to qualitatively relate changes in the flexibility of active site regions, induced either by a temperature change or by the introduction of mutations, to experimentally observed losses of the enzyme function. As for applications, the results demonstrate that exploiting the principle of corresponding states not only allows for successful thermostability optimization but also for guiding experiments in order to improve enzyme activity in protein engineering.

Homology model-based virtual screening for GPCR ligands using docking and target-biased scoring.

The current study investigates the combination of two recently reported techniques for the improvement of homology model-based virtual screening for G-protein coupled receptor (GPCR) ligands. First, ligand-supported homology modeling was used to generate receptor models that were in agreement with mutagenesis data and structure-activity relationship information of the ligands. Second, interaction patterns from known ligands to the receptor were applied for scoring and rank ordering compounds from a virtual library using ligand-receptor interaction fingerprint-based similarity (IFS). Our approach was evaluated in retrospective virtual screening experiments for antagonists of the metabotropic glutamate receptor (mGluR) subtype 5. The results of our approach were compared to the results obtained by conventional scoring functions (Dock-Score, PMF-Score, Gold-Score, ChemScore, and FlexX-Score). The IFS lead to significantly higher enrichment rates, relative to the competing scoring functions. Though using a target-biased scoring approach, the results were not biased toward the chemical classes of the reference structures. Our results indicate that the presented approach has the potential to serve as a general setup for successful structure-based GPCR virtual screening.

Maximum common binding modes (MCBM): consensus docking scoring using multiple ligand information and interaction fingerprints.

Improving the scoring functions for small molecule-protein docking is a highly challenging task in current computational drug design. Here we present a novel consensus scoring concept for the prediction of binding modes for multiple known active ligands. Similar ligands are generally believed to bind to their receptor in a similar fashion. The presumption of our approach was that the true binding modes of similar ligands should be more similar to each other compared to false positive binding modes. The number of conserved (consensus) interactions between similar ligands was used as a docking score. Patterns of interactions were modeled using ligand receptor interaction fingerprints. Our approach was evaluated for four different data sets of known cocrystal structures (CDK-2, dihydrofolate reductase, HIV-1 protease, and thrombin). Docking poses were generated with FlexX and rescored by our approach. For comparison the CScore scoring functions from Sybyl were used, and consensus scores were calculated thereof. Our approach performed better than individual scoring functions and was comparable to consensus scoring. Analysis of the distribution of docking poses by self-organizing maps (SOM) and interaction fingerprints confirmed that clusters of docking poses composed of multiple ligands were preferentially observed near the native binding mode. Being conceptually unrelated to commonly used docking scoring functions our approach provides a powerful method to complement and improve computational docking experiments.

Improving binding mode predictions by docking into protein-specifically adapted potential fields.

The development of a protein-specifically adapted objective function for docking is described. Structural and energetic information about known protein-ligand complexes is exploited to tailor knowledge-based potentials using a "reverse", protein-based CoMFA-type (=AFMoC) approach. That way, effects due to protein flexibility and information about multiple solvation schemes can be implicitly incorporated. Compared to the application of AFMoC for binding affinity predictions, a Shannon entropy based column filtering of the descriptor matrix and the capping of adapted repulsive potentials within the binding site have turned out to be crucial for the success of this method. The new developed approach (AFMoC(obj)) was validated on a data set of 66 HIV-1 protease inhibitors, for which experimental structural information was available. Convincingly, for ligands with up to 20 rotatable bonds, in more than 75% of all cases a binding mode below 2 A rmsd has been identified on the first scoring rank when AFMoC(obj)-based potentials were used as the objective function in AutoDock. With respect to nonadapted DrugScore or AutoDock fields, the binding mode prediction accuracy was significantly improved by 14%. Noteworthy, very similar results were obtained for training and test set compounds, demonstrating the strength and robustness of this method. Implications of our findings for binding affinity predictions and its usage in virtual screening are further discussed.