Quantum molecular dynamics simulation of proton transfer in cytochrome c oxidase

Proton transfer/translocation is studied in cytochrome c oxidase (CcO) by a combination of quantum mechanics (QM) for the transferring protons and classical molecular dynamics (MD) for the protein and solvent. The possibility of a glutamate, Glu286 in the Rhodobacter sphaeroides numbering scheme, acting as a rely point for proton translocation is investigated. The MD finds a hydrogen-bonded cycle of two waters and the carboxylate oxygens of Glu286. The possibility of protonating Glu286 to form neutral GluH is studied and we find that, as experimentally inferred, this glutamate can spend most of its time as GluH. Since translocation relies on the presence of water chains within CcO channels, MD is used to assess their formation. Glu286 and Mg(2+) can be connected by continuous hydrogen-bonded chains that are robust, though transient, and the protein appears spongy above (toward the outer membrane) the Mg(2+). In contrast, the D-channel spanning Asp132, close to the inner membrane surface, to Glu286, forms water chains that are much sparser and do not continuously connect these residues. Rather, there are chains spanning Glu286 to the vicinity of Asn140, and other more robust and ramified water structures that connect Asp132 with waters close to Asn140.     

Elucidation of the protonation states of the catalytic residues in mtKasA: implications for inhibitor design

KasA (β-ketoacyl ACP synthase I) is involved in the biosynthetic pathway of mycolic acids, an essential component of the cell wall in Mycobacterium tuberculosis. It was shown that KasA is essential for the survival of the pathogen and thus could serve as a new drug target for the treatment of tuberculosis. The active site of KasA was previously characterized by X-ray crystallography. However, questions regarding the protonation state of specific amino acids, the orientation of the histidine groups within the active site, and additional conformers being accessible at ambient temperatures remain open and have to be addressed prior to the design of new inhibitors. We investigate the active site of KasA in this work by means of structural motifs and relative energies. Molecular dynamics (MD) simulations, free energy perturbation computations, and calculations employing the hybrid quantum mechanics/molecular mechanics (QM/MM) method made it possible to determine the protonation status and reveal important details about the catalytic mechanism of KasA. Additionally, we can rationalize the molecular basis for the acyl-transfer activity in the H311A mutant. Our data strongly suggest that inhibitors should be able to inhibit different protonation states because the enzyme can switch easily between a zwitterionic and neutral state.     

A computational study of the deacylation mechanism of human butyrylcholinesterase

To investigate the mechanism of the deacylation reaction in the active site of human butyrylcholinesterase (BuChE), we carried out quantum mechanical (QM) calculations on cluster models of the active site built from a crystallographic structure. The models consisted of the substrate butyrate moiety, the catalytic triad of residues (Ser198, Glu325, and His438), the "oxy-anion hole" (Gly116, Gly117, and Ala199), the side chain of Glu197, four water molecules, the side chain of Ser225, and the peptide linkage between Val321 and Asn322. Analyses of the equilibrium geometries, electronic properties, and energies of the QM models gave insights into the catalytic mechanism. In addition, the QM calculations provided the data required to build a molecular mechanics representation of the reactive BuChE region that was employed in molecular dynamics simulations followed by molecular-mechanics-Poisson-Boltzmann (MM-PB) calculations. Subsequently, we combined the QM energies with average MM-PB energies to estimate the free energy of the reactive structures in the enzyme. The rate-determining step corresponds to the formation of a tetrahedral intermediate with a free-energy barrier of approximately 14.0 kcal/mol. The modulation of the BuChE activity, exerted by either neutral molecules (glycerol, GOL) or a second butyrylcholine (CHO) molecule bound to the cation-pi site, does not involve any significant allosteric effect. Interestingly, the presence of GOL or CHO stabilizes a product complex formed between a butyric acid molecule and BuChE. These results are in consonance with the crystallographic structure of BuChE, in which the catalytic Ser198 interacts with a butyric fragment, while the cation-pi site is occupied by one GOL molecule.     

