Alchemical absolute protein-ligand binding free energies for drug design.

The recent advances in relative protein-ligand binding free energy calculations have shown the value of alchemical methods in drug discovery. Accurately assessing absolute binding free energies, although highly desired, remains a challenging endeavour, mostly limited to small model cases. Here, we demonstrate accurate first principles based absolute binding free energy estimates for 128 pharmaceutically relevant targets. We use a novel rigorous method to generate protein-ligand ensembles for the ligand in its decoupled state. Not only do the calculations deliver accurate protein-ligand binding affinity estimates, but they also provide detailed physical insight into the structural determinants of binding. We identify subtle rotamer rearrangements between apo and holo states of a protein that are crucial for binding. When compared to relative binding free energy calculations, obtaining absolute binding free energies is considerably more challenging in large part due to the need to explicitly account for the protein in its apo state. In this work we present several approaches to obtain apo state ensembles for accurate absolute ΔG calculations, thus outlining protocols for prospective application of the methods for drug discovery.

The anti-protozoal drug pentamidine blocks KIR2.x-mediated inward rectifier current by entering the cytoplasmic pore region of the channel.

BACKGROUND AND PURPOSE: Pentamidine is a drug used in treatment of protozoal infections. Pentamidine treatment may cause sudden cardiac death by provoking cardiac arrhythmias associated with QTc prolongation and U-wave alterations. This proarrhythmic effect was linked to inhibition of hERG trafficking, but not to acute block of ion channels contributing to the action potential. Because the U-wave has been linked to the cardiac inward rectifier current (I(K1)), we examined the action and mechanism of pentamidine-mediated I(K1) block. EXPERIMENTAL APPROACH: Patch clamp measurements of I(K1) were made on cultured adult canine ventricular cardiomyocytes, K(IR)2.1-HEK293 cells and K(IR)2.x inside-out patches. Pentamidine binding to cytoplasmic amino acid residues of K(IR)2.1 channels was studied by molecular modelling. KEY RESULTS: Pentamidine application (24 h) decreased I(K1) in cultured canine cardiomyocytes and K(IR)2.1-HEK293 cells under whole cell clamp conditions. Pentamidine inhibited I(K1) in K(IR)2.1-HEK293 cells 10 min after application. When applied to the cytoplasmic side under inside-out patch clamp conditions, pentamidine block of I(K1) was acute (IC(50)= 0.17 microM). Molecular modelling predicted pentamidine-channel interactions in the cytoplasmic pore region of K(IR)2.1 at amino acids E224, D259 and E299. Mutation of these conserved residues to alanine reduced pentamidine block of I(K1). Block was independent of the presence of spermine. K(IR)2.2, and K(IR)2.3 based I(K1) was also sensitive to pentamidine blockade. CONCLUSIONS AND IMPLICATIONS: Pentamidine inhibits cardiac I(K1) by interacting with three negatively charged amino acids in the cytoplasmic pore region. Our findings may provide new insights for development of specific I(K1) blocking compounds.

A plate reader-based method for cell water permeability measurement.

Cell volume and water permeability measurements in cultured mammalian cells are typically conducted under a light microscope. Many of the employed approaches are time consuming and not applicable to a study of confluent epithelial cell monolayers. We present here an adaptation of a calcein-quenching-based approach for a plate reader. A standard curve of fluorescence intensities at equilibrium has been recorded, following a shift from 285 mosmol/kgH(2)O to a series of altered extracellular osmolyte concentrations, ranging from final concentrations of 185 to 585 mosmol/kgH(2)O, by changing buffer d-mannitol concentrations. Similarly, according average cell volumes have been measured in suspension in a Coulter counter (particle-sizing device). Based on these measurements, we have derived an equation that facilitates the modeling of cell volume changes based on fluorescence intensity changes. We have utilized the method to study the role of a carboxyl-terminus aquaporin (AQP)-2 phosphorylation site, which is known to affect AQP2 membrane trafficking, in heterologous type I Madin-Darby canine kidney cells. We find that water permeability in cells expressing phosphorylation site mutants was in the following order: AQP2-S256D > AQP2 wild-type > AQP2-S256A. We propose that the method can be applied to study AQP function and more generally to study cell volume changes in adherent cell lines. Furthermore, it should be adaptable for AQP inhibitor screening in chemical compound libraries.

