New-generation amber united-atom force field.

We have developed a new-generation Amber united-atom force field for simulations involving highly demanding conformational sampling such as protein folding and protein-protein binding. In the new united-atom force field, all hydrogens on aliphatic carbons in all amino acids are united with carbons except those on Calpha. Our choice of explicit representation of all protein backbone atoms aims at minimizing perturbation to protein backbone conformational distributions and to simplify development of backbone torsion terms.

Chelate effect in cyclodextrin dimers: a computational (MD, MM/PBSA, and MM/GBSA) study.

The complexation of an adamantyl-phosphate derivative with one beta-cyclodextrin, with two beta-cyclodextrins, and with two beta-cyclodextrins dimerized with a disulfide bridge was studied by computational methods (MD, MM/PBSA, and MM/GBSA) to analyze and rationalize the chelate effect. Although this effect is usually explained by invoking favorable entropy contribution due to the preorganization of the ligand, it has been determined experimentally that in this case it is enthalpy-driven. The computational results are in accord with this finding, although the entropy contribution due to the solvent structural organization around the complex is crucial for the final estimates of the free energy of complexation.

Automatic atom type and bond type perception in molecular mechanical calculations.

In molecular mechanics (MM) studies, atom types and/or bond types of molecules are needed to determine prior to energy calculations. We present here an automatic algorithm of perceiving atom types that are defined in a description table, and an automatic algorithm of assigning bond types just based on atomic connectivity. The algorithms have been implemented in a new module of the AMBER packages.

Hierarchical database screenings for HIV-1 reverse transcriptase using a pharmacophore model, rigid docking, solvation docking, and MM-PB/SA.

In this work, an efficient strategy was presented to search drug leads for human immunodeficiency virus type 1 reverse transcriptase (HIV-1 RT) using hierarchical database screenings, which included a pharmacophore model, multiple-conformation rigid docking, solvation docking, and molecular mechanics-Poisson-Boltzmann/surface area (MM-PB/SA) sequentially. Encouraging results were achieved in searching a refined available chemical directory (ACD) database: the enrichment factor after the first three filters was estimated to be 25-fold; the hit rate for all the four filters was predicted to be 41% in a control test using 37 known HIV-1 non-nucleoside reverse transcriptase inhibitors; 10 out of 30 promising solvation-docking hits had MM-PB/SA binding free energies better than -6.8 kcal/mol and the best one, HIT15, had -17.

Prediction of pKa shifts in proteins using a combination of molecular mechanical and continuum solvent calculations.

The prediction of pKa shifts of ionizable groups in proteins is of great relevance for a number of important biological phenomena. We present an implementation of the MM-GBSA approach, which combines molecular mechanical (MM) and generalized Born (GB) continuum solvent energy terms, to the calculation of pKa values of a panel of nine proteins, including 69 individual comparisons with experiment. While applied so far mainly to the calculation of biomolecular binding free energies, we show that this method can also be used for the estimation of protein pKa shifts, with an accuracy around 1 pKa unit, even for strongly shifted residues.

Development and testing of a general amber force field.

We describe here a general Amber force field (GAFF) for organic molecules. GAFF is designed to be compatible with existing Amber force fields for proteins and nucleic acids, and has parameters for most organic and pharmaceutical molecules that are composed of H, C, N, O, S, P, and halogens. It uses a simple functional form and a limited number of atom types, but incorporates both empirical and heuristic models to estimate force constants and partial atomic charges.

Enhanced ab initio protein folding simulations in Poisson-Boltzmann molecular dynamics with self-guiding forces.

We have investigated the sampling efficiency in molecular dynamics with the PB implicit solvent when self-guiding forces are added. Compared with a high-temperature dynamics simulation, the use of self-guiding forces in room-temperature dynamics is found to be rather efficient as measured by potential energy fluctuation, gyration radius fluctuation, backbone RMSD fluctuation, number of unique clusters, and distribution of low RMSD structures over simulation time. Based on the enhanced sampling method, we have performed ab initio folding simulations of two small proteins, betabetaalpha1 and villin headpiece.

Investigation of neuraminidase-substrate recognition using molecular dynamics and free energy calculations.

Development of the new generation of therapeutics against the influenza viral coat protein neuraminidase is a response to the continuing threat of influenza epidemics. A variety of structurally similar compounds have been reported that vary greatly in their ability to inhibit neuraminidase, a critical enzyme that cleaves sialic acid and promotes virion release. To determine how neuraminidase exhibits this wide range of affinities with structurally similar compounds, molecular dynamic simulations, coupled with free energy calculations, were used to determine the binding components of a series of neuraminidase inhibitors.

Direct hydroxide attack is a plausible mechanism for amidase antibody 43C9.

Direct hydroxide attack on the scissile carbonyl of the substrate has been suggested as a likely mechanism for esterase antibodies elicited by phosphonate haptens, which mimic the transition states for the alkaline hydrolysis of esters.1 The unique amidase activity of esterase antibody 43C9 has been attributed to nucleophilic attack by an active-site histidine residue.2 Yet, the active site of 43C9 is strikingly similar to those of other esterase antibodies, particularly 17E8.

Closing of the nucleotide pocket of kinesin-family motors upon binding to microtubules.