QM/MM study of the active site of free papain and of the NMA-papain complex.

Hybrid quantum mechanical/molecular mechanical (QM/MM) calculations using restricted and unrestricted Hartree-Fock and B3LYP ab initio (QM) and Amber force field (MM), respectively, have been applied to study the catalytic site of papain in both free and substrate bonded forms. Ab initio geometry optimizations have been performed for the active site of papain and the N-methyl-acetamide (NMA)-papain complex within the molecular mechanical treatment of the protein environment. A covalent tetrahedral intermediate structure could be obtained only when the amide N atom of the substrate molecule was protonated through a proton transfer from the His-159 in the catalytic site. Our results support the previous assumption that a proton transfer from His-159 to the amide N atom of the substrate occurs prior to or concerted with the nucleophilic attack of the Cys-25 sulfur atom to the carbonyl group of the substrate. The electron correlation effect will reduce the proton transfer barrier. Therefore, this proton transfer can be easily observed in the B3LYP/6-31G* calculations. The HF/6-31G* method overestimates the reaction barrier against this proton transfer. The sulfur atom of Cys-25 and the imidazole ring of His-159 are found to be coplanar in the free form of the enzyme. However, the rotation of the imidazole ring of His-159 was observed during the formation of the tetrahedral intermediate. Without the papain environment, the coplanar thiolate-imidazolium ion pair RS-...ImH+ is much less stable than the neutral form of RSH....Im. Within the protein environment, however, the thiolate-imidazolium ion pair becomes more stable than its neutral form by 4.1 and 0.4 kcal/mol in HF/6-31G* and B3LYP/6-31G* calculations, respectively. The barrier of proton transfer from S-H group of Cys-25 to the imidazole ring of His-159 was reduced from 22.0 kcal/mol to 15.2 kcal/mol by the protein environment in HF/6-31G* calculations. This barrier is found to be much smaller (2.5 kcal/mol) in B3LYP/6-31G* calculations.     

Solvent Dielectric Effect and Side Chain Mutation on the Structural Stability of Burkholderia cepacia Lipase Active Site: A Quantum Mechanical/Molecular Mechanics Study

Quantum mechanical and molecular dynamics methods were used to analyze the structure and stability of neutral and zwitterionic configurations of the extracted active site sequence from a Burkholderia cepacia lipase, histidyl-seryl-glutamin (His86-Ser87-Gln88) and its mutated form, histidyl-cysteyl-glutamin (His86-Cys87-Gln88) in vacuum and different solvents. The effects of solvent dielectric constant, explicit and implicit water molecules and side chain mutation on the structure and stability of this sequence in both neutral and zwitterionic forms are represented. The quantum mechanics computations represent that the relative stability of zwitterionic and neutral configurations depends on the solvent structure and its dielectric constant. Therefore, in vacuum and the considered non-polar solvents, the neutral form of the interested sequences is more stable than the zwitterionic form, while their zwitterionic form is more stable than the neutral form in the aqueous solution and the investigated polar solvents in most cases. However, on the potential energy surfaces calculated, there is a barrier to proton transfer from the positively charged ammonium group to the negatively charged carboxylat group or from the ammonium group to the adjacent carbonyl oxygen and or from side chain oxygen and sulfur to negatively charged carboxylat group. Molecular dynamics simulations (MD) were also performed by using periodic boundary conditions for the zwitterionic configuration of the hydrated molecules in a box of water molecules. The obtained results demonstrated that the presence of explicit water molecules provides the more compact structures of the studied molecules. These simulations also indicated that side chain mutation and replacement of sulfur with oxygen leads to reduction of molecular flexibility and packing.     

Ionization properties of mixed lipid membranes: a Gouy-Chapman model of the electrostatic-hydrogen bond switch