The 5A structure of heterologously expressed plant aquaporin SoPIP2;1.

SoPIP2;1 is one of the major integral proteins in spinach leaf plasma membranes. In the Xenopus oocyte expression system its water channel activity is regulated by phosphorylation at the C terminus and in the first cytosolic loop. To assess its structure, SoPIP2;1 was heterologously expressed in Pichia pastoris as a His-tagged protein and in the non-tagged form. Both forms were reconstituted into 2D crystals in the presence of lipids. Tubular crystals and double-layered crystalline sheets of non-tagged SoPIP2;1 were observed and analyzed by cryo-electron microscopy. Crystalline sheets were highly ordered and diffracted electrons to a resolution of 2.96A. High-resolution projection maps of tilted specimens provided a 3D structure at 5A resolution. Superposition of the SoPIP2;1 potential map with the atomic model of AQP1 demonstrates the generally well conserved overall structure of water channels. Differences concerning the extracellular loop A explain the particular crystal contacts between oppositely oriented membrane sheets of SoPIP2;1 2D crystals, and may have a function in rapid volume changes observed in stomatal guard cells or mesophyll protoplasts. This crystal packing arrangement provides access to the phosphorylated C terminus as well as the loop B phosphorylation site for studies of channel gating.

Progress in the analysis of membrane protein structure and function.

Structural information on membrane proteins is sparse, yet they represent an important class of proteins that is encoded by about 30% of all genes. Progress has primarily been achieved with bacterial proteins, but efforts to solve the structure of eukaryotic membrane proteins are also increasing. Most of the structures currently available have been obtained by exploiting the power of X-ray crystallography. Recent results, however, have demonstrated the accuracy of electron crystallography and the imaging power of the atomic force microscope. These instruments allow membrane proteins to be studied while embedded in the bi-layer, and thus in a functional state. The low signal-to-noise ratio of cryo-electron microscopy is overcome by crystallizing membrane proteins in a two-dimensional protein-lipid membrane, allowing its atomic structure to be determined. In contrast, the high signal-to-noise ratio of atomic force microscopy allows individual protein surfaces to be imaged at sub-nanometer resolution, and their conformational states to be sampled. This review summarizes the steps in membrane protein structure determination and illuminates recent progress.

Water permeation across biological membranes: mechanism and dynamics of aquaporin-1 and GlpF.

"Real time" molecular dynamics simulations of water permeation through human aquaporin-1 (AQP1) and the bacterial glycerol facilitator GlpF are presented. We obtained time-resolved, atomic-resolution models of the permeation mechanism across these highly selective membrane channels. Both proteins act as two-stage filters: Conserved fingerprint [asparagine-proline-alanine (NPA)] motifs form a selectivity-determining region; a second (aromatic/arginine) region is proposed to function as a proton filter. Hydrophobic regions near the NPA motifs are rate-limiting water barriers. In AQP1, a fine-tuned water dipole rotation during passage is essential for water selectivity. In GlpF, a glycerol-mediated "induced fit" gating motion is proposed to generate selectivity for glycerol over water.

A refined structure of human aquaporin-1.

A refined structure of the human water channel aquaporin-1 is presented. The model rests on the high resolution X-ray structure of the homologous bacterial glycerol transporter GlpF, electron crystallographic data at 3.8 A resolution and a multiple sequence alignment of the aquaporin superfamily. The crystallographic R and free R values (36.7% and 37.8%) for the refined structure are significantly lower than for previous models. Improved geometry and enhanced stability in molecular dynamics simulations demonstrate a significant improvement of the aquaporin-1 structure. Comparison with previous aquaporin-1 models shows significant differences, not only in the loop regions, but also in the core of the water channel.

Essential dynamics of reversible peptide folding: memory-free conformational dynamics governed by internal hydrogen bonds.