We have used adenosine diphosphate analogs containing electron paramagnetic resonance (EPR) spin moieties and EPR spectroscopy to show that the nucleotide-binding site of kinesin-family motors closes when the motor.diphosphate complex binds to microtubules. Structural analyses demonstrate that a domain movement in the switch 1 region at the nucleotide site, homologous to domain movements in the switch 1 region in the G proteins [heterotrimeric guanine nucleotide-binding proteins], explains the EPR data.

Free energy calculations for theophylline binding to an RNA aptamer: Comparison of MM-PBSA and thermodynamic integration methods.

We have applied the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method (J. Srinivasan, T. E.

The intramolecular mechanism for the second proton transfer in triosephosphate isomerase (TIM): a QM/FE approach.

The intramolecular mechanism we earlier proposed [Alagona, G.; Desmeules, P.; Ghio, C.

pKa, MM, and QM studies of mechanisms of beta-lactamases and penicillin-binding proteins: acylation step.

The acylation step of the catalytic mechanism of beta-lactamases and penicillin-binding proteins (PBPs) has been studied with various approaches. The methods applied range from molecular dynamics (MD) simulations to multiple titration calculations using the Poisson-Boltzmann approach to quantum mechanical (QM) methods. The mechanism of class A beta-lactamases was investigated in the greatest detail.

Molecular dynamics simulation study of the negative correlation in antibody AZ28-catalyzed oxy-cope rearrangement.

The oxy-Cope rearrangement reaction in the antibody AZ28 is investigated using ab initio molecular orbital calculations and molecular mechanical molecular dynamics simulations. This antibody, AZ28, is known as one of the few systems where the mature catalytic antibody shows a negative correlation between the transition state analogue (TSA) binding affinity and the catalytic rate of the oxy-Cope rearrangement compared to the germ line catalytic antibody. The ab initio optimized structure shows that the transition state structure has a more planar configuration than the TSA.

Simulating proteins at constant pH: An approach combining molecular dynamics and Monte Carlo simulation.

For the structure and function of proteins, the pH of the solution is one of the determining parameters. Current molecular dynamics (MD) simulations account for the solution pH only in a limited way by keeping each titratable site in a chosen protonation state. We present an algorithm that generates trajectories at a Boltzmann distributed ensemble of protonation states by a combination of MD and Monte Carlo (MC) simulation.

Conformational preferences of substituted prolines in the collagen triple helix.

Researchers have recently questioned the role hydroxylated prolines play in stabilizing the collagen triple helix. To address these issues, we have developed new molecular mechanics parameters for the simulation of peptides containing 4(R)-fluoroproline (Flp), 4(R)-hydroxyproline (Hyp), and 4(R)-aminoproline (Amp). Simulations of peptides based on these parameters can be used to determine the components that stabilize hydroxyproline over proline in the triple helix.

Computational alanine scanning of the 1:1 human growth hormone-receptor complex.

The MM-PBSA (Molecular Mechanics-Poisson-Boltzmann surface area) method was applied to the human Growth Hormone (hGH) complexed with its receptor to assess both the validity and the limitations of the computational alanine scanning approach. A 400-ps dynamical trajectory of the fully solvated complex was simulated at 300 K in a 101 A x 81 A x 107 A water box using periodic boundary conditions. Long-range electrostatic interactions were treated with the particle mesh Ewald (PME) summation method.

Molecular dynamics and free energy analyses of cathepsin D-inhibitor interactions: insight into structure-based ligand design.

In this study, we compare the calculated and experimental binding free energies for a combinatorial library of inhibitors of cathepsin D (CatD), an aspartyl protease. Using a molecular dynamics (MD)-based, continuum solvent method (MM-PBSA), we are able to reproduce the experimental binding affinity for a set of seven inhibitors with an average error of ca. 1 kcal/mol and a correlation coefficient of 0.

Theoretical and experimental studies of biotin analogues that bind almost as tightly to streptavidin as biotin.

We have used a newly developed qualitative computational approach, PROFEC (Pictorial Representation of Free Energy Changes), to visualize the areas of the ligand biotin where modifications of its structure might lead to tighter binding to the protein streptavidin. The PROFEC analysis, which includes protein flexibility and ligand solvation/desolvation, led to the suggestion that the pro-9R hydrogen atom of biotin, which is in alpha-position to the CO(2)(-) group, might be changed to a larger group and lead to better binding with streptavidin and avidin. Free energy calculations supported this suggestion and predicted that the methyl analogue should bind approximately 3 kcal/mol more tightly to streptavidin, with this difference coming exclusively from the relative desolvation free energy of the ligand.

A classical and ab initio study of the interaction of the myosin triphosphate binding domain with ATP.

We used classical molecular mechanics (MM) simulations and quantum mechanical (QM) structural relaxations to examine the active site of myosin when bound to ATP. Two conformations of myosin have been determined by x-ray crystallography. In one, there is no direct interaction between switch 2 and the nucleotide (open state).

Elucidating the Origin of Conformational Energy Differences in Substituted 1,3-Dioxanes: A Combined Theoretical and Experimental Study.

13C NMR spectroscopy, ab initio quantum mechanics, and molecular mechanics have been used to investigate the trans-4-(trifluoromethyl)-2,2,6-trimethyl-1,3-dioxane chair/twist-boat equilibrium. The molecular mechanics calculations were based upon the MM3 and AMBER force fields. A 6-31G basis set was used for the ab initio calculations, and MP2 correlation corrections were applied.