The dissociation state of phosphatidic acid (PA) in a lipid bilayer is governed by the competition of proton binding and formation of a hydrogen bond through a mechanism termed the electrostatic-hydrogen bond switch. This mechanism has been suggested to provide the basis for the specific recognition of PA by proteins. Even in bare lipid bilayers the electrostatic-hydrogen bond switch is present if the membrane contains lipids like phosphatidylethanolamine that act as hydrogen bond donors. In this case, the dissociation state (pK(a)) of PA depends strongly on membrane composition. In the present work we incorporate the electrostatic-hydrogen bond switch mechanism into the Gouy-Chapman model for a membrane that is composed of PA and a hydrogen bond-donating zwitterionic lipid. To this end, our model integrates into the Gouy-Chapman approach a recently suggested electrostatic model for zwitterionic lipids. Hydrogen bond formation is incorporated phenomenologically as an additional non-electrostatic interaction between the phosphomonoester headgroup of PA and the zwitterionic lipid headgroup. We express the energetics of the composite membrane in terms of a free energy functional whose minimization leads to a modified non-linear Poisson-Boltzmann equation that we have solved numerically. Our calculations focus on the influence of the membrane environment on the dissociation state of PA. This influence can be expressed as a shift of the second pK(a) of PA, which we calculate as function of membrane composition, following experimental observation. The shift is large and negative if PA is the minor component in the membrane, and it changes over four pH units as function of the mole fraction of PA in the membrane. In contrast, the shift of the second pK(a) of PA remains small and is always positive if the zwitterionic lipid is unable to act as hydrogen bond donor. Hence, we find that the electrostatic-hydrogen bond switch mechanism regulates the dissociation state of PA with much greater sensitivity than would be possible based on a pure electrostatic regulation through the membrane potential.     

Chromophore protonation state controls photoswitching of the fluoroprotein asFP595

Fluorescent proteins have been widely used as genetically encodable fusion tags for biological imaging. Recently, a new class of fluorescent proteins was discovered that can be reversibly light-switched between a fluorescent and a non-fluorescent state. Such proteins can not only provide nanoscale resolution in far-field fluorescence optical microscopy much below the diffraction limit, but also hold promise for other nanotechnological applications, such as optical data storage. To systematically exploit the potential of such photoswitchable proteins and to enable rational improvements to their properties requires a detailed understanding of the molecular switching mechanism, which is currently unknown. Here, we have studied the photoswitching mechanism of the reversibly switchable fluoroprotein asFP595 at the atomic level by multiconfigurational ab initio (CASSCF) calculations and QM/MM excited state molecular dynamics simulations with explicit surface hopping. Our simulations explain measured quantum yields and excited state lifetimes, and also predict the structures of the hitherto unknown intermediates and of the irreversibly fluorescent state. Further, we find that the proton distribution in the active site of the asFP595 controls the photochemical conversion pathways of the chromophore in the protein matrix. Accordingly, changes in the protonation state of the chromophore and some proximal amino acids lead to different photochemical states, which all turn out to be essential for the photoswitching mechanism. These photochemical states are (i) a neutral chromophore, which can trans-cis photoisomerize, (ii) an anionic chromophore, which rapidly undergoes radiationless decay after excitation, and (iii) a putative fluorescent zwitterionic chromophore. The overall stability of the different protonation states is controlled by the isomeric state of the chromophore. We finally propose that radiation-induced decarboxylation of the glutamic acid Glu215 blocks the proton transfer pathways that enable the deactivation of the zwitterionic chromophore and thus leads to irreversible fluorescence. We have identified the tight coupling of trans-cis isomerization and proton transfers in photoswitchable proteins to be essential for their function and propose a detailed underlying mechanism, which provides a comprehensive picture that explains the available experimental data. The structural similarity between asFP595 and other fluoroproteins of interest for imaging suggests that this coupling is a quite general mechanism for photoswitchable proteins. These insights can guide the rational design and optimization of photoswitchable proteins.     

Acid-triggered membrane insertion of Pseudomonas exotoxin A involves an original mechanism based on pH-regulated tryptophan exposure