A principal component analysis has been applied on equilibrium simulations of a beta-heptapeptide that shows reversible folding in a methanol solution. The analysis shows that the configurational space contains only three dense sub-states. These states of relatively low free energy correspond to the "native" left-handed helix, a partly helical intermediate, and a hairpin-like structure. The collection of unfolded conformations form a relatively diffuse cloud with little substructure. Internal hydrogen-bonding energies were found to correlate well with the degree of folding. The native helical structure folds from the N terminus; the transition from the major folding intermediate to the native helical structure involves the formation of the two most C-terminal backbone hydrogen bonds. A four-state Markov model was found to describe transition frequencies between the conformational states within error limits, indicating that memory-effects are negligible beyond the nanosecond time-scale. The dominant native state fluctuations were found to be very similar to unfolding motions, suggesting that unfolding pathways can be inferred from fluctuations in the native state. The low-dimensional essential subspace, describing 69% of the collective atomic fluctuations, was found to converge at time-scales of the order of one nanosecond at all temperatures investigated, whereas folding/unfolding takes place at significantly longer time-scales, even above the melting temperature.

The fold of human aquaporin 1.

The fold of human aquaporin 1 is determined from cryo-electron microscopic data at 4.5 A resolution. The monomeric structure consists of two transmembrane triple helices arranged around a pseudo-2-fold axis connected by a long flexible extracellular loop. Each triplet contains between its second and third helix a functional loop containing the highly conserved fingerprint NPA motif. These functional loops are assumed to fold inwards between the two triplets, thereby forming the heart of the water channel. The helix topology was determined from the directionality pattern of each of the six transmembrane helices with respect to the membrane, together with constraints defined by the sequence and atomic force microscopy data. The directionality of the helices was determined by collecting the best-fitting orientations resulting from a search through the three-dimensional experimental map for a large number of alpha-helical fragments. Tests on cryo-electron crystallographic bacteriorhodopsin data suggest that our method is generally applicable to determine the topology of helical proteins for which only medium-resolution electron microscopy data are available.

A kinetic model for the internal motions of proteins: diffusion between multiple harmonic wells.

The dynamics of collective protein motions derived from Molecular Dynamics simulations have been studied for two small model proteins: initiation factor I and the B1 domain of Protein G. First, we compared the structural fluctuations, obtained by local harmonic approximations in different energy minima, with the ones revealed by large scale molecular dynamics (MD) simulations. It was found that a limited set of harmonic wells can be used to approximate the configurational fluctuations of these proteins, although any single harmonic approximation cannot properly describe their dynamics. Subsequently, the kinetics of the main (essential) collective protein motions were characterized. A dual-diffusion behavior was observed in which a fast type of diffusion switches to a much slower type in a typical time of about 1-3 ps. From these results, the large backbone conformational fluctuations of a protein may be considered as "hopping" between multiple harmonic wells on a basically flat free energy surface.

Conformational changes in the chaperonin GroEL: new insights into the allosteric mechanism.

Conformational changes are known to play a crucial role in the function of the bacterial GroE chaperonin system. Here, results are presented from an essential dynamics analysis of known experimental structures and from computer simulations of GroEL using the CONCOORD method. The results indicate a possible direct form of inter-ring communication associated with internal fluctuations in the nucleotide-binding domains upon nucleotide and GroES binding that are involved in the allosteric mechanism of GroEL. At the level of conformational transitions in entire GroEL rings, nucleotide-induced structural changes were found to be distinct and in principle uncoupled from changes occurring upon GroES binding. However, a coupling is found between nucleotide-induced conformational changes and GroES-mediated transitions, but only in simulations of GroEL double rings, and not in simulations of single rings. This provides another explanation for the fact that GroEL functions a double ring system.

Domain motions in bacteriophage T4 lysozyme: a comparison between molecular dynamics and crystallographic data.

