Computer Simulation Studies of the Fully Solvated Wild-type and Mutated GnRH in Extended and β-turn Conformations

Abstract. The conformational preference of the gonadotropin-releasing hormone (GnRH) and its Lys-8 mutant, studied earlier with a continuum model, was revisited using an explicit solvent model and thermodynamic integration to calculate the solvent's contribution to the conformationdependence of its free energy. In addition, the Proximity Criterion was used to further analyze the effects of conformational changes

Predicting Protein Interactions by Brownian Dynamics Simulations

We present a newly adapted Brownian-Dynamics (BD)-based protein docking method for predicting native protein complexes. The approach includes global BD conformational sampling, compact complex selection, and local energy minimization. In order to reduce the computational costs for energy evaluations, a shell-based grid force field was developed to represent the receptor protein and solvation effects. The performance of this BD protein docking approach has been evaluated on a test set of 24 crystal protein complexes. Reproduction of experimental structures in the test set indicates the adequate conformational sampling and accurate scoring of this BD protein docking approach. Furthermore, we have developed an approach to account for the flexibility of proteins, which has been successfully applied to reproduce the experimental complex structure from the structure of two unbounded proteins. These results indicate that this adapted BD protein docking approach can be useful for the prediction of protein-protein interactions

Interfacial water as a “hydration fingerprint” in the noncognate complex of BamHI. Biophys

ABSTRACT The molecular code of specific DNA recognition by proteins as a paradigm in molecular biology remains an unsolved puzzle primarily because of the subtle interplay between direct protein-DNA interaction and the indirect contribution from water and ions. Transformation of the nonspecific, low affinity complex to a specific, high affinity complex is accompanied by the release of interfacial water molecules. To provide insight into the conversion from the loose to the tight form, we characterized the structure and energetics of water at the protein-DNA interface of the BamHI complex with a noncognate sequence and in the specific complex. The fully hydrated models were produced with Grand Canonical Monte Carlo simulations. Proximity analysis shows that water distributions exhibit sequence dependent variations in both complexes and, in particular, in the noncognate complex they discriminate between the correct and the star site. Variations in water distributions control the number of water molecules released from a given sequence upon transformation from the loose to the tight complex as well as the local entropy contribution to the binding free energy. We propose that interfacial waters can serve as a ''hydration fingerprint'' of a given DNA sequence

Calculation of the Free Energy and Cooperativity of Protein Folding

Calculation of the free energy of protein folding and delineation of its pre-organization are of foremost importance for understanding, predicting and designing biological macromolecules. Here, we introduce an energy smoothing variant of parallel tempering replica exchange Monte Carlo (REMS) that allows for efficient configurational sampling of flexible solutes under the conditions of molecular hydration. Its usage to calculate the thermal stability of a model globular protein, Trp cage TC5b, achieves excellent agreement with experimental measurements. We find that the stability of TC5b is attained through the coupled formation of local and non-local interactions. Remarkably, many of these structures persist at high temperature, concomitant with the origin of native-like configurations and mesostates in an otherwise macroscopically disordered unfolded state. Graph manifold learning reveals that the conversion of these mesostates to the native state is structurally heterogeneous, and that the cooperativity of their formation is encoded largely by the unfolded state ensemble. In all, these studies establish the extent of thermodynamic and structural pre-organization of folding of this model globular protein, and achieve the calculation of macromolecular stability ab initio, as required for ab initio structure prediction, genome annotation, and drug design

Computer Aided Identification of Small Molecules Disrupting uPAR/α5β1- Integrin Interaction: A New Paradigm for Metastasis Prevention

Disseminated dormant cancer cells can resume growth and eventually form overt metastases, but the underlying molecular mechanism responsible for this change remains obscure. We previously established that cell surface interaction between urokinase receptor (uPAR) and alpha5beta1-integrin initiates a sequel of events, involving MAPK-ERK activation that culminates in progressive cancer growth. We also identified the site on uPAR that binds alpha5beta1-integrin. Disruption of uPAR/integrin interaction blocks ERK activation and forces cancer cells into dormancy.Using a target structure guided computation docking we identified 68 compounds from a diversity library of 13,000 small molecules that were predicted to interact with a previously identified integrin-binding site on uPAR. Of these 68 chemical hits, ten inhibited ERK activation in a cellular assay and of those, 2 compounds, 2-(Pyridin-2-ylamino)-quinolin-8-ol and, 2,2'-(methylimino)di (8-quinolinol) inhibited ERK activation by disrupting the uPAR/integrins interaction. These two compounds, when applied in vivo, inhibited ERK activity and tumor growth and blocked metastases of a model head and neck carcinoma.We showed that interaction between two large proteins (uPAR and alpha5beta1-integrin) can be disrupted by a small molecule leading to profound downstream effects. Because this interaction occurs in cells with high uPAR expression, a property almost exclusive to cancer cells, we expect a new therapy based on these lead compounds to be cancer cell specific and minimally toxic. This treatment, rather than killing disseminated metastatic cells, should induce a protracted state of dormancy and prevent overt metastases