Biomedical subjects

Milan Hodoscek

Publications and source records attributed to Milan Hodoscek.

13 recordsLinked to original sources

PubMed2006

Thrombin inhibitors with novel P1 binding pocket functionality: free energy of binding analysis.

The high incidence of thrombembolic diseases justifies the development of new antithrombotics. The search for a direct inhibitor has resulted in the synthesis of a considerable number of low molecular weight molecules that inhibit human alpha-thrombin potently. However, efforts to develop an orally active drug remain in progress as the most active inhibitors with a highly basic P1 moiety exhibit an unsatisfactory bioavailability profile. In our previous work we solved several X-ray structures of human alpha-thrombin in complexes with (1) novel bicyclic arginine mimetics attached to the glycylproline amide and pyridinone acetamide scaffold and (2) inhibitors with a novel aza scaffold and with charged or neutral P1 moieties. In the present contribution, we correlate the structures of the complex between these inhibitors and the protein with the calculated free energy of binding. The energy of solvation was calculated using the Poisson-Boltzmann approach. In particular, the requirements for successful recognition of an inhibitor at the protein's active site pocket S1 are discussed.

Antithrombins↗

PubMed2006

Quantum mechanical investigation of the atmospheric reaction CH3O2 + NO.

The important stationary points on the potential energy surface of the reaction CH(3)O(2) + NO have been investigated using ab initio and density functional theory techniques. The optimizations were carried out at the B3LYP/6-311++G(d,p) and MP2/6-311++G(d,p) levels of theory while the energetics have been refined using the G2MP2, G3//B3LYP, and CCSD(T) methodologies. The calculations allow the proper characterization of the transition state barriers that determine the fate of the nascent association conformeric minima of methyl peroxynitrite. The main products, CH(3)O + NO(2), are formed through either rearrangement of the trans-conformer to methyl nitrate and its subsequent dissociation or via the breaking of the peroxy bond of the cis-conformer to CH(3)O + NO(2) radical pair. The important consequences of the proposed mechanism are (a) the allowance on energetic grounds for nitrate formation parallel to radical propagation under favorable external conditions and (b) the confirmation of the conformational preference of the homolytic cleavage of the peroxy bond, discussed in previous literature.

Air Pollutants↗

PubMed2006

New force field for calcium binding sites in annexin-membrane complexes.

For accurate classical molecular dynamics (MD) simulations of the calcium mediated bound complexes of annexin and membrane we have developed new force-field parameters correctly describing the interaction of the Ca ion with its environment. We have used quantum chemical calculations to investigate the potential energy surface experienced by the Ca ion within the three different binding sites found in domain 1 of annexin V (ANX V/1). Based on these calculations we were able to quantify the charge polarization of atoms within the binding sites, and to determine the geometry and force constants of harmonic restraints between the Ca ion and its coordinating oxygen atoms. Harmonic restraints were introduced to compensate for the deviations between the quantum mechanical potential energy surface and that of the classical force field. Our analysis has shown that using the refined force field for the Ca binding sites enables long-time MD simulations that conserve the initial structure of ANX V/1 significantly better than MD simulations using the standard force field.

Annexin A5↗

PubMed2005

Electronic spectrum and photodissociation of ClONO in comparison to BrONO.

A theoretical study of the low-lying singlet and triplet states of ClONO is presented. Calculations of excitation energies and oscillator strengths are reported using multireference configuration interaction, MRD-CI, methods with the cc-pVDZ + sp basis set. The calculations predict the dominant transition, 4(1)A' <-- 1(1)A', at 5.70 eV. The transition 2(1)A' <-- 1(1)A', at 4.44 eV, with much lower intensity nicely matches the experimental absorption maximum observed around 290 nm (4.27 eV). The potential energy curves for both states are found to be highly repulsive along the Cl-O coordinate implying that direct and fast dissociation to the Cl + NO2 products will occur. Photodissociation along the N-O coordinate is less likely because of barriers on the order of 0.3 eV for low-lying excited states. A comparison between the calculated electronic energies related to the two dominant excited states of ClONO and BrONO indicates that the transitions lie about 0.6 eV higher if bromine is replaced by chlorine. The stratospheric chemistry implications of ClONO and BrONO are discussed.

Journal Article↗

PubMed2005

Molecular modeling of cobalt(II) hyaluronate.

Structural data for complexes of hyaluronic acid and 3d metals(II) of the fourth group of the periodic table are lacking. A combined QM/MM method was used to solve the structure of the first coordination sphere around the cobalt(II) ion. Some available experimental data were compared with the results obtained via computation and were found to be in good agreement. Our results open the way for using molecular modeling to solve the structure of other metal(II) hyaluronates.