A comparison of a series of extended molecular dynamics (MD) simulations of bacteriophage T4 lysozyme in solvent with X-ray data is presented. Essential dynamics analyses were used to derive collective fluctuations from both the simulated trajectories and a distribution of crystallographic conformations. In both cases the main collective fluctuations describe domain motions. The protein consists of an N- and C-terminal domain connected by a long helix. The analysis of the distribution of crystallographic conformations reveals that the N-terminal helix rotates together with either of these two domains. The main domain fluctuation describes a closure mode of the two domains in which the N-terminal helix rotates concertedly with the C-terminal domain, while the domain fluctuation with second largest amplitude corresponds to a twisting mode of the two domains, with the N-terminal helix rotating concertedly with the N-terminal domain. For the closure mode, the difference in hinge-bending angle between the most open and most closed X-ray structure along this mode is 49 degrees. In the MD simulation that shows the largest fluctuation along this mode, a rotation of 45 degrees was observed. Although the twisting mode has much less freedom than the closure mode in the distribution of crystallographic conformations, experimental results suggest that it might be functionally important. Interestingly, the twisting mode is sampled more extensively in all MD simulations than it is in the distribution of X-ray conformations.

Structural homology of the central conserved region of the attachment protein G of respiratory syncytial virus with the fourth subdomain of 55-kDa tumor necrosis factor receptor.

The attachment protein G of respiratory syncytial virus (RSV) has a modular architecture. The ectodomain of the protein comprises a small folded conserved region which is bounded by two mucin-like regions. In this study, a sequence and structural homology is described between this central conserved region of RSV-G and the fourth subdomain of the 55-kDa tumor necrosis factor receptor (TNFr). The three-dimensional structures of RSV-G and human TNFr were previously determined with NMR spectroscopy and X-ray crystallography, respectively. The C-terminal part of both subdomains fold into a cystine noose connected by two cystine bridges with the same spacing between cysteine residues and the same topology. Although a general structural similarity is observed, there are differences in secondary structure and other structural features. Molecular Dynamics calculations show that the BRSV-G NMR structure of the cystine noose is stable and that the TNFr crystal structure of the cystine noose drifts towards the BRSV-G NMR structure in the simulated solution environment. By homology modelling a model was built for the unresolved N-terminal part of the central conserved region of RSV-G. The functions for both protein domains are not known but the structural similarity of both protein domains suggests a similar function. Although the homology suggests that the cystine noose of RSV-G may interfere with the antiviral and apoptotic effect of TNF, the biological activity remains to be proven.

Protein dynamics derived from clusters of crystal structures.

A method is presented to mathematically extract concerted structural transitions in proteins from collections of crystal structures. The "essential dynamics" procedure is used to filter out small-amplitude fluctuations from such a set of structures; the remaining large conformational changes describe motions such as those important for the uptake/release of substrate/ligand and in catalytic reactions. The method is applied to sets of x-ray structures for a number of proteins, and the results are compared with the results from essential dynamics as applied to molecular dynamics simulations of those proteins. A significant degree of similarity is found, thereby providing a direct experimental basis for the application of such simulations to the description of large concerted motions in proteins.

Prediction of protein conformational freedom from distance constraints.

A method is presented that generates random protein structures that fulfil a set of upper and lower interatomic distance limits. These limits depend on distances measured in experimental structures and the strength of the interatomic interaction. Structural differences between generated structures are similar to those obtained from experiment and from MD simulation. Although detailed aspects of dynamical mechanisms are not covered and the extent of variations are only estimated in a relative sense, applications to an IgG-binding domain, an SH3 binding domain, HPr, calmodulin, and lysozyme are presented which illustrate the use of the method as a fast and simple way to predict structural variability in proteins. The method may be used to support the design of mutants, when structural fluctuations for a large number of mutants are to be screened. The results suggest that motional freedom in proteins is ruled largely by a set of simple geometric constraints.

An extended sampling of the configurational space of HPr from E. coli.

Recently, we developed a method (Amadei et al., J. Biomol. Str. Dyn. 13: 615-626; de Groot et al., J. Biomol. Str. Dyn. 13: 741-751, 1996) to obtain an extended sampling of the configurational space of proteins, using an adapted form of molecular dynamics (MD) simulations, based on the essential dynamics (ED) (Amadei et al., Proteins 17:412-425, 1993) method. In the present study, this ED sampling technique is applied to the histidine-containing phosphocarrier protein HPr from Escherichia coli. We find a cluster of conformations that is an order of magnitude larger than that found for a usual MD simulation of comparable length. The structures in this cluster are geometrically and energetically comparable to NMR structures. Moreover, on average, this large cluster satisfies nearly all NMR-derived distance restraints.