Exposure to low endosomal pH during internalization of Pseudomonas exotoxin A (PE) triggers membrane insertion of its translocation domain. This process is a prerequisite for PE translocation to the cytosol where it inactivates protein synthesis. Although hydrophobic helices enable membrane insertion of related bacterial toxins such as diphtheria toxin, the PE translocation domain is devoid of hydrophobic stretches and the structural features triggering acid-induced membrane insertion of PE are not known. Here we have identified a molecular device that enables PE membrane insertion. This process is promoted by exposure of a key tryptophan residue. At neutral pH, this Trp is buried in a hydrophobic pocket closed by the smallest alpha-helix of the translocation domain. Upon acidification, protonation of the Asp that is the N-cap residue of the helix leads to its destabilization, enabling Trp side chain insertion into the endosome membrane. This tryptophan-based membrane insertion system is surprisingly similar to the membrane-anchoring mechanism of human annexin-V and could be used by other proteins as well.     

Theoretical studies of the ATP hydrolysis mechanism of myosin

The ATP hydrolysis mechanism of myosin was studied using quantum chemical (QM) and molecular dynamics calculations. The initial model compound for QM calculations was constructed on the basis of the energy-minimized structure of the myosin(S1dc)-ATP complex, which was determined by molecular mechanics calculations. The result of QM calculations suggested that the ATP hydrolysis mechanism of myosin consists of a single elementary reaction in which a water molecule nucleophilically attacked gamma-phosphorus of ATP. In addition, we performed molecular dynamics simulations of the initial and final states of the ATP hydrolysis reaction, that is, the myosin-ATP and myosin-ADP.Pi complexes. These calculations revealed roles of several amino acid residues (Lys185, Thr186, Ser237, Arg238, and Glu459) in the ATPase pocket. Lys185 maintains the conformation of beta- and gamma-phosphate groups of ATP by forming the hydrogen bonds. Thr186 and Ser237 are coordinated to a Mg(2+) ion, which interacts with the phosphates of ATP and therefore contributes to the stabilization of the ATP structure. Arg238 and Glu459, which consisted of the gate of the ATPase pocket, retain the water molecule acting on the hydrolysis at the appropriate position for initiating the hydrolysis.     

Passive Diffusion of Ciprofloxacin and its Metalloantibiotic: A Computational and Experimental study

Fluoroquinolones (FQ) are antibiotics widely used in clinical practise, but the development of bacterial resistance to these drugs is currently a critical public health problem. In this context, ternary copper complexes of FQ (CuFQPhen) have been studied as a potential alternative. In this study, we compared the passive diffusion across the lipid bilayer of one of the most used FQ, ciprofloxacin (Cpx), and its ternary copper complex, CuCpxPhen, that has shown previous promising results regarding antibacterial activity and membrane partition. A combination of spectroscopic studies and molecular dynamics simulations were used and two different model membranes tested: one composed of anionic phospholipids, and the other composed of zwitterionic phospholipids. The obtained results showed a significantly higher membrane permeabilization activity, larger partition, and a more favourable free energy landscape for the permeation of CuCpxPhen across the membrane, when compared to Cpx. Furthermore, the computational results indicated a more favourable translocation of CuCpxPhen across the anionic membrane, when compared to the zwitterionic one, suggesting a higher specificity towards the former. These findings are important to decipher the influx mechanism of CuFQPhen in bacterial cells, which is crucial for the ultimate use of CuFQPhen complexes as an alternative to FQ to tackle multidrug-resistant bacteria.     

Exploring the RING-Catalyzed Ubiquitin Transfer Mechanism by MD and QM/MM Calculations

