Predicting relative binding affinities of non-peptide HIV protease inhibitors with free energy perturbation calculations

The relative binding free energies in HIV protease of haloperidol thioketal (THK) and three of its derivatives were examined with free energy calculations. THK is a weak inhibitor (IC50 = 15 M) for which two cocrystal structures with HIV type 1 proteases have been solved [Rutenber, E. et al., J. Biol. Chem., 268 (1993) 15343]. A THK derivative with a phenyl group on C2 of the piperidine ring was expected to be a poor inhibitor based on experiments with haloperidol ketal and its 2- phenyl derivative (Caldera, P., personal communication). Our calculations predict that a 5-phenyl THK derivative, suggested based on examination of the crystal structure, will bind significantly better than THK. Although there are large error bars as estimated from hysteresis, the calculations predict that the 5-phenyl substituent is clearly favored over the 2-phenyl derivative as well as the parent compound. The unfavorable free energies of solvation of both phenyl THK derivatives relative to the parent compound contributed to their predicted binding free energies. In a third simulation, the change in binding free energy for 5-benzyl THK relative to THK was calculated. Although this derivative has a lower free energy in the protein, its decreased free energy of solvation increases the predicted G(bind) to the same range as that of the 2-phenyl derivative.     

Classical Reactive Molecular Dynamics Implementations: State of the Art

Reactive molecular dynamics (RMD) implementations equipped with force field approaches to simulate both the time evolution as well as chemical reactions of a broad class of materials are reviewed herein. We subdivide the RMD approaches developed during the last decade as well as older ones already reviewed in 1995 by Srivastava and Garrison and in 2000 by Brenner into two classes. The methods in the first RMD class rely on the use of a reaction cutoff distance and employ a sudden transition from the educts to the products. Due to their simplicity these methods are well suited to generate equilibrated atomistic or material‐specific coarse‐grained polymer structures. In connection with generic models they offer useful qualitative insight into polymerization reactions. The methods in the second RMD class are based on empirical reactive force fields and implement a smooth and continuous transition from the educts to the products. In this RMD class, the reactive potentials are based on many‐body or bond‐order force fields as well as on empirical standard force fields, such as CHARMM, AMBER or MM3 that are modified to become reactive. The aim with the more sophisticated implementations of the second RMD class is the investigation of the reaction kinetics and mechanisms as well as the evaluation of transition state geometries. Pure or hybrid ab initio, density functional, semi‐empirical, molecular mechanics, and Monte Carlo methods for which no time evolution of the chemical systems is achieved are excluded from the present review. So are molecular dynamics techniques coupled with quantum chemical methods for the treatment of the reactive regions, such as Car–Parinello molecular dynamics.     

Development of a parallel molecular dynamics code on SIMD computers: Algorithm for use of pair list criterion

In recent years several implementations of molecular dynamics (MD) codes have been reported on multiple instruction multiple data (MIMD) machines. However, very few implementations of MD codes on single instruction multiple data (SIMD) machines have been reported. The difficulty in using pair lists of nonbonded interactions is the major problem with MD codes for SIMD machines, such that, generally, the full connectivity computation has been used. We present an algorithm, the global cut-off algorithm (GCA), which permits the use of pair lists on SIMD machines. GCA is based on a probabilistic approach and requires the cut-off condition to be simultaneously verified on all nodes of the machine. The MD code used was taken from the GROMOS package; only the routines involved in the pair lists and in the computation of nonbonded interactions were rewritten for a parallel architecture. The remaining calculations were performed on the host computer. The algorithm has been tested on Quadrics computers for configurations of 32, 128, and 512 processors and for systems of 4000, 8000, 15,000, and 30,000 particles. Quadrics was developed by Istituto Nazionale di Fisica Nucleare (INFN) and marketed by Alenia Spazio. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 685–694, 1998     

Calculation of free energy landscapes: A histogram reweighted metadynamics approach

We present an efficient method for the calculation of free energy landscapes. Our approach involves a history‐dependent bias potential, which is evaluated on a grid. The corresponding free energy landscape is constructed via a histogram reweighting procedure a posteriori. Because of the presence of the bias potential, it can be also used to accelerate rare events. In addition, the calculated free energy landscape is not restricted to the actual choice of collective variables and can in principle be extended to auxiliary variables of interest without further numerical effort. The applicability is shown for several examples. We present numerical results for the alanine dipeptide and the Met‐Enkephalin in explicit solution to illustrate our approach. Furthermore, we derive an empirical formula that allows the prediction of the computational cost for the ordinary metadynamics variant in comparison with our approach, which is validated by a dimensionless representation. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Free Energy Calculations using a Swarm‐Enhanced Sampling Molecular Dynamics Approach

Free energy simulations are an established computational tool in modelling chemical change in the condensed phase. However, sampling of kinetically distinct substates remains a challenge to these approaches. As a route to addressing this, we link the methods of thermodynamic integration (TI) and swarm‐enhanced sampling molecular dynamics (sesMD), where simulation replicas interact cooperatively to aid transitions over energy barriers. We illustrate the approach by using alchemical alkane transformations in solution, comparing them with the multiple independent trajectory TI (IT‐TI) method. Free energy changes for transitions computed by using IT‐TI grew increasingly inaccurate as the intramolecular barrier was heightened. By contrast, swarm‐enhanced sampling TI (sesTI) calculations showed clear improvements in sampling efficiency, leading to more accurate computed free energy differences, even in the case of the highest barrier height. The sesTI approach, therefore, has potential in addressing chemical change in systems where conformations exist in slow exchange.     

Excited‐State Deactivation of Branched Phthalocyanine Compounds

The excited‐state relaxation dynamics and chromophore interactions in two phthalocyanine compounds (bis‐ and trisphthalocyanines) are studied by using steady‐state and femtosecond transient absorption spectral measurements, where the excited‐state energy‐transfer mechanism is explored. By exciting phthalocyanine compounds to their second electronically excited states and probing the subsequent relaxation dynamics, a multitude of deactivation pathways are identified. The transient absorption spectra show the relaxation pathway from the exciton state to excimer state and then back to the ground state in bisphthalocyanine (bis‐Pc). In trisphthalocyanine (tris‐Pc), the monomeric and dimeric subunits are excited and the excitation energy transfers from the monomeric vibrationally hot S1 state to the exciton state of a pre‐associated dimer, with subsequent relaxation to the ground state through the excimer state. The theoretical calculations and steady‐state spectra also show a face‐to‐face conformation in bis‐Pc, whereas in tris‐Pc, two of the three phthalocyanine branches form a pre‐associated face‐to‐face dimeric conformation with the third one acting as a monomeric unit; this is consistent with the results of the transient absorption experiments from the perspective of molecular structure. The detailed structure–property relationships in phthalocyanine compounds is useful for exploring the function of molecular aggregates in energy migration of natural photosynthesis systems.