The consistency of large concerted motions in proteins in molecular dynamics simulations.

A detailed investigation is presented into the effect of limited sampling time and small changes in the force field on molecular dynamics simulations of a protein. Thirteen independent simulations of the B1 IgG-binding domain of streptococcal protein G were performed, with small changes in the simulation parameters in each simulation. Parameters studied included temperature, bond constraints, cut-off radius for electrostatic interactions, and initial placement of hydrogen atoms. The essential dynamics technique was used to reveal dynamic differences between the simulations. Similar essential dynamics properties were found for all simulations, indicating that the large concerted motions found in the simulations are not particularly sensitive to small changes in the force field. A thorough investigation into the stability of the essential dynamics properties as derived from a molecular dynamics simulation of a few hundred picoseconds is provided. Although the definition of the essential modes of motion has not fully converged in these short simulations, the subspace in which these modes are confined is found to be reproducible.

Conformational analysis of retinoids and restriction of their dynamics by retinoid-binding proteins.

An exhaustive sampling of the configurational space of all-trans retinol using a 0.1 microsecond molecular-dynamics simulation is presented. The essential dynamics technique is used to describe the conformational changes in retinol using only three degrees of freedom. The different conformational states of retinol are analysed, and differences in free energy are calculated. The essential dynamics description allows a detailed comparison of free retinol and retinoids bound to retinoid-binding proteins and opens new possibilities in the small-molecule docking field. The dynamics of retinoids when complexed with their binding proteins are restricted, and they are forced into strained conformations. A "spring' model for retinoid binding is proposed. This model is extended to a hypothesis for retinoid binding to visual pigments and bacteriorhodopsin.

Bending of the calmodulin central helix: a theoretical study.

The crystal structure of calcium-calmodulin (CaM) reveals a protein with a typical dumbbell structure. Various spectroscopic studies have suggested that the central linker region of CaM, which is alpha-helical in the crystal structure, is flexible in solution. In particular, NMR studies have indicated the presence of a flexible backbone between residues Lys 77 and Asp 80. This flexibility is related directly to the function of the protein because it enables the N- and C-terminal domains of the protein to move toward each other and bind to the CaM-binding domain of a target protein. We have investigated the flexibility of the CaM central helix by a variety of computational techniques: molecular dynamics (MD) simulations, normal mode analysis (NMA), and essential dynamics (ED) analysis. Our MD results reproduce the experimentally determined location of the bend in a simulation of only the CaM central helix, indicating that the bending point is an intrinsic property of the alpha-helix, for which the remainder of the protein is not important. Interestingly, the modes found by the ED analysis of the MD trajectory are very similar to the lowest frequency modes from the NM analysis and to modes found by an ED analysis of different structures in a set of NMR structures. Electrostatic interactions involving residues Arg 74 and Asp 80 seem to be important for these bending motions and unfolding, which is in line with pH-dependent NMR and CD studies.

Toward an exhaustive sampling of the configurational spaces of the two forms of the peptide hormone guanylin.

The recently introduced Essential Dynamics sampling method is extended such that an exhaustive sampling of the available (backbone) configurational space can be achieved. From an initial Molecular Dynamics simulation an approximated definition of the essential subspace is obtained. This subspace is used to direct subsequent simulations by means of constraint forces. The method is applied to the peptide hormone guanylin, solvated in water, of which the structure was determined recently. The peptide exists in two forms and for both forms, an extensive sampling was produced. The sampling algorithm fills the available space (of the essential coordinates used in the procedure) at a rate that is approximately six to seven times larger than that for traditional Molecular Dynamics. The procedure does not cause any significant perturbation, which is indicated by the fact that free Molecular Dynamics simulations started at several places in the space defined by the Essential Dynamics sample that complete space. Moreover, analyses of the average free Molecular Dynamics step have shown that nowhere except close to the edge of the available space, there are regions where the system shows a drift in a particular direction. This result also shows that in principle, the essential subspace is a constant free energy surface, with well-defined and steep borders, in which the system moves diffusively. In addition, a comparison between two independent essential dynamics sampling runs, of one form of the peptide, shows that the obtained essential subspaces are virtually identical.