Carbohydrate Conformation↗

PubMed2003

Unfolding of the cold shock protein studied with biased molecular dynamics.

The cold shock protein from Bacillus caldolyticus is a small beta-barrel protein that folds in a two-state mechanism. For the native protein and for several mutants, a wealth of experimental data are available on stability and folding, so that it is an optimal system to study this process. We compare data from unfolding simulations (trajectories of 5 and up to 12 ns) obtained with a bias potential at room temperature and from unbiased thermal unfolding simulations with experimental data. The unfolding patterns derived from the trajectories starting from different native-like conformations and subject to different unfolding conditions agree. The transition state found in the simulations of unfolding is close to the native structure in agreement with experiment. Moreover, a lower value of the free energy barrier of unfolding was found for the mutant R3E than for the mutant E46A and the native protein, as indicated by experimental data. The first unfolding event involves the three-stranded beta-sheet whose decomposition corresponds to the transition state. In contrast to conclusions drawn from experiments, we found that the two-stranded beta-strand forms the most stable substructure, which decomposes very late in the unfolding process. However, assuming that this structure forms very early in the folding process, our findings would not contradict the experiments but require a different interpretation of them.

Bacterial Proteins↗

PubMed2003

Hydrogen bonding interactions between aldose reductase complexed with NADP(H) and inhibitor tolrestat studied by molecular dynamics simulations and binding assay.

Molecular dynamics simulations and binding affinity studies have been carried out in order to probe the effect of the charge state of His110 and cofactor NADPH on the binding affinity of the potent inhibitor tolrestat to aldose reductase (ALR2) complexed with either NADPH or NADP(+). Molecular dynamics simulations of ALR2-NADP(+)-tolrestat indicate that the carboxylate group of tolrestat forms a hydrogen bond with Tyr48 and His110 of ALR2 regardless of the charge state of His110. In the case of ALR2-NADPH-tolrestat, the H-bonding pattern is significantly different from that of ALR2-NADP(+)-tolrestat, in that Tyr48 does not H-bond to tolrestat. The binding affinity of tolrestat to ALR2 complexed with either NADPH or NADP(+) is comparable and pH-dependent. Based on the H-bonding interactions seen in computer simulations, it is proposed that the cationic moiety at the active site of ALR2-NADP(+) and ALR2-NADPH that interacts with the carboxylate of tolrestat is NADP(+) and His110, respectively. The residue that gives rise to the pH-dependent binding of tolrestat to ALR2-NADP(+) and ALR2-NADPH has been identified as Tyr48 and His110, respectively.

Aldehyde Reductase↗

PubMed2002

CROW for large scale macromolecular simulations.

