PubMed Health⌕ Search

Biomedical subjects

Markus Meuwly

Publications and source records attributed to Markus Meuwly.

At least 19 recordsLinked to original sources
PubMed2007

Energetics and dynamics in MbCN: CN--vibrational relaxation from molecular dynamics simulations.

The dynamics of the cyanide anion bound to sperm-whale myoglobin is investigated using atomistic simulations. With density-functional theory, a 2D potential energy surface for the cyanide-heme complex is calculated. Two deep minima with a stabilization energy of approximately 50 kcal/mol corresponding to two different binding orientations (Fe-CN and Fe-NC) of the ligand are found. The Fe-CN conformation is favored over Fe-NC by several kcal/mol. Mixed quantum mechanics/molecular mechanics calculations show that the binding orientation affects the bond strength of the ligand, with a significantly different bond length and a 25 cm-1 shift in the fundamental CN-frequency. For the molecular dynamics (MD) simulations, a 3-center fluctuating charge model for the Fe-CN unit is developed that captures polarization and ligand-metal charge transfer. Stability arguments based on the energetics around the active site and the CN- frequency shifts suggest that the Fe-CN conformation with epsilon-protonation of His epsilon 64 are most likely, which is in agreement with experiment. Both equilibrium and nonequilibrium MD simulations are carried out to investigate the relaxation time scale and possible relaxation pathways in bound MbCN. The nonequilibrium MD simulations with a vibrationally excited ligand reveal that vibrational relaxation takes place on a time scale of hundreds of picoseconds within the active site. This finding supports the hypothesis that the experimentally observed relaxation rate (3.6 ps) reflects the repopulation of the electronic ground state.

Animals↗
PubMed2006

Importance of individual side chains for the stability of a protein fold: computational alanine scanning of the insulin monomer.

A new computational approach is proposed to probe the importance of residue side chains for the stability of a protein fold. Computational mutations to estimate protein stability (CMEPS) is based on the notion that the binding free energy corresponding to the complexation of a given side chain, considered as a "pseudo-ligand" of the wild type protein, reflects the importance of this side chain to the thermodynamic stability of the protein. The contribution of a particular side chain to the folding energy is estimated according to the molecular mechanics-generalized born surface area MM-GBSA approach, using a single molecular dynamics simulation trajectory of the wild type protein. CMEPS is a first principles method which does not contain any adjustable parameter that could be fitted to experimental data. The approach is first validated for Barnase and the B1 domain of protein L, for which a correlation coefficient R = 0.73, between experimental and CMEPS calculated DeltaDeltaG values, is found and then applied to the insulin monomer. In the present application, CMEPS replaces each amino acid by an alanine residue. Therefore, most mutations lead to cavities in the protein. From this the change in stability can be correlated with increased cavity volume. For insulin, this correlation is very similar compared with data previously analyzed for T4 lysozyme from an experiment for buried apolar side chains. There, the increased cavity volume has been related to the hydrophobic effect. However, since CMEPS uses the energetics in terms of electrostatic and van der Waals interactions (and not the hydrophobic effect which is difficult to relate to physical interactions), it is possible to study the effect of mutations of polar and solvent accessible side chains. According to CMEPS, residues Leu A16, Tyr A19, Leu B11, Leu B15, and Arg B22 are most important for the stability of the monomeric insulin fold. This is in agreement with experimental data. As a consequence, mutation of these residues may lead to misfolded and inactive insulin analogues.

Alanine↗
PubMed2006

How inaccuracies in protein structure models affect estimates of protein-ligand interactions: computational analysis of HIV-I protease inhibitor binding.

The influence of possible inaccuracies that can arise during homology modeling of protein structures used for ligand binding studies were investigated with the molecular mechanics generalized Born surface area (MM-GBSA) method. For this, a family of well-characterized HIV-I protease-inhibitor complexes was used. Validation of MM-GBSA led to a correlation coefficient ranging from 0.72 to 0.93 between calculated and experimental binding free energies DeltaG. All calculated DeltaG values were based on molecular dynamics simulations with explicit solvent. Errors introduced into the protein structure through misplacement of side-chains during rotamer modeling led to a correlation coefficient between DeltaG(calc) and DeltaG(exp) of 0.75 compared with 0.90 for the correctly placed side chains. This is in contrast to homology models for members of the retroviral protease family with template structures ranging in sequence identity between 32% and 51%. For these protein models, the correlation coefficients vary between 0.84 and 0.87, which is considerably closer to the original protein (0.90). It is concluded that HIV-I low sequence identity with the template structure still allows creating sufficiently reliable homology models to be used for ligand-binding studies, although placement of the rotamers is a critical step during the modeling.

