Charles H. Robert

Community-wide assessment of protein-interface modeling suggests improvements to design methodology

The CAPRI (Critical Assessment of Predicted Interactions) and CASP (Critical Assessment of protein Structure Prediction) experiments have demonstrated the power of community-wide tests of methodology in assessing the current state of the art and spurring progress in the very challenging areas of protein docking and structure prediction. We sought to bring the power of community-wide experiments to bear on a very challenging protein design problem that provides a complementary but equally fundamental test of current understanding of protein-binding thermodynamics. We have generated a number of designed protein-protein interfaces with very favorable computed binding energies but which do not appear to be formed in experiments, suggesting that there may be important physical chemistry missing in the energy calculations. A total of 28 research groups took up the challenge of determining what is missing: we provided structures of 87 designed complexes and 120 naturally occurring complexes and asked participants to identify energetic contributions and/or structural features that distinguish between the two sets. The community found that electrostatics and solvation terms partially distinguish the designs from the natural complexes, largely due to the nonpolar character of the designed interactions. Beyond this polarity difference, the community found that the designed binding surfaces were, on average, structurally less embedded in the designed monomers, suggesting that backbone conformational rigidity at the designed surface is important for realization of the designed function. These results can be used to improve computational design strategies, but there is still much to be learned; for example, one designed complex, which does form in experiments, was classified by all metrics as a nonbinder.

Pore formation kinetics in membranes, determined from the release of marker molecules out of liposomes or cells

We discuss the efflux of entrapped marker material from liposomes or cells through pores in the membrane, being monitored by the time course of a certain signal F (e.g., fluorescence emission). This is expressed in terms of an appropriate normalized function of time, the so-called efflux function E(t). Under conditions frequently encountered in practice the measured E(t) can be easily related to the forward rate of pore formation if the liposomes/cells are monodisperse in size. In the basic case of a time-independent rate law it turns out that E(t) must be single exponential. Deviations from such a simple functional behavior might be due to a fairly broad distribution of liposome/cell sizes and/or a more complicated pore formation mechanism. A relevant evaluation of original data is demonstrated making use of experimental results obtained with small unilamellar lipid vesicles where pores are induced by the antibiotic peptide alamethicin. This includes the application of a general method to eliminate the effect of a given liposome/cell size distribution.

Normal mode analysis as a prerequisite for drug design: Application to matrix metalloproteinases inhibitors

We demonstrate the utility of normal mode analysis in correctly predicting the binding modes of inhibitors in the active sites of matrix metalloproteinases (MMPs). We show the accuracy in predicting the positions of MMP‐3 inhibitors is strongly dependent on which structure is used as the target, especially when it has been energy minimized. This dependency can be overcome by using intermediate structures generated along one of the normal modes previously calculated for a given target. These results may be of prime importance for further in silico drug discovery.

New Insights into the Allosteric Mechanism of Human Hemoglobin from Molecular Dynamics Simulations

It is still difficult to obtain a precise structural description of the transition between the deoxy T-state and oxy R-state conformations of human hemoglobin, despite a large number of experimental studies. We used molecular dynamics with the Path Exploration with Distance Constraints (PEDC) method to provide new insights into the allosteric mechanism at the atomic level, by simulating the T-to-R transition. The T-state molecule in the absence of ligands was seen to have a natural propensity for dimer rotation, which nevertheless would be hampered by steric hindrance in the joint region. The binding of a ligand to the alpha subunit would prevent such hindrance due to the coupling between this region and the alpha proximal histidine, and thus facilitate completion of the dimer rotation. Near the end of this quaternary transition, the switch region adopts the R conformation, resulting in a shift of the beta proximal histidine. This leads to a sliding of the beta-heme, the effect of which is to open the beta-hemes distal side, increasing the accessibility of the Fe atom and thereby the affinity of the protein. Our simulations are globally consistent with the Perutz strereochemical mechanism.

Odorant Binding and Conformational Dynamics in the Odorant-binding Protein

In mammals, the olfactory epithelium secretes odorant-binding proteins (OBPs), which are lipocalins found freely dissolved in the mucus layer protecting the olfactory neurons. OBPs may act as passive transporters of predominantly hydrophobic odorant molecules across the aqueous mucus layer, or they may play a more active role in which the olfactory neuronal receptor recognizes the OBP-ligand complex. To better understand the molecular events accompanying the initial steps in the olfaction process, we have performed molecular dynamics studies of rat and pig OBPs with the odorant molecule thymol. These calculations provide an atomic level description of conformational changes and pathway intermediates that remain difficult to study directly. A series of eight independent molecular dynamics trajectories of rat OBP permitted the observation of a consensus pathway for ligand unbinding and the calculation of the potential of mean force (PMF) along this path. Titration microcalorimetry confirmed the specific binding of thymol to this protein with a strong hydrophobic component. In both rat and pig OBPs we observed lipocalin strand pair opening in the presence of ligand, consistent with potential roles of these proteins in olfactive receptor recognition.

Kinetics of pore-mediated release of marker molecules from liposomes or cells