CROW (Columns and Rows Of Workstations - http://www.sicmm.org/crow/) is a parallel computer cluster based on the Beowulf (http://www.beowulf.org/) idea, modified to support a larger number of processors. Its architecture is based on point-to-point network architecture, which does not require the use of any network switching equipment in the system. Thus, the cost is lower, and there is no degradation in network performance even for a larger number of processors.

Computational Biology↗

PubMedDate not supplied

Transition state structure and energetics of the N(2)O + X (X = Cl,Br) reactions.

The structural and vibrational properties of the transition state of the N(2)O + X (X = Cl,Br) reactions have been characterized by ab initio methods using density functional theory. We have employed Becke's hybrid functional (B3LYP), and transition state optimizations were performed with 6-31G(d), 6-311G(2d,2p), 6-311+G(3d,2p), and 6-311+G(3df,2p) basis sets. For the chlorine atom reaction the coupled-cluster method (CCSD(T)) with 6-31G(d) basis set was also used. All calculations resulted in transition state structures with a planar cis arrangement of atoms for both reactions. The geometrical parameters of transition states at B3LYP are very similar, and the reaction coordinates involve mainly the breaking of the N-O bond. At CCSD(T)/6-31G(d) level a contribution of the O-Cl forming bond is also observed in the reaction coordinate. In addition, several highly accurate ab initio composite methods of Gaussian-n (G1, G2, G3), their variations (G2(MP2), G3//B3LYP), and complete basis set (CBS-Q, CBS-Q//B3LYP) series of models were applied to compute reaction energetics. All model chemistries predict exothermic reactions. The G3 and G2 methods result in the smallest deviations from experiment, 1.8 and 0 kcal mol(-1), for the enthalpies of reaction for N(2)O reaction with chlorine and bromine, respectively. The G3//B3LYP and G1 methods perform best among the composite methods in predicting energies of the transition state, with a deviation of 1.9 and 3.0 kcal mol(-1), respectively, in the activation energies for the above processes. However, the B3LYP/6-311+G(3df,2p) method gives smaller deviations of 0.4 and -1.0 kcal mol(-1), respectively. The performance of the methodologies applied in predicting transition state energies was analyzed.

Journal Article↗

PubMedDate not supplied

Enthalpy of the gas-phase CO2 + Mg reaction from ab initio total energies.

Various highly accurate ab initio composite methods of Gaussian-n (G1, G2, G3), their variations (G2(MP2), G3(MP2), G3//B3LYP, G3(MP2)//B3LYP), and complete basis set (CBS-Q, CBS-Q//B3LYP) series of models were applied to compute reaction enthalpies of the ground-state reaction of CO2 with Mg. All model chemistries predict highly endothermic reactions, with DeltaH(298) = 63.6-69.7 kcal x mol(-1). The difference between the calculated reaction enthalpies and the experimental value, evaluated with recommended experimental standard enthalpies of formation for products and reactants, is more than 20 kcal x mol(-1) for all methods. This difference originates in the incorrect experimental enthalpy of formation of gaseous MgO given in thermochemical databases. When the theoretical formation enthalpy for MgO calculated by a particular method is used, the deviation is reduced to 1.3 kcal x mol(-1). The performance of the methodologies used to calculate the heat of this particular reaction and the enthalpy of formation of MgO are discussed.

Journal Article↗

PubMedDate not supplied

Improving the performance of molecular dynamics simulations on parallel clusters.

In this article a procedure is derived to obtain a performance gain for molecular dynamics (MD) simulations on existing parallel clusters. Parallel clusters use a wide array of interconnection technologies to connect multiple processors together, often at different speeds, such as multiple processor computers and networking. It is demonstrated how to configure existing programs for MD simulations to efficiently handle collective communication on parallel clusters with processor interconnections of different speeds.

Journal Article↗

PubMedDate not supplied

Quantum mechanical study of the potential energy surface of the ClO + NO2 reaction.

In the study of the reaction pathways of the ClO + NO2 reaction including reliable structures of the reactants, products, intermediates, and transition states as well as energies the MP2/6-311G(d), B3LYP/6-311G(d), and G2(MP2) methods have been employed. Chlorine nitrate, ClONO2, is formed by N-O association without an entrance barrier and is stabilized by 29.8 kcal mol(-1). It can undergo either a direct 1,3 migration of Cl or OCl rotation to yield an indistinguishable isomer. The corresponding barriers are 45.8 and 7.1 kcal mol(-1), respectively. ClONO2 can further decompose into NO3 + Cl with an endothermicity of 46.4 kcal mol(-1). The overall endothermicity of the NO2 + ClO --> NO3 + Cl reaction is calculated to be 16.6 kcal mol(-1). The formation of cis-perp and trans-perp conformer of chlorine preoxynitrite, ClOONO(cp) and ClOONO(tp), are exothermic by 5.4 and 3.8 kcal mol(-1), respectively. Calculations on the possible reaction pathways for the isomerization of ClOONO to ClONO2 showed that the activation barriers are too high to account for appreciable nitrate formation from peroxynitrite isomerization. All quoted relative energies are related to G2(MP2) calculations.

Journal Article↗

PubMedDate not supplied

Development and validation of empirical force field parameters for netropsin.

The netropsin molecule preferentially binds to the four consecutive A.T base pairs of the DNA minor groove and could therefore inhibit the expression of specific genes. The understanding of its binding on a molecular level is indispensable for computer-aided design of new antitumor agents. This knowledge could be obtained via molecular dynamics (MD) and docking simulations, but in this case appropriate force field parameters for the netropsin molecule should be explicitly defined. Our parametrization was based on the results of quantum chemical calculations. The resulting set of parameters was able to reproduce bond lengths, bond angles, torsional angles of the ab initio minimized geometry within 0.03 A, 3 deg and 5 deg, respectively, and its vibrational frequencies with a relative error of 4.3% for low and 2.8% for high energy modes. To show the accuracy of the developed parameters we calculated an IR spectrum of the netropsin molecule using MD simulation and found it to be in good agreement with the experimental one. Finally, we performed a 10 ns long MD simulation of the netropsin-DNA complex immersed in explicit water. The overall complex conformation remained stable at all times, and its secondary structure was well retained.

Algorithms↗