Ubiquitylation is a universal mechanism for controlling cellular functions. A large family of ubiquitin E3 ligases (E3) mediates Ubiquitin (Ub) modification. To facilitate Ub transfer, RING E3 ligases bind both the substrate and ubiquitin E2 conjugating enzyme (E2) linked to Ub via a thioester bond to form a catalytic complex. The mechanism of Ub transfer catalyzed by RING E3 remains elusive. By employing a combined computational approach including molecular modeling, molecular dynamics (MD) simulations, and quantum mechanics/molecular mechanics (QM/MM) calculations, we characterized this catalytic mechanism in detail. The three-dimensional model of dimeric RING E3 ligase RNF4 RING, E2 ligase UbcH5A, Ub and the substrate SUMO2 shows close contact between the substrate and Ub transfer catalytic center. Deprotonation of the substrate lysine by D117 on UbcH5A occurs with almost no energy barrier as calculated by MD and QM/MM calculations. Then, the side chain of the activated lysine gets close to the thioester bond via a conformation change. The Ub transfer pathway begins with a nucleophilic addition that forms an oxyanion intermediate of a 4.23 kcal/mol energy barrier followed by nucleophilic elimination, resulting in a Ub modified substrate by a 5.65 kcal/mol energy barrier. These results provide insight into the mechanism of RING-catalyzed Ub transfer guiding the discovery of Ub system inhibitors.     

Conformational preferences of hypermodified nucleoside lysidine (k2C) occurring at "wobble" position in anticodon loop of tRNA(Ile)

Conformational preferences of hypermodified nucleoside, 4-amino-2-(N(6)-lysino)-1-(beta-D-ribofuranosyl) pyrimidinium (Lysidine or 2-lysyl cytidine), usually designated as k(2)C, have been investigated theoretically by the quantum chemical perturbative configuration interaction with localized orbitals (PCILO) method. The zwitterionic, non-zwitterionic, neutral, and tautomeric forms have been studied. Automated geometry optimization using molecular mechanics force field (MMFF), semi-empirical quantum chemical PM3, and ab initio molecular orbital Hartree-Fock SCF quantum mechanical calculations have also been made to compare the salient features. The predicted most stable conformations of zwitterionic, non-zwitterionic, neutral, and tautomeric form are such that in each of these molecules the orientation of lysidine moiety (R) is trans to the N(1) of cytidine. The preferred base orientation is anti (chi = 3 degrees ) and the lysine substituent folds back toward the ribose ring. This results in hydrogen bonding between the carboxyl oxygen O(12a) of lysine moiety and the 2'-hydroxyl group of ribose sugar. In all these four forms of lysidine O(12a)...H-C(9) and O(12b)...H-N(11) interactions provide stability to respective stable conformers. Watson-Crick base pairing of lysidine with A is feasible only with the tautomeric form of usual anti oriented lysidine. This can help in recognition of AUA codon besides in avoiding misrecognition of AUG.     

Conformational analysis of Amphotericin B

Within a theoretical approach to the problem of antifungal action of Amphotericin B (AmB), a conformational analysis of the neutral and zwitterionic form of this antibiotic in vacuo was performed by the MM2P and AM1 methods. The analysis was carried out with regard to the mutual orientation of the macrolidic and glycosidic fragments of the molecule, which is defined by the phi and psi steric angles. This orientation defines the overall shape of the molecule and is postulated to be important for the antifungal action of the drug. As a result of the MM2P calculations, phi, psi steric energy and population maps were prepared. Several conformers were found on these maps but only two of them (one each for the zwitterionic and the neutral forms of the antibiotic) were previously observed experimentally for isolated molecules. Our other calculated conformers were not observed experimentally but we propose that they may also appear in the AmB channel structure. The results of our conformational analysis were compared with experimental NMR data (nuclear Overhauser effects between selected hydrogen atoms) obtained previously. New structural information obtained for AmB in the present work will be useful for building a molecular model of AmB-target interactions as well as for designing new derivatives of AmB.     

An Alternative Mechanism for the Methylation of Phosphoethanolamine Catalyzed by Plasmodium falciparum Phosphoethanolamine Methyltransferase