We present a general mathematical treatment of marker efflux from liposomes or cells mediated by pore formation with the idea of using it in practice to obtain basic information about the underlying rates and mechanism. The approach encompasses permeation of molecules through any kind of pore-like defects in a cell membrane as they are induced by the action of some external agent. The approach broadens an earlier treatment to the more realistic general case in which a distribution of pore lifetimes must be taken into account. We derive a theoretical retention function describing the amount of marker remaining in the cells, formulated in terms of the pore activation and inactivation kinetics. The phenomenological efflux function evaluated directly from experimental data, is directly comparable with this retention function so long as the experimental signal is linearly related to the marker concentration. With the use of self-quenching dyes the relationship between signal and concentration is not, in general, linear, so that a more complicated treatment may be required. Even for these dyes, however, linearity holds under the frequently encountered condition of all-or-none release of dye from vesicles, a condition that can itself be verified experimentally.

A combined study of aggregation, membrane affinity and pore activity of natural and modified melittin

The pore activity of melittin and several chemically modified derivatives has been investigated using conductance measurements on planar lipid bilayers and marker release from small unilamellar vesicles. The modifications included N-terminal formylation, acetylation, succinylation and modification of the tryptophan residue. All of the compounds showed bilayer permeabilizing properties, though quantitative differences were evident. These comprised changes in the voltage dependence of the conductance, in the single-pore kinetics, in the concentration of aqueous peptide required to induce a given pore activity and in the apparent molecularity reflected by the power law of its concentration dependence. A strong tendency for disrupting bilayers was not always correlated with strong pore activity. For a better understanding of these results, measurements of pore activity were complemented by studying the aggregation behavior in solution and the water-membrane partition equilibrium. Modifications of charged residues gave rise to significant changes in the aggregation properties, had virtually no influence on the partition coefficient. The latter decreased strongly, however, as a result of tryptophan modification. Analysis of the isotherms was consistent with the assumption that the arginine residues in melittin do not contribute very much to charge accumulation at the immediate membrane/water interface.

How does heparin prevent the pH inactivation of cathepsin B? Allosteric mechanism elucidated by docking and molecular dynamics

BackgroundCathepsin B (catB) is a promising target for anti-cancer drug design due to its implication in several steps of tumorigenesis. catB activity and inhibition are pH-dependent, making it difficult to identify efficient inhibitor candidates for clinical trials. In addition it is known that heparin binding stabilizes the enzyme in alkaline conditions. However, the molecular mechanism of stabilization is not well understood, indicating the need for more detailed structural and dynamic studies in order to clarify the influence of pH and heparin binding on catB stability.ResultsOur pKa calculations of catB titratable residues revealed distinct protonation states under different pH conditions for six key residues, of which four lie in the crucial interdomain interface. This implies changes in the overall charge distribution at the catB surface, as revealed by calculation of the electrostatic potential. We identified two basic surface regions as possible heparin binding sites, which were confirmed by docking calculations. Molecular dynamics (MD) of both apo catB and catB-heparin complexes were performed using protonation states for catB residues corresponding to the relevant acidic or alkaline conditions. The MD of apo catB at pH 5.5 was very stable, and presented the highest number and occupancy of hydrogen bonds within the inter-domain interface. In contrast, under alkaline conditions the enzymes overall flexibility was increased: interactions between active site residues were lost, helical content decreased, and domain separation was observed as well as high-amplitude motions of the occluding loop – a main target of drug design studies. Essential dynamics analysis revealed that heparin binding modulates large amplitude motions promoting rearrangement of contacts between catB domains, thus favoring the maintenance of helical content as well as active site stability.ConclusionsThe results of our study contribute to unraveling the molecular events involved in catB inactivation in alkaline pH, highlighting the fact that protonation changes of few residues can alter the overall dynamics of an enzyme. Moreover, we propose an allosteric role for heparin in the regulation of catB stability in such a manner that the restriction of enzyme flexibility would allow the establishment of stronger contacts and thus the maintenance of overall structure.

Free Energy Profiles along Consensus Normal Modes Provide Insight into HIV-1 Protease Flap Opening

Describing biological macromolecular energetics from computer simulations can pose major challenges, and often necessitates enhanced conformational sampling. We describe the calculation of conformational free-energy profiles along carefully chosen collective coordinates: consensus normal modes, developed recently as robust alternatives to conventional normal modes. In an application to the HIV-1 protease, we obtain efficient sampling of significant flap opening movements governing inhibitor binding from relatively short simulations, in close correspondence with experimental results.

Side-chain rotamer transitions at protein-protein interfaces†

We compare the changes in side chain conformations that accompany the formation of protein–protein complexes, in residues forming either the interface or the remainder of the solvent‐accessible surface of the proteins in the Docking Benchmark 3.0. We find that the interface residues undergo significantly more changes than other surface residues, and these changes are more likely to convert them from a high‐energy torsion angle state to a lower‐energy one than the reverse. Moreover, in both the unbound proteins and the complexes, the interface residues are more frequently found to be in a high‐energy torsion angle state than the noninterface residues. As these differences exist before the binding step, they may be relevant to specificity and help in identifying binding sites for docking predictions. Proteins 2010.