Amino Acid Sequence↗
PubMed2006

Allosteric control of cyclic di-GMP signaling.

Cyclic di-guanosine monophosphate is a bacterial second messenger that has been implicated in biofilm formation, antibiotic resistance, and persistence of pathogenic bacteria in their animal host. Although the enzymes responsible for the regulation of cellular levels of c-di-GMP, diguanylate cyclases (DGC) and phosphodiesterases, have been identified recently, little information is available on the molecular mechanisms involved in controlling the activity of these key enzymes or on the specific interactions of c-di-GMP with effector proteins. By using a combination of genetic, biochemical, and modeling techniques we demonstrate that an allosteric binding site for c-di-GMP (I-site) is responsible for non-competitive product inhibition of DGCs. The I-site was mapped in both multi- and single domain DGC proteins and is fully contained within the GGDEF domain itself. In vivo selection experiments and kinetic analysis of the evolved I-site mutants led to the definition of an RXXD motif as the core c-di-GMP binding site. Based on these results and based on the observation that the I-site is conserved in a majority of known and potential DGC proteins, we propose that product inhibition of DGCs is of fundamental importance for c-di-GMP signaling and cellular homeostasis. The definition of the I-site binding pocket provides an entry point into unraveling the molecular mechanisms of ligand-protein interactions involved in c-di-GMP signaling and makes DGCs a valuable target for drug design to develop new strategies against biofilm-related diseases.

Allosteric Site↗
PubMed2006

New experimental and theoretical results for argon broadening and shift of HCO+ rotational lines.

An experimental and theoretical study of the pressure broadening and the pressure shift of three HCO(+) rotational lines (j=4<--3, 5<--4 and 6<--5) perturbed by collisions with Ar is presented. The measurements are carried out at 77 K and are compared to close-coupling calculations performed on an accurate potential energy surface for the Ar-HCO(+) interaction extending from small to very large separations between the ion and the perturber. For the pressure broadening, agreement between experiment and theory is satisfactory for both close-coupling and semiclassical calculations. For the pressure shift, however, close-coupling calculations are superior. The results agree with experiment in sign and order of magnitude, while semiclassical calculations are inaccurate for the shift of the presently studied lines because they neglect the contribution of strong collisions.

Journal Article↗
PubMed2005

Structures and dynamics of protonated ammonia clusters.

The structures and infrared spectra of protonated ammonia clusters NH(4+)(NH3)n, for n < or = 8, are investigated using density functional-theory (DFT) calculations and semiempirical DFT/molecular dynamics simulations. For n < 5 the clusters are found to be mostly stable up to 100 K, while the larger clusters (n > or = 5) isomerize. Temperature effects are taken into account by performing ab initio molecular dynamics simulations with the computationally tractable self-consistent charges density functional tight-binding method. The infrared spectra at 10 K for the most stable isomers for n = 3-8 compare qualitatively with predissociation experiments, and using a common scaling factor almost quantitative agreement is found. For n > or = 6 the notion of multiple isomers present under the experimental conditions is supported. Of the 13 stable structures for n = 8 only three are found to survive at 100 K. All other clusters isomerize. Cluster structures are inferred from the analysis of the cumulative radial distribution function of the ammonia molecules surrounding the NH(4+) core. The infrared spectra are found to be typical for the structure of the clusters, which should help to relate the experimentally measured infrared spectra to the number and identity of the contributing isomers. For clusters that reorganize to a more stable isomer during the dynamics, the infrared spectrum is generally similar to that of the stable isomer itself. The clusters are found to preferably form globular structures, although chain-like arrangements are also among the low-energy configurations.

Journal Article↗
PubMed2005

Studying reactive processes with classical dynamics: rebinding dynamics in MbNO.