The phosphobase methylation pathway catalyzed by the phosphoethanolamine methyltransferase in Plasmodium falciparum (PfPMT), the malaria parasite, offers an attractive target for anti-parasitic drug development. PfPMT methylates phosphoethanolamine (pEA) to phosphocholine for use in membrane biogenesis. Quantum mechanics and molecular mechanics (QM/MM) calculations tested the proposed reaction mechanism for methylation of pEA involving the previously identified Tyr-19–His-132 dyad, which indicated an energetically unfavorable mechanism. Instead, the QM/MM calculations suggested an alternative mechanism involving Asp-128. The reaction coordinate involves the stepwise transfer of a proton to Asp-128 via a bridging water molecule followed by a typical S_n2-type methyl transfer from S-adenosylmethionine to pEA. Functional analysis of the D128A, D128E, D128Q, and D128N PfPMT mutants shows a loss of activity with pEA but not with the final substrate of the methylation pathway. X-ray crystal structures of the PfPMT-D128A mutant in complex with S-adenosylhomocysteine and either pEA or phosphocholine reveal how mutation of Asp-128 disrupts a hydrogen bond network in the active site. The combined QM/MM, biochemical, and structural studies identify a key role for Asp-128 in the initial step of the phosphobase methylation pathway in Plasmodium and provide molecular insight on the evolution of multiple activities in the active site of the PMT.     

Mena–GRASP65 interaction couples actin polymerization to Golgi ribbon linking

In mammalian cells, the Golgi reassembly stacking protein 65 (GRASP65) has been implicated in both Golgi stacking and ribbon linking by forming trans-oligomers through the N-terminal GRASP domain. Because the GRASP domain is globular and relatively small, but the gaps between stacks are large and heterogeneous, it remains puzzling how GRASP65 physically links Golgi stacks into a ribbon. To explore the possibility that other proteins may help GRASP65 in ribbon linking, we used biochemical methods and identified the actin elongation factor Mena as a novel GRASP65-binding protein. Mena is recruited onto the Golgi membranes through interaction with GRASP65. Depleting Mena or disrupting actin polymerization resulted in Golgi fragmentation. In cells, Mena and actin were required for Golgi ribbon formation after nocodazole washout; in vitro, Mena and microfilaments enhanced GRASP65 oligomerization and Golgi membrane fusion. Thus Mena interacts with GRASP65 to promote local actin polymerization, which facilitates Golgi ribbon linking.     

Reversible intracellular translocation of KRas but not HRas in hippocampal neurons regulated by Ca2+/calmodulin

The Ras/MAPK pathway regulates synaptic plasticity and cell survival in neurons of the central nervous system. Here, we show that KRas, but not HRas, acutely translocates from the plasma membrane (PM) to the Golgi complex and early/recycling endosomes in response to neuronal activity. Translocation is reversible and mediated by the polybasic-prenyl membrane targeting motif of KRas. We provide evidence that KRas translocation occurs through sequestration of the polybasic-prenyl motif by Ca2+/calmodulin (Ca2+/CaM) and subsequent release of KRas from the PM, in a process reminiscent of GDP dissociation inhibitor-mediated membrane recycling of Rab and Rho GTPases. KRas translocation was accompanied by partial intracellular redistribution of its activity. We conclude that the polybasic-prenyl motif acts as a Ca2+/CaM-regulated molecular switch that controls PM concentration of KRas and redistributes its activity to internal sites. Our data thus define a novel signaling mechanism that differentially regulates KRas and HRas localization and activity in neurons.     

Insights into the structure and dynamics of the dinuclear zinc beta-lactamase site from Bacteroides fragilis

