Enhancing the accuracy of virtual screening: molecular dynamics with quantum-refined force fields.

A methodology aimed at improving the accuracy of current docking-scoring procedures is proposed, and validated through detailed tests of its performance in predicting the activity of HIV-1 protease inhibitors. This methodology is based on molecular dynamics simulations using a force field whose effective charges are refined by means of a novel procedure that relies on quantum-mechanical calculations and preserves the internal consistency of the parameterization scheme.

The binding mode of progesterone to its receptor deduced from molecular dynamics simulations.

An unambiguous understanding of the binding mode of human progesterone to its receptor still eludes experimental search. According to the X-ray structure of the ligand-binding domain, only one (O3) of the two keto groups at the ligand ends (O3 and O20) should play a role. This result is in conflict with chemical intuition and the results of site-directed mutagenesis experiments. Herein, we report classical molecular dynamics simulations that reveal the dynamic nature of the binding in solution, elucidate the reasons why X-ray studies failed to determine the role of O20, and clarify the effects of the mutations. The predictive power of the force field is ensured by the consistent introduction of a first-principles representation of the ligand.

Why do divalent metal ions either promote or inhibit enzymatic reactions? The case of BamHI restriction endonuclease from combined quantum-classical simulations.

Divalent metal ions are essential to many enzymatic reactions involving nucleic acids, but their critical and specific role still needs to be uncovered. Restriction endonucleases are a prominent group of such metal-requiring enzymes. Large scale accurate simulations of Mg- and Ca-BamHI elucidate the mechanism of the catalytic reaction leading to DNA cleavage and show that it involves the concerted action of two metal ions and water molecules. It is also established that what is decisive for the dramatically different behavior of magnesium (a cocatalyst) and calcium (an inhibitor) are kinetic factors and not the properties of the prereactive states of the enzymes. A new perspective is opened for the understanding of the functional role of metal ions in biological processes.

Studies in the assessment of folding quality for protein modeling and structure prediction.

A diagnostic for assessing the quality of a fold has been developed to which further criteria can be progressively added. The goal is to create a measure that can follow the status of a protein structure in a simulation or modeling process, when the answer (the experimental structure) is not known in advance, rather than simply reject deliberate misfolds. This places greater emphasis on the need to study, and calibrate against, marginal cases, i.e., unusual native structures, incomplete structures, partially erroneous X-ray structures, good models, poor models, and the effect of cofactors. The first three terms introduced in the diagnostic are appropriate core-forming properties or noncore properties of residues in relation to tertiary structure, appropriate neighboring structure density for each residue in relation to tertiary structure, and secondary structure consistency. While the method emerges as a useful simulation analysis tool, we find a need for further fine-tuning to diminish sensitivity to minor conformational changes that retain essential features of the fold, balanced against the need to obtain a more sensitive response when a conformational change involves less physically meaningful interatomic interactions. This dual utility is difficult to obtain: the investigation highlights some of the issues. Initial attempts to obtain it have led to terms in the diagnostic that are admittedly complex: simplifications must also be explored.