A new surface-crossing algorithm suitable for describing bond-breaking and bond-forming processes in molecular dynamics simulations is presented. The method is formulated for two intersecting potential energy manifolds which dissociate to different adiabatic states. During simulations, crossings are detected by monitoring an energy criterion. If fulfilled, the two manifolds are mixed over a finite number of time steps, after which the system is propagated on the second adiabat and the crossing is carried out with probability one. The algorithm is extensively tested (almost 0.5 mus of total simulation time) for the rebinding of NO to myoglobin. The unbound surface (Fe...NO) is represented using a standard force field, whereas the bound surface (Fe-NO) is described by an ab initio potential energy surface. The rebinding is found to be nonexponential in time, in agreement with experimental studies, and can be described using two time constants. Depending on the asymptotic energy separation between the manifolds, the short rebinding timescale is between 1 and 9 ps, whereas the longer timescale is about an order of magnitude larger. NO molecules which do not rebind within 1 ns are typically found in the Xenon-4 pocket, indicating the high affinity of NO to this region in the protein.

Algorithms↗
PubMed2005

Potential energy surface and molecular dynamics of MbNO: existence of an unsuspected FeON minimum.

Ligands such as CO, O(2), or NO are involved in the biological function of myoglobin. Here we investigate the energetics and dynamics of NO interacting with the Fe(II) heme group in native myoglobin using ab initio and molecular dynamics simulations. At the global minimum of the ab initio potential energy surface (PES), the binding energy of 23.4 kcal/mol and the Fe-NO structure compare well with the experimental results. Interestingly, the PES is found to exhibit two minima: There exists a metastable, linear Fe-O-N minimum in addition to the known, bent Fe-N-O global minimum conformation. Moreover, the T-shaped configuration is found to be a saddle point, in contrast to the corresponding minimum for NO interacting with Fe(III). To use the ab initio results for finite temperature molecular dynamics simulations, an analytical function was fitted to represent the Fe-NO interaction. The simulations show that the secondary minimum is dynamically stable up to 250 K and has a lifetime of several hundred picoseconds at 300 K. The difference in the topology of the heme-NO PES from that assumed previously (one deep, single Fe-NO minimum) suggests that it is important to use the full PES for a quantitative understanding of this system. Why the metastable state has not been observed in the many spectroscopic studies of myoglobin interacting with NO is discussed, and possible approaches to finding it are outlined.

Iron↗
PubMed2005

Study of the insulin dimerization: binding free energy calculations and per-residue free energy decomposition.

A calculation of the binding free energy for the dimerization of insulin has been performed using the molecular mechanics-generalized Born surface area approach. The calculated absolute binding free energy is -11.9 kcal/mol, in approximate agreement with the experimental value of -7.2 kcal/mol. The results show that the dimerization is mainly due to nonpolar interactions. The role of the hydrogen bonds between the 2 monomers appears to give the direction of the interactions. A per-atom decomposition of the binding free energy has been performed to identify the residues contributing most to the self association free energy. Residues B24-B26 are found to make the largest favorable contributions to the dimerization. Other residues situated at the interface between the 2 monomers were found to make favorable but smaller contributions to the dimerization: Tyr B16, Val B12, and Pro B28, and to an even lesser extent, Gly B23. The energy decomposition on a per-residue basis is in agreement with experimental alanine scanning data. The results obtained from a single trajectory (i.e., the dimer trajectory is also used for the monomer analysis) and 2 trajectories (i.e., separate trajectories are used for the monomer and dimer) are similar.

Crystallography, X-Ray↗
PubMed2005

Free-energy barriers in MbCO rebinding.

The rebinding of CO to myoglobin (Mb) from locations around the active site is studied using a combination of molecular dynamics and stochastic simulations for native and L29F mutant Mb. The interaction between the dissociated ligand and the protein environment is described by the recently developed fluctuating three-point charge model for the CO molecule. Umbrella sampling along trajectories, previously found to sample the binding site (B) and the Xe4 pocket, is used to construct free-energy profiles for the ligand escape. On the basis of the Smoluchowski equation, the relaxation of different initial population distributions is followed in space and time. For native Mb at room temperature, the calculated rebinding times are in good agreement with experimental values and give an inner barrier of 4.3 kcal/mol between the docking site B (Mb...CO) and the A state (bound MbCO), compared to an effective barrier, Heff, of 4.5 kcal/mol and barriers into the majority conformation A1 and the minority conformation A3 of 2.4 and 4.3 kcal/mol, respectively. In the case of the L29F mutant, the free-energy surface is flatter and the dynamics is much more rapid. As was found in experiment, escape to the Xe4 pocket is facile for L29F whereas, for native Mb, the barriers to this site are larger. At lower temperatures, the rebinding dynamics is delayed by orders of magnitude also due to increased barriers between the docking sites.

Binding Sites↗
PubMed2004