Herein, we report quantum chemical calculations and molecular dynamics (MD) simulations of the dinuclear form of the Bacteroides fragilis zinc beta-lactamase. We studied four different configurations which differ in the protonation state of the Asp103 residue and in the presence or absence of a Zn1-OH-Zn2 bridge. The flexibility of the Zn1-OH-Zn2 bridge was studied by means of quantum mechanical (QM) calculations on cluster models while the relative stabilities of the different configurations were estimated from QM linear scaling calculations on the enzyme. Contacts between important residues (Cys104, Asp69, Lys185, etc.), the solvation of the zinc ions, and the conformation of the active site beta-hairpin loop were characterized by the MD analyses. The influence of the buried sodium ion close to the Zn2 position was investigated by carrying out a secondary simulation where the sodium ion was replaced with an internal water molecule. The comparative structural analyses among the different MD trajectories augmented with energetic calculations have demonstrated that the B. fragilis protein efficiently binds the internal Na(+) ion observed crystallographically. Moreover, we found that when Asp103 is unprotonated, a rigid Zn1-OH-Zn2 bridge results, while for neutral Asp103, a fluctuating Zn1-Zn2 distance was possible via the breaking and formation of the Zn1-OH-Zn2 bridge. The mechanistic implications of these observations are discussed in detail.     

Characterization of a catalytic ligand bridging metal ions in phosphodiesterases 4 and 5 by molecular dynamics simulations and hybrid quantum mechanical/molecular mechanical calculations

Cyclic nucleotide phosphodiesterases (PDEs) constitute a large superfamily of enzymes regulating concentrations of intracellular second messengers cAMP and cGMP through PDE-catalyzed hydrolysis. Although three-dimensional x-ray crystal structures of PDE4 and PDE5 have been reported, it is uncertain whether a critical, second bridging ligand (BL2) in the active site is H2O or HO- because hydrogen atoms cannot be determined by x-ray diffraction. The identity of BL2 is theoretically determined by performing molecular dynamics simulations and hybrid quantum mechanical/molecular mechanical (QM/MM) calculations, for the first time, on the protein structures resolved by x-ray diffraction. The computational results confirm our previous suggestion, which was based on QM calculations on a simplified active site model, that BL2 in PDE4 should be HO-, rather than H2O, serving as the nucleophile to initialize the catalytic hydrolysis of cAMP. The molecular dynamics simulations and QM/MM calculations on PDE5 demonstrate for the first time that the BL2 in PDE5 should also be HO- rather than H2O as proposed in recently published reports on the x-ray crystal structures, which serves as the nucleophile to initialize the PDE5-catalyzed hydrolysis of cGMP. These fundamental structural insights provide a rational basis for future structure-based drug design targeting PDEs.     

Molecular dynamics simulations of Leu-enkephalin in water and DMSO.

The structure of Leu-enkephalin (L-Enk) and Met-enkephalin (M-Enk) have frequently been studied, in particular by nuclear magnetic resonance spectroscopy. After more than 20 years of research, it was concluded that enkephalins have no preferred structure in aqueous solution, but that they may have in other solvents. We have performed molecular dynamics simulations of zwitterionic L-Enk in water, and zwitterionic as well as neutral L-Enk dimethyl sulfoxide (DMSO). In water the peptide is very flexible, although there seems to be a preference for compact conformations. In DMSO, the peptide forms a clear salt bridge in the zwitterionic form, but has no preferred conformation in the neutral form. This difference in conformation may provide an explanation for measurements in DMSO in which multiple conformations were found to exist. In this paper we introduce a new formulation for a dihedral angle autocorrelation function, and apply it to study side-chain dynamics in L-Enk. We find that the side-chain dynamics of the large Tyr and Phe residues cannot be adequately sampled in 2.0-ns simulations, while this does seem to be possible for the smaller Leu side chain.     

Cover Picture

Fibroblast Growth Factor 2 (FGF2) is a potent cell survival factor involved in tumor induced angiogenesis. FGF2 is secreted from cells through an unconventional secretory mechanism that is based upon direct translocation across the plasma membrane. The molecular mechanism underlying this process depends on a surprisingly small set of trans‐acting factors and cis‐elements in FGF2. With the reconstitution of the core mechanism of FGF2 membrane translocation using purified components, the molecular basis of FGF2 secretion from tumor cells has been revealed.