The VIZIER project: preparedness against pathogenic RNA viruses

Life-threatening RNA viruses emerge regularly, and often in an unpredictable manner. Yet, the very few drugs available against known RNA viruses have sometimes required decades of research for development. Can we generate preparedness for outbreaks of the, as yet, unknown viruses? The VIZIER (VIral enZymes InvolvEd in Replication) (http://www.vizier-europe.org/) project has been set-up to develop the scientific foundations for countering this challenge to society. VIZIER studies the most conserved viral enzymes (that of the replication machinery, or replicases) that constitute attractive targets for drug-design. The aim of VIZIER is to determine as many replicase crystal structures as possible from a carefully selected list of viruses in order to comprehensively cover the diversity of the RNA virus universe, and generate critical knowledge that could be efficiently utilized to jump-start research on any emerging RNA virus. VIZIER is a multidisciplinary project involving (i) bioinformatics to define functional domains, (ii) viral genomics to increase the number of characterized viral genomes and prepare defined targets, (iii) proteomics to express, purify, and characterize targets, (iv) structural biology to solve their crystal structures, and (v) pre-lead discovery to propose active scaffolds of antiviral molecules.

Computational inhibitor design against malaria plasmepsins

Plasmepsins are aspartic proteases involved in the degradation of the host cell hemoglobin that is used as a food source by the malaria parasite. Plasmepsins are highly promising as drug targets, especially when combined with the inhibition of falcipains that are also involved in hemoglobin catabolism. In this review, we discuss the mechanism of plasmepsins I-IV in view of the interest in transition state mimetics as potential compounds for lead development. Inhibitor development against plasmepsin II as well as relevant crystal structures are summarized in order to give an overview of the field. Application of computational techniques, especially binding affinity prediction by the linear interaction energy method, in the development of malarial plasmepsin inhibitors has been highly successful and is discussed in detail. Homology modeling and molecular docking have been useful in the current inhibitor design project, and the combination of such methods with binding free energy calculations is analyzed.

Ions and blockers in potassium channels: insights from free energy simulations

Potassium ion channels enable efficient and selective permeation of K+ ions across nonpolar biological membranes. Here we review the results of recent free energy calculations related to the permeation of monovalent cations through K+ channels and to the channel inhibition by blocker compounds. In particular, the progress in computational studies of the bacterial KcsA channel is discussed.

The linear interaction energy method for predicting ligand binding free energies

An overview of the simplified linear interaction energy (LIE) method for calculation of ligand binding free energies is given. This method is based on force field estimations of the receptor-ligand interactions and thermal conformational sampling. A notable feature is that the binding energetics can be predicted by considering only the intermolecular interactions between the ligand and receptor. The approximations behind this approach are examined and different parametrizations of the model are discussed. In general, LIE type of methods appear particularly useful for computational drug lead optimization.

K(+)/Na(+) selectivity of the KcsA potassium channel from microscopic free energy perturbation calculations

Microscopic molecular dynamics free energy perturbation calculations of the K(+)/Na(+) selectivity in the KcsA potassium channel, based on its experimental three-dimensional structure, are reported. The relative binding free energies for K(+) and Na(+) in the most relevant ion occupancy states of the four-site selectivity filter are calculated. The previously proposed mechanism for ion permeation through the KcsA channel is predicted, in agreement with available experimental data, to have a significant selectivity for K(+) over Na(+). The calculations also show that the individual 'binding site' selectivities are generally not additive and the doubly loaded states of the filter thus display cooperative effects. The only site that is not K(+) selective is that which is located at the entrance to the internal water cavity, suggesting the possibility that internal Na(+) could block outward currents.

Computational analysis of binding of P1 variants to trypsin

The binding of P1 variants of bovine pancreatic trypsin inhibitor (BPTI) to trypsin has been investigated by means of molecular dynamics simulations. The specific interaction formed between the amino acid at the primary binding (P1) position of the binding loop of BPTI and the specificity pocket of trypsin was estimated by use of the linear interaction energy (LIE) method. Calculations for 13 of the naturally occurring amino acids at the P1 position were carried out, and the results obtained were found to correlate well with the experimental binding free energies. The LIE calculations rank the majority of the 13 variants correctly according to the experimental association energies and the mean error between calculated and experimental binding free energies is only 0.38 kcal/mole, excluding the Glu and Asp variants, which are associated with some uncertainties regarding protonation and the possible presence of counter-ions. The three-dimensional structures of the complex with three of the P1 variants (Asn, Tyr, and Ser) included in this study have not at present been solved by any experimental techniques and, therefore, were modeled on the basis of experimental data from P1 variants of similar size. Average structures were calculated from the MD simulations, from which specific interactions explaining the broad variation in association energies were identified. The present study also shows that explicit treatment of the complex water-mediated hydrogen bonding network at the protein-protein interface is of crucial importance for obtaining reliable binding free energies. The successful reproduction of relative binding energies shows that this type of methodology can be very useful as an aid in rational design and redesign of biologically active macromolecules.

Sensitivity of an empirical affinity scoring function to changes in receptor-ligand complex conformations

A combination of empirical scoring and conformational sampling for ligand binding affinity prediction is examined. The behaviour of a scoring function with respect to the sensitivity to conformational changes is investigated using ensembles of structures generated by molecular dynamics simulation. The correlation between the calculated score and the coordinate deviation from the experimental structure is clear for the complex of arabinose with arabinose-binding protein, which is dominated by hydrogen bond interactions, while the score calculated for the hydrophobic complex between retinol and retinol binding protein is rather insensitive to ligand conformational changes. For typical ensembles of structures generated by molecular dynamics at 300 K, the variation of the calculated score is considerably smaller than that of the underlying molecular mechanics interaction energies.

Design, synthesis, computational prediction, and biological evaluation of ester soft drugs as inhibitors of dihydrofolate reductase from Pneumocystis carinii

A series of lipophilic soft drugs structurally related to the nonclassical dihydrofolate reductase (DHFR) inhibitors trimetrexate and piritrexim have been designed, synthesized, and evaluated in DHFR assays, with special emphasis on the inhibition of P. carinii DHFR. The best inhibitors, encompassing an ester bond in the bridge connecting the two aromatic systems, were approximately 10 times less potent than trimetrexate and piritrexim. The metabolites were designed to be poor inhibitors. Furthermore, molecular dynamics simulations of three ligands in complex with DHFR from Pneumocystis carinii and from the human enzyme were conducted in order to better understand the factors determining the selectivity. A correct ranking of the relative inhibition of DHFR was achieved utilizing the linear interaction energy method. The soft drugs are intended for local administration. One representative ester was selected for a pharmacokinetic study in rats where it was found to undergo fast metabolic degradation to the predicted inactive metabolites.

Electrostatic effects play a central role in cold adaptation of trypsin

Organisms that live in constantly cold environments have to adapt their metabolism to low temperatures, but mechanisms of enzymatic adaptation to cold environments are not fully understood. Cold active trypsin catalyses reactions more efficiently and binds ligands more strongly in comparison to warm active trypsin. We have addressed this issue by means of comparative free energy calculations studying the binding of positively charged ligands to two trypsin homologues. Stronger inhibition of the cold active trypsin by benzamidine and positively charged P1-variants of BPTI is caused by rather subtle electrostatic effects. The different affinity of benzamidine originates solely from long range interactions, while the increased binding of P1-Lys and -Arg variants of BPTI is attributed to both long and short range effects that are enhanced in the cold active trypsin compared to the warm active counterpart. Electrostatic interactions thus provide an efficient strategy for cold adaptation of trypsin.

The catalytic mechanism of protein tyrosine phosphatases revisited

Experimental and theoretical studies of the catalytic mechanism in protein tyrosine phosphatases and dual specific phosphatases are reviewed. The structural properties of these enzymes contributing to the efficient rate enhancement of phosphate monoester hydrolysis have been established during the last decade. There are, however, uncertainties in the interpretation of available experimental data that make the commonly assumed reaction mechanism somewhat doubtful. Theoretical calculations as well as analysis of crystal structures point towards an alternative interpretation of the ionisation state in the reactive complex.

Mechanisms of tetraethylammonium ion block in the KcsA potassium channel

We report results from automated docking and microscopic molecular dynamics simulations of the tetraethylammonium (TEA) complexes with KcsA. Binding modes and energies for TEA binding at the external and internal sides of the channel pore are examined utilising the linear interaction energy method. Effects of the channel ion occupancy (based on our previous results for the ion permeation mechanisms) on the binding energies are considered. Calculations show that TEA forms stable complexes at both the external and internal entrances of the selectivity filter. Furthermore, the effects of the Y82V mutation are evaluated and the results show, in agreement with experimental data, that the mutant has a significantly reduced binding affinity for TEA at the external binding site, which is attributed to stabilising hydrophobic interactions between the ligand and the tyrosines.

Computational modelling of inhibitor binding to human thrombin

Thrombin is an essential protein involved in blood clot formation and an important clinical target, since disturbances of the coagulation process cause serious cardiovascular diseases such as thrombosis. Here we evaluate the performance of a molecular dynamics based method for predicting the binding affinities of different types of human thrombin inhibitors. For a series of eight ligands the method ranks their relative affinities reasonably well. The binding free energy difference between high and low affinity representatives in the test set is quantitatively reproduced, as well as the stereospecificity for a chiral inhibitor. The original parametrisation of this linear interaction energy method requires the addition of a constant energy term in the case of thrombin. This yields a mean unsigned error of 0.68 kcal/mol for the absolute binding free energies. This type of approach is also useful for elucidating three-dimensional structure-activity relationships in terms of microscopic interactions of the ligands with the solvated enzyme.

Computational predictions of binding affinities to dihydrofolate reductase: synthesis and biological evaluation of methotrexate analogues

The relative binding affinities to human dihydrofolate reductase of four new potential antifolates, containing ester linkages between the two aromatic systems, were estimated by free energy perturbation simulations. The ester analogue, predicted to exhibit the highest binding affinity to human dihydrofolate reductase, and a reference ester (more structurally related to methotrexate) were synthesized. As deduced from the measured IC(50) values, the calculated ranking of the ligands was correct although a greater difference in affinity was indicated by the experimental measurements. Among the new antifolates the most potent inhibitor exhibited a similar pharmacokinetic profile to methotrexate but lacked activity in a complex antiarthritic model in rat in vivo.

A computational study of ion binding and protonation states in the KcsA potassium channel

We report results from microscopic molecular dynamics and free energy perturbation simulations of the KcsA potassium channel based on its experimental atomic structure. Conformational properties of selected amino acid residues as well as equilibrium positions of K(+) ions inside the selectivity filter and the internal water cavity are examined. Positions three and four (counting from the extracellular site) in the experimental structure correspond to distinctly separate binding sites for K(+) ions inside the selectivity filter. The protonation states of Glu71 and Asp80, which are close to each other and to the selectivity filter, as well as K(+) binding energies are determined using free energy perturbation calculations. The Glu71 residue which is buried inside a protein cavity is found to be most stable in the neutral form while the solvent exposed Asp80 is ionized. The channel altogether exothermically binds up to three ions, where two of them are located inside the selectivity filter and one in the internal water cavity. Ion permeation mechanisms are discussed in relation to these results.

Computer simulation of primary kinetic isotope effects in the proposed rate-limiting step of the glyoxalase I catalyzed reaction

The proposed rate-limiting step of the glyoxalase I catalyzed reaction is the proton abstraction from the C1 carbon of the substrate by Glu(172). Here we examine primary kinetic isotope effects and the influence of quantum dynamics on this process by computer simulations. The calculations utilize the empirical valence bond method in combination with the molecular dynamics free energy perturbation technique and path integral simulations. For the enzyme-catalyzed reaction a H/D kinetic isotope effect of 5.0 +/- 1. 3 is predicted in reasonable agreement with the experimental result of about 3. Furthermore, the magnitude of quantum mechanical effects is found to be very similar for the enzyme reaction and the corresponding uncatalyzed process in solution, in agreement with other studies. The problems associated with attaining the required accuracy in order for the present approach to be useful as a diagnostic tool for the study of enzyme reactions are also discussed.

Ion permeation mechanism of the potassium channel

Ion-selective channels enable the specific permeation of ions through cell membranes and provide the basis of several important biological functions; for example, electric signalling in the nervous system. Although a large amount of electrophysiological data is available, the molecular mechanisms by which these channels can mediate ion transport remain a significant unsolved problem. With the recently determined crystal structure of the representative K+ channel (KcsA) from Streptomyces lividans, it becomes possible to examine ion conduction pathways on a microscopic level. K+ channels utilize multi-ion conduction mechanisms, and the three-dimensional structure also shows several ions present in the channel. Here we report results from molecular dynamics free energy perturbation calculations that both establish the nature of the multiple ion conduction mechanism and yield the correct ion selectivity of the channel. By evaluating the energetics of all relevant occupancy states of the selectivity filter, we find that the favoured conduction pathway involves transitions only between two main states with a free difference of about 5 kcal mol(-1). Other putative permeation pathways can be excluded because they would involve states that are too high in energy.

Prediction of a ligand-induced conformational change in the catalytic core of Cdc25A

The cell cycle control phosphatases Cdc25 are dual specificity phosphatases that dephosphorylate both phosphothreonine and phosphotyrosine residues on their substrate proteins. The determination of the apo-protein structure of Cdc25A revealed that this enzyme has a completely different fold compared to all other phosphatases crystallised to date. The conformation of the active site residues does not seem very suitable for catalysis in this unliganded structure. We have studied some structural features of the Cdc25A apo-structure and a modelled Cdc25A-ligand complex by molecular dynamics simulations. The simulations predict a conformational change in the peptide backbone of the complex, which is not observed in the apo-structure. This ligand-induced conformational change yields a structure that is similar to other protein tyrosine phosphatase-ligand complexes that have been crystallised. The change in conformation takes place in the position between a serine and a glutamic acid residue in the phosphate binding loop. We suggest that this type of conformational change is an important molecular switch in the catalytic process.

Long-range electrostatic effects on peptide folding

The reversible folding/unfolding of a short peptide in solution is studied by molecular dynamics simulations. The effects of long-range electrostatic interactions are examined and found to be important both for the equilibrium between folded and unfolded states and the dynamics of the folding process. The neglect of long-range electrostatics leads to an increased population of unfolded states and increased structural fluctuations. When such interactions are taken into account, the peptide unfolds and folds to the experimentally determined structure several times during a 25 ns trajectory, with approximately equal populations of folded and unfolded states in the neighborhood of its proposed melting temperature. The effect of using spherical boundary conditions rather than periodic ones does not appear to have any major effect on the folding dynamics.

Mechanism of substrate dephosphorylation in low Mr protein tyrosine phosphatase

Substrate dephosphorylation by the low molecular weight protein tyrosine phosphatases proceeds via nucleophilic substitution at the phosphorous atom yielding a cysteinyl phosphate intermediate. However, several questions regarding the exact reaction mechanism remain unanswered. Starting from the crystal structure of the enzyme we study the energetics of this reaction, using the empirical valence bond method in combination with molecular dynamics and free energy perturbation simulations. The free energy profiles of two mechanisms corresponding to different protonation states of the reacting groups are examined along stepwise and concerted pathways. The activation barriers calculated relative to the enzyme-substrate complex are very similar for both monoanionic and dianionic substrates, but taking the substrate binding step into account shows that hydrolysis of monoanionic substrates is strongly favored by the enzyme, because a dianionic substrate will not bind when the reacting cysteine is ionized. The calculated activation barrier for dephosphorylation of monoanionic phenyl phosphate according to this novel mechanism is 14 kcal mol(-1), which is in good agreement with experimental data. Proteins 1999;36:370-379.

Computational modeling of the rate limiting step in low molecular weight protein tyrosine phosphatase

Hydrolysis of the phosphoenzyme intermediate is the second and rate limiting step of the reaction catalyzed by the protein tyrosine phosphatases (PTPs). The cysteinyl phosphate thioester bond is cleaved by nucleophilic displacement where an active site water molecule attacks the phosphorus atom. Starting from the crystal structure of the low molecular weight PTP, we study the energetics of this reaction utilizing the empirical valence bond method in combination with molecular dynamics and free energy perturbation simulations. The reactions of the wild-type as well as the D129A and C17S mutants are modeled. For the D129A mutant, which lacks the general acid/base residue Asp-129, an alternative reaction mechanism is proposed. The calculated activation barriers are in all cases in good agreement with experimental reaction rates. The present results together with earlier computational and experimental work now provide a detailed picture of the complete reaction mechanism in many PTPs. The key role played by the structurally invariant signature motif in stabilizing a double negative charge is reflected by its control of the energetics of both transition states and the reaction intermediate.

Energetics of the proposed rate-determining step of the glyoxalase I reaction

The proposed rate-limiting step of the reaction catalyzed by glyoxalase I is the proton abstraction from the C1 carbon atom of the substrate by a glutamate residue, resulting in a high-energy enolate intermediate. This proton transfer reaction was modelled using molecular dynamics and free energy perturbation simulations, with the empirical valence bond method describing the potential energy surface of the system. The calculated rate constant for the reaction is approximately 300-1500 s(-1) with Zn2+, Mg2+ or Ca2+ bound to the active site, which agrees well with observed kinetics of the enzyme. Furthermore, the results imply that the origin of the catalytic rate enhancement is mainly associated with enolate stabilization by the metal ion.

Mechanistic alternatives in phosphate monoester hydrolysis: what conclusions can be drawn from available experimental data?

Phosphate monoester hydrolysis reactions in enzymes and solution are often discussed in terms of whether the reaction pathway is associative or dissociative. Although experimental results for solution reactions have usually been considered as evidence for the second alternative, a closer thermodynamic analysis of observed linear free energy relationships shows that experimental information is consistent with the associative, concerted and dissociative alternatives.

Q: a molecular dynamics program for free energy calculations and empirical valence bond simulations in biomolecular systems

A new molecular dynamics program for free energy calculations in biomolecular systems is presented. It is principally designed for free energy perturbation simulations, empirical valence bond calculations, and binding affinity estimation by linear interaction energy methods. Evaluation of ligand-binding selectivity and free energy profiles for nucleophile activation in two protein tyrosine phosphatases as well as absolute binding affinity estimation for a lysine-binding protein are given as examples.

Computation of affinity and selectivity: binding of 2,4-diaminopteridine and 2,4-diaminoquinazoline inhibitors to dihydrofolate reductases

Binding energy calculations for complexes of mutant and wild-type human dihydrofolate reductases with 2,4-diaminopteridine and 2,4-diaminoquinazoline inhibitors are reported. Quantitative insight into binding energetics of these molecules is obtained from calculations based on force field energy evaluation and thermal sampling by molecular dynamics simulations. The calculated affinity of methotrexate for wild-type and mutant enzymes is reasonably well reproduced. Truncation of the methotrexate glutamate tail results in a loss of affinity by several orders of magnitude. No major difference in binding strength is predicted between the pteridines and the quinazolìnes, while the N-methyl group present in methotrexate appears to confer significantly stronger binding. The recent improvement, which is used here, of our linear interaction energy method for binding affinity prediction, as well as problems with treating charged and flexible ligands are discussed. This approach should be suitable in a drug discovery context for prediction of binding energies of new inhibitors prior to their synthesis, when some information about the binding mode is available.

Ligand binding affinity prediction by linear interaction energy methods

A recent method for estimating ligand binding affinities is extended. This method employs averages of interaction potential energy terms from molecular dynamics simulations or other thermal conformational sampling techniques. Incorporation of systematic deviations from electrostatic linear response, derived from free energy perturbation studies, into the absolute binding free energy expression significantly enhances the accuracy of the approach. This type of method may be useful for computational prediction of ligand binding strengths, e.g., in drug design applications.

Cyclic HIV-1 protease inhibitors derived from mannitol: synthesis, inhibitory potencies, and computational predictions of binding affinities

Ten C2-symmetric cyclic urea and sulfamide derivatives have been synthesized from L-mannonic gamma-lactone and D-mannitol. The results of experimental measurement of their inhibitory potencies against HIV-1 protease were compared to calculated free energies of binding derived from molecular dynamics (MD) simulations. The compounds were selected, firstly, to enable elucidation of the role of stereochemistry for binding affinity (1a-d) and, secondly, to allow evaluation of the effects of variation in the link to the P1 and P1' phenyl groups on affinity (1a and 2-5). Thirdly, compounds with hydrogen bond-accepting or-donating groups attached to the phenyl groups in the P2 and P2' side chains (6 and 7) were selected. Binding free energies were estimated by a linear response method, whose predictive power for estimating binding affinities from MD simulations was demonstrated.

Energetics of nucleophile activation in a protein tyrosine phosphatase

The nucleophilic attack by cysteine 12 in the low-molecular-weight protein tyrosine phosphatase is believed to be carried out by the thiolate anion form of this residue. We here study the energetics of proton transfer between the thiol group of cysteine 12 and a substrate phosphate oxygen atom, to examine the effects of the enzymic environment on the stability of the thiolate nucleophile. This is done by molecular dynamics and free energy perturbation simulations, utilizing the empirical valence bond method to describe the potential surface of the system. The calculations show that the protein environment significantly stabilizes the thiolate ion, thereby setting the stage for the nucleophilic attack. We compare these results with those from further simulations of a mutant enzyme, and demonstrate the importance of serine 19 in thiolate stabilization.

Computer simulation of the triosephosphate isomerase catalyzed reaction

A major challenge for theoretical simulation methods is the calculation of enzymic reaction rates directly from the three-dimensional protein structure together with some idea of the chemical reaction mechanism. Here, we report the evaluation of a complete free energy profile for all the elementary steps of the triosephosphate isomerase catalyzed reaction using such an approach. The results are compatible with available experimental data and also suggest which of the possible reaction intermediates is kinetically observable. In addition to previously identified catalytic residues, the simulations show that a crystallographically observed active site water molecule plays an important role during catalysis and an intersubunit interaction that could explain the low activity of the monomeric enzyme is also observed. The calculations clearly demonstrate the important catalytic effects associated with stabilization of charged high energy intermediates and reduction of reorganization energy, which are likely to be general principles of enzyme catalyzed charge transfer and separation reactions.

Estimation of binding free energies for HIV proteinase inhibitors by molecular dynamics simulations

Absolute binding free energies for three inhibitors of HIV-1 proteinase were estimated from molecular dynamics simulations by a recently reported linear approximation procedure. The results were in fairly good agreement with experimental binding data. Two of the inhibitors were very similar and, for comparison, their relative free energies of binding were also calculated by free energy perturbation methods, giving virtually the same result. Effects of cut-off radii and charge states of the protein model were examined. The effects of pH on binding of one of the inhibitors were predicted.

Sugar recognition by a glucose/galactose receptor. Evaluation of binding energetics from molecular dynamics simulations

A new theoretical method for free energy calculations is used to compute the absolute binding constants for beta-D-glucose and methyl-beta-D-galactoside to the periplasmic glucose/galactose receptor from Salmonella typhimurium. The computer simulation results agree well with available experimental data and make it possible to assess the sources of both the high affinity as well as the specificity for glucose. It was found that the major contribution to the binding energy comes from electrostatic interactions and particularly hydrogen bonds of the charge-dipole type. We also predict the structure of the complex with methyl-galactoside as this has not yet been experimentally determined.

A new method for predicting binding affinity in computer-aided drug design

A new semi-empirical method for calculating free energies of binding from molecular dynamics (MD) simulations is presented. It is based on standard thermodynamic cycles and on a linear approximation of polar and non-polar free energy contributions from the corresponding MD averages. The method is tested on a set of endothiapepsin inhibitors and found to give accurate results both for absolute as well as relative free energies.

A molecular model for the retinol binding protein-transthyretin complex

A three-dimensional model for the complex between human serum retinol binding protein and transthyretin (formerly named prealbumin) is presented. The model was obtained by interactive rigid-body computer graphics docking and the characterization of the molecular surfaces in terms of fractal dimension. Available experimental data, as well as results from molecular dynamics calculations, support the proposed model.

Computer simulation of the initial proton transfer step in human carbonic anhydrase I

The initial water proteolysis step in the proton transfer "half-reaction" of human carbonic anhydrase I is simulated using the empirical valence bond method in combination with free energy perturbation molecular dynamics calculations. A free energy profile for the enzyme catalysed reaction and the corresponding pKa associated with ionization of the zinc-bound water is calculated. The obtained pKa value of 7 to 8 appears to be in good agreement with experimental observations and the calculated rate constant for this step is also compatible with kinetic data. The simulations clearly emphasize the important electrostatic effect associated with the catalytic zinc ion.

Computer simulations of enzymatic reactions: examination of linear free-energy relationships and quantum-mechanical corrections in the initial proton-transfer step of carbonic anhydrase

Computer simulation approaches can provide a powerful tool for correlating the structure of enzymes with their catalytic activity. One of the most effective ways of simulating enzymatic reactions is provided by the empirical valence bond method. The general applicability of this method has been demonstrated in several enzymatic reactions and it is reexamined here in a study of the initial proton-transfer step in the catalytic reaction of carbonic anhydrase. The simulations produce a rate constant which is in agreement with the observed kinetic data and emphasizes the importance of the electrostatic effect associated with the catalytic zinc ion. The calculations are also used to examine the validity of linear free-energy relationships (LFERs) in enzyme catalysis and to evaluate quantum-mechanical corrections of the calculated rate constant. It is found that LFERs are valid in the present case and it is argued that this reflects the fact that the protein responds linearly to the development of electrostatic forces during the reaction. It is concluded that the present approach can be used to augment experimental studies in establishing the general validity of LFERs. It is noted, however, that such relationships are much more valid for transitions between different resonance structures than for transitions between reactants and product states.

Enzymes work by solvation substitution rather than by desolvation

Considerable attention has recently been drawn to the hypothesis that enzymes catalyze their reactions by displacing solvent and creating an environment similar to the gas phase for the reacting substrates. This "desolvation hypothesis" is reexamined in this paper by defining a common reference energy for reactions in various environments. It is argued that consistent attempts to describe the actual energetics of enzymatic reactions, taking either gas phase or solution as a reference, would contradict the above hypothesis. That is, the enzyme does remove water molecules from its substrate, but substitutes these molecules for another polar environment (namely, its active site). By taking amide hydrolysis as an example, we use experimentally estimated solvation energies and analyze the reaction profile in the gas phase, in solution, and in enzyme active sites. We show that the gas-phase reaction is characterized by an enormous activation barrier (associated with forming the charged nucleophile from neutral fragments), although the nucleophilic attack is essentially barrierless. On the other hand, the enzyme and solution reactions are found to have similar reaction profiles, with a lower activation barrier for the enzymatic reaction. Presumably, the fact that previous analyses of this problem did not involve the construction of the relevant thermodynamic cycles (and quantitative estimates of the corresponding solvation energies) led to the desolvation hypothesis. Our conclusion is that enzyme active sites provide specific polar environments that do not resemble the gas phase but that are designed for electrostatic stabilization of ionic transition states and that "solvate" these states more than water does.

Energetics of ion permeation through membrane channels. Solvation of Na+ by gramicidin A

Calculations of the solvation energetics for a Na+ ion inside the Gramicidin A channel and in water are presented. The protein dipoles Langevin dipoles (PDLD) method is used to obtain an electrostatic free energy profile for ion permeation through the channel. To gauge the quality of the PDLD results the solvation free energy of a Na+ ion in water and in the center of the channel is also calculated using free energy perturbation (FEP) simulations. The effect of the polarisability of the surrounding lipid membrane is taken into account by representing the membrane by a large grid of polarisable point dipoles. The two methods give similar solvation energies in the interior of the channel and these are less than 5 kcal/mol above the solvation free energy for Na+ in water, in good agreement with experimental data on the activation barriers for ion permeation. It appears that the problems associated with previous calculations of energy profiles in membrane channels can be overcome by a consistent treatment of all the relevant electrostatic contributions. In particular, we find that the induced dipoles of the membrane and the protein contributes with approximately 10 kcal/mol to the solvation energy inside the channel and can therefore not be discarded in a realistic description of ion solvation in the Gramicidin channel.

Calculations of free energy profiles for the staphylococcal nuclease catalyzed reaction

Calculations of the free energy profile for the first two (rate-limiting) steps of the staphylococcal nuclease catalyzed reaction are reported. The calculations are based on the empirical valence bond method in combination with free energy perturbation molecular dynamics simulations. The calculated activation free energy is in good agreement with experimental kinetic data, and the catalytic effect of the enzyme is reproduced without any arbitrary adjustment of parameters. The enormous reduction of the activation barrier (relative to the reference reaction in water) appears to be largely associated with the strong electrostatic effect of the Ca2+ ion and the two arginine residues in the active site. This favorable electrostatic environment reduces the cost of the general-base catalysis step by almost 15 kcal/mol (by stabilizing the OH- nucleophile) and then stabilizes the developing negative charge on the 5'-phosphate group in the second step of the reaction by about 19 kcal/mol. The basic features of the originally postulated enzyme mechanism (Cotton et al., 1979) are found to be compatible with the observed activation free energy. However, the proposed modification of the mechanism (Sepersu et al., 1987), in which Arg 87 interacts only with the pentacoordinated transition state, is supported by the simulations. Further calculations on the D21E mutant also give results in good agreement with kinetic data.

A molecular dynamics study of the C-terminal fragment of the L7/L12 ribosomal protein. II. Effects of intermolecular interactions on structure and dynamics

The crystallographic dimer of the C-terminal fragment (CTF) of the L7/L12 ribosomal protein has been subjected to molecular dynamics (MD) simulations. A 90 picosecond (ps) trajectory for the protein dimer, 19 water molecules and two counter ions has been calculated at constant temperature. Effects of intermolecular interactions on the structure and dynamics have been studied. The exact crystallographic symmetry is lost and the atomic fluctuations differ from one monomer to the other. The average MD structure is more stable than the X-ray one, as judged by accessible surface area and energy calculations. Crystal (non-dimeric) interactions have been simulated in another 40 ps trajectory by using harmonic restraints to represent intermolecular hydrogen bonds. The conformational changes with respect to the X-ray structure are then virtually suppressed. The unrestrained dimer trajectory has been scanned for cooperative motions involving secondary structure elements. The intrinsic collective motions of the monomer are transmitted via intermolecular contacts to the dimer structure. The existence of a stable dimeric form of CTF, resembling the crystallographic one, has been documented. At the cost of fairly small energy expenditure the dimer has considerable conformational flexibility. This flexibility may endow the dimer with some functional potential as an energy transducer.

Aspects of model building applied to the C-terminal domain of the L12 protein from chloroplast ribosomes: a molecular dynamics study

A 170 picosecond molecular dynamics trajectory has been calculated starting from a model-built structure of chloroplast CTF. Local conformational changes occur during the equilibration period. Thereafter, a dynamically stable structure is attained. The conformational changes involve a turn connecting two structural subdomains which has an amino acid insertion and several substitutions with respect to the E. coli sequence. Potential energy minimisation alone fails to detect such a change. The overall folding and atomic positional fluctuations are very similar to those found in MD simulations of the E. coli molecule. The combined use of computer graphics based model building and MD calculations has lead to a thermally stable putative structure for the chloroplast CTF.

Molecular dynamics simulations of the holo and apo forms of retinol binding protein. Structural and dynamical changes induced by retinol removal

The effects of removing retinol from the X-ray structure of holo-retinol binding protein are studied using the molecular dynamics technique. Structural and dynamical properties emerging from an 80 ps simulation of the apo form, for which no crystallographic structure is available, are compared with the results of a 70 ps trajectory of the holo-protein. Dynamical stationarity is attained after roughly 30 ps, and the resulting average structure is proposed as a reasonable model of the apo-protein. Conformational changes are observed for the loops at the beta-barrel entrance during the non-equilibrium part of the apo-trajectory. Tryptophan labelling experiments and retinoid reconstitution experiments point towards this part of the molecule as being involved in prealbumin binding. Structural changes in this region may therefore explain the differences in prealbumin affinity between the apo and holo forms. Furthermore, a change in the position of the alpha-helix, corresponding to a pivot around its C terminus, is observed for the apo-protein. The resulting conformation of the alpha-helix is found to be similar to that in apo-beta-lactoglobulin, which also can bind retinol and for which a crystal structure exists. The results from the holo simulation are compared to the crystallographic data and show good agreement. The dynamics of the secondary and tertiary structural elements are analysed and compared for the two forms. The beta-barrel is found to be extremely cooperative in its atomic motions in both simulations, and the top and bottom beta-sheets perform collective fluctuations with respect to each other in the low-frequency limit of the simulations. The dynamics of the alpha-helical region presents clear differences between the two forms; while the holo-protein has a well-defined spectrum for the longitudinal stretching mode, the apo form displays a fairly large bending of the alpha-helix at several points of the trajectory.

Structural changes in retinol binding protein induced by retinol removal. A molecular dynamics study

Relationships between structure and function for retinol binding protein (RBP) are elucidated with help of a 2.0 A resolution X-ray structure of the holo-protein and an average molecular dynamics (MD) structure of the apo-form. Comparisons between MD simulations of both the apo- and holo-forms with the X-ray holo-structure show conformational changes in apo-RBP that may be functionally significant. The average three dimensional structure obtained for apo-RBP is compared to the related protein apo-beta-lactoglobulin. Available biochemical information is consistent with structure/function relationships derived here.

A molecular dynamics study of the C-terminal fragment of the L7/L12 ribosomal protein. Secondary structure motion in a 150 picosecond trajectory

A 150 picosecond molecular dynamics computer simulation of the C-terminal fragment of the L7/L12 ribosomal protein from Escherichia coli is reported. The molecular dynamics results are compared with the available high-resolution X-ray data in terms of atomic positions, distances and positional fluctuations. Good agreement is found between the molecular dynamics results and the X-ray data. The form and parameters of the interaction potential energy function and the procedures for deriving it are discussed. Some current misunderstandings concerning the ways of evaluating the efficiency of molecular dynamics algorithms and of application of bond-length constraints in protein simulations are cleared up. The 150 picosecond trajectory has been scanned in a search for correlated motions within and between secondary structure elements. The beta-strands have diffusional stretching modes, and uncorrelated transversal displacements. The dynamic analysis of alpha-helices shows a variety of features. The atomic fluctuations differ between the helix ends; this effect reflects long time-scale motions. Two alpha-helices, alpha A and alpha C, show diffusive longitudinal stretching modes. The third helix, alpha B, has a correlated asymmetric longitudinal stretching; the N-terminal part dominates this behaviour. Furthermore, alpha B presents a librational motion with respect to the other parts of the molecule with a frequency of approximately 5 cm-1. This motion is coupled to helix stretching. Interestingly, the regions of highly conserved residues contain the most mobile parts of the molecule.

The three-dimensional structure of retinol-binding protein

The complex of retinol with its carrier protein, retinol-binding protein (RBP) has been crystallized and its three-dimensional structure determined using X-ray crystallography. Its most striking feature is an eight-stranded up-and-down beta barrel core that completely encapsulates the retinol molecule. The retinol molecule lies along the axis of the barrel with the beta-ionone ring innermost and the tip of the isoprene tail close to the surface.