Ligand dynamics in myoglobin: calculation of infrared spectra for photodissociated NO.

We present molecular dynamics simulations of the photodissociated state of MbNO performed at 300 K using a fluctuating charge model for the nitric oxide (NO) ligand. After dissociation, NO is observed to remain mainly in the centre of the distal haem pocket, although some movement towards the primary docking site and the xenon-4 pocket can be seen. We calculate the NO infrared spectrum for the photodissociated ligand within the haem pocket and find a narrow peak in the range 1915-1922 cm(-1). The resulting blue shift of 1 to 8 cm(-1) compared to gas-phase NO is much smaller than the red shifts calculated and observed for carbon monoxide (CO) in Mb. A small splitting, due to NO in the xenon-4 pocket, is also observed. At lower temperatures, the spectra and conformational space explored by the ligand remain largely unchanged, but the electrostatic interactions with residue His64 become increasingly significant in determining the details of the ligand orientation within the distal haem pocket. The investigation of the effect of the L29F mutation reveals significant differences between the behaviour of NO and that of CO, and suggests a coupling between the ligand and the protein dynamics due to the different ligand dipole moments.

Amino Acid Substitution↗
PubMed2004

A combined experimental and computational study of dihydrido(phosphinooxazoline)iridium complexes.

The reaction of a [(PHOX)Ir(COD)](+) complex (COD = 1,5-cyclooctadiene) with dihydrogen was studied by NMR spectroscopy (PHOX = chiral phosphinooxazoline ligand). A single [(PHOX)Ir(H)(2)(COD)](+) isomer was formed as the primary product at -40 degrees C in THF. Subsequent reaction with H(2) at -40 to 0 degrees C led to a mixture of two diastereomeric [(PHOX)Ir(H)(2)(solvent)(2)](+) complexes with concomitant loss of cyclooctane. The stereochemistry of the three hydride complexes could be assigned from the NMR data. The structures and energies of the observed hydride complexes and the possible stereoisomers were calculated using density functional theory. The substantial energy differences (up to 39 kcal/mol) between the various stereoisomers demonstrate the strong influence of the chiral ligand. The observed stereoselective formation of dihydride complexes can be explained by steric effects of the PHOX ligand combined with a strong electronic influence of the coordinating N and P atoms, favoring addition of a hydride trans to the Ir-N bond.

Journal Article↗
PubMed2004

A comparison of the dynamic behavior of monomeric and dimeric insulin shows structural rearrangements in the active monomer.

Molecular dynamics (MD) simulations (5-10ns in length) and normal mode analyses were performed for the monomer and dimer of native porcine insulin in aqueous solution; both starting structures were obtained from an insulin hexamer. Several simulations were done to confirm that the results obtained are meaningful. The insulin dimer is very stable during the simulation and remains very close to the starting X-ray structure; the RMS fluctuations calculated from the MD simulation agree with the experimental B-factors. Correlated motions were found within each of the two monomers; they can be explained by persistent non-bonded interactions and disulfide bridges. The correlated motions between residues B24 and B26 of the two monomers are due to non-bonded interactions between the side-chains and backbone atoms. For the isolated monomer in solution, the A chain and the helix of the B chain are found to be stable during 5ns and 10ns MD simulations. However, the N-terminal and the C-terminal parts of the B chain are very flexible. The C-terminal part of the B chain moves away from the X-ray conformation after 0.5-2.5ns and exposes the N-terminal residues of the A chain that are thought to be important for the binding of insulin to its receptor. Our results thus support the hypothesis that, when monomeric insulin is released from the hexamer (or the dimer in our study), the C-terminal end of the monomer (residues B25-B30) is rearranged to allow binding to the insulin receptor. The greater flexibility of the C-terminal part of the beta chain in the B24 (Phe-->Gly) mutant is in accord with the NMR results. The details of the backbone and side-chain motions are presented. The transition between the starting conformation and the more dynamic structure of the monomers is characterized by displacements of the backbone of Phe B25 and Tyr B26; of these, Phe B25 has been implicated in insulin activation.

Amino Acid Sequence↗
PubMed2004

Selective hyperfine excitation of N2H+ by He: potential energy surface, cross sections, and propensity rules.

We present potential energy surfaces for the He-N2H+ system adiabatically corrected for the zero-point motion along the intermolecular stretching vibrations v1 = 0 and v1 = 1. The potentials are extended to shorter He-N2H+ separations which makes them useful for scattering calculations. Close coupling calculations of the spinless S matrices for the rotational excitation of N2H+ by He are presented, and recoupling techniques to obtain collisional excitation cross sections between the N2H+ hyperfine levels are used. The propensity rules between hyperfine levels are investigated for the case where two nuclear spins are involved. It is found that the only well defined propensity rule is DeltaF = DeltaF1 = Deltaj and that calculations are required in order to obtain the relative intensities of the two-spin hyperfine cross sections.

Journal Article↗
PubMed2004

Double proton transfer in the isolated and DNA-embedded guanine-cytosine base pair.

The energetics and dynamics of double proton transfer (DPT) is investigated theoretically for the Watson-Crick conformation of the guanine-cytosine (GC) base pair. Using semiempirical density functional theory the isolated and DNA-embedded GC pair is considered. Differences in the energetics and dynamics of DPT thus addresses the question of how relevant studies of isolated base pairs are for the understanding of processes occurring in DNA. Two-dimensional potential energy surfaces involving the transferring hydrogen atoms and the proton donors and acceptors are presented for both systems. The DPT reaction is accompanied by a contraction of the distance between the two bases with virtually identical energetic barriers being 18.8 and 18.7 kcal/mol for the isolated and DNA-embedded system, respectively. However, the transition state for DPT in the DNA-embedded GC pair is offset by 0.1 A to larger N-H separation compared to the isolated GC pair. Using activated ab initio molecular dynamics, DPT is readily observed for the isolated base pair with a minimal amount of 21.4 kcal/mol of initial average kinetic energy along the DPT normal mode vector. On a time scale of approximately 100 fs DPT has occurred and the excess energy is redistributed. For the DNA-embedded GC pair considerably more kinetic energy is required (30.0 kcal/mol) for DPT and the process is completed within one hydrogen vibration. The relevance of studies of isolated base pairs and base pair analogs in regard of reactions or properties involving DNA is discussed.

Base Pairing↗
PubMed2004

Investigation of glucose binding sites on insulin.

Possible insulin binding sites for D-glucose have been investigated theoretically by docking and molecular dynamics (MD) simulations. Two different docking programs for small molecules were used; Multiple Copy Simultaneous Search (MCSS) and Solvation Energy for Exhaustive Docking (SEED) programs. The configurations resulting from the MCSS search were evaluated with a scoring function developed to estimate the binding free energy. SEED calculations were performed using various values for the dielectric constant of the solute. It is found that scores emphasizing non-polar interactions gave a preferential binding site in agreement with that inferred from recent fluorescence and NMR NOESY experiments. The calculated binding affinity of -1.4 to -3.5 kcal/mol is within the measured range of -2.0 +/- 0.5 kcal/mol. The validity of the binding site is suggested by the dynamical stability of the bound glucose when examined with MD simulations with explicit solvent. Alternative binding sites were found in the simulations and their relative stabilities were estimated. The motions of the bound glucose during molecular dynamics simulations are correlated with the motions of the insulin side chains that are in contact with it and with larger scale insulin motions. These results raise the question of whether glucose binding to insulin could play a role in its activity. The results establish the complementarity of molecular dynamics simulations and normal mode analyses with the search for binding sites proposed with small molecule docking programs.

Binding Sites↗
PubMed2004

On the influence of semirigid environments on proton transfer along molecular chains.

The dynamics of proton transfer along ammonia chains (chemical composition N(x)H(+)(3x+1), x=2, 4, and 6) in a constraining environment is investigated by ab initio molecular dynamics simulations. A carbon nanotube of defined length and diameter is used as an idealized constraining environment such that the ammonia chain is forced to maintain its quasilinear geometry. It is found that, although the energetics of proton transport shows considerable energetic barriers, proton translocation along the wire is possible at finite temperature for all chain lengths studied. The proton transport involves rotational reorientation of the proton-carrying ammonia molecule. High level ab initio calculations (MP2/aug-cc-pVTZ) yield barriers for internal rotation of 9.1 kcal/mol for NH(4) (+)-NH(3) and 11.7 kcal/mol for OH(3) (+)-OH(2), respectively. The infrared spectrum calculated from the dipole-dipole autocorrelation function shows distinct spectral features in the regions (2000-3000 cm(-1)) where the NHN proton transfer mode is expected to absorb. Assigning moderate opposite total charges between 0.002 and 0.2e to the carbon atoms at the end caps of the nanotube leads to a considerable speedup of the proton transfer.

Journal Article↗