The presence of an aromatic residue (Trp, Phe, Tyr) facing the nonpolar face of galactose is a common feature of galactose-specific lectins. The interactions such as those between the C-H groups of galactose and the pi-electron cloud of aromatic residues have been characterized as weak hydrogen bonds between soft acids and soft bases, largely governed by dispersive and charge transfer interactions. An analysis of the binding sites of several galactose-specific lectins revealed that the spatial position-orientation of galactose relative to the binding site aromatic residue varies substantially. The effect of variations in position-orientations of galactose on the interaction energies of galactose-aromatic residue complexes has not been determined so far. In view of this, MP2/6-311G++** calculations were performed on galactose- and glucose-aromatic residue analogue complexes in eight position-orientations. The results show that the strength of the C-H...pi interactions in galactose-aromatic residue complexes is comparable to that of a hydrogen bond. Rather than the type of aromatic residue, the position-orientation of the saccharide appears to be more critical in determining the strength of their interactions. Earlier studies have found the binding site aromatic residue to be critical, but its role was not clear. This study shows that the aromatic residue is important for discriminating galactose from glucose, in addition to its contribution to binding energy.

An aromatic amino acid is present in the binding site of a number of sugar binding proteins. The interaction of the saccharide with the aromatic residue is determined by their relative position as well as orientation. The position-orientation of the saccharide relative to the aromatic residue was found to vary in different sugar-binding proteins. In the present study, interaction energies of the complexes of galactose (Gal) and of glucose (Glc) with aromatic residue analogs have been calculated by ab initio density functional (U-B3LYP/ 6-31G**) theory. The position-orientations of the saccharide with respect to the aromatic residue observed in various Gal-, Glc-, and mannose-protein complexes were chosen for the interaction energy calculations. The results of these calculations show that galactose can interact with the aromatic residue with similar interaction energies in a number of position-orientations. The interaction energy of Gal-aromatic residue analog complex in position-orientations observed for the bound saccharide in Glc/Man-protein complexes is comparable to the Glc-aromatic residue analog complex in the same position-orientation. In contrast, there is a large variation in interaction energies of complexes of Glc-and of Gal-with the aromatic residue analog in position-orientations observed in Gal-protein complexes. Furthermore, the conformation wherein the O6 atom is away from the aromatic residue is preferred for the exocyclic -CH 2 OH group in Galaromatic residue analog complexes. The implications of these results for saccharide binding in Gal-specific proteins and the possible role of the aromatic amino acid to ensure proper positioning and orientation of galactose in the binding site have been discussed.

No abstract

HIV-1 protease is an essential enzyme in the life cycle of the HIV-1 virus. The conformational dynamics of the flap region of the protease is critical for the ligand binding mechanism, as well as for the catalytic activity. The monoclonal antibody F11.2.32 raised against HIV-1 protease inhibits its activity on binding. We have studied the conformational dynamics of protease in its free, inhibitor ritonavir and antibody bound forms using molecular dynamics simulations. We find that upon Ab binding to the epitope region (residues 36-46) of protease, the overall flexibility of the protease is decreased including the flap region and the active site, which is similar to the decrease in flexibility observed by inhibitor binding to the protease. This suggests an allosteric mechanism to inhibit protease activity. Further, the protease mutants G40E and G40R are known to have decreased activity and were also subjected to MD simulations. We find that the loss of flexibility in the mutants is similar to that observed in the protease bound to the Ab/inhibitor. these insights highlight the role played by dynamics in the function of the protease and how control of flexibility through Ab binding and site specific mutations can inhibit protease activity.The HIV-1 protease belongs to the family of aspartyl protease. It cleaves the newly synthesized polyproteins, which is the vital step to create the mature protein components of an infectious HIV-1 virus 1 . Thus, HIV-1 protease is an essential enzyme in the life-cycle of the HIV-1 virus and a potential target for the structure-based drug design. There are several commercial drugs available in the market against the HIV-1 protease. The success rate of these drugs is low due to the occurrence of drug-resistant mutations in the HIV-1 protease 2 . For example, the thermodynamic integration (TI) and MD simulation studies by Chen and his group suggested that there is a change in the shape and conformation of the binding pocket upon certain drug-resistant mutations and this consequently reduces the binding affinity of inhibitor 3 . Also, conventional method of drug delivery targets active sites and many enzymes with related function may have very similar active sites. This may cause adverse side-effects. Therefore, there is a requirement to generate new generation drugs that can function away from active-site and these are allosteric drugs 4-7 . To identify a new target on HIV-1 protease, there is a need to understand the complete structure and dynamics of HIV-1 protease to inhibits its enzymatic activity. The HIV-1 protease is a homodimeric enzyme with each monomer comprising of 99 amino-acid residues 8 . The active site (residues Asp25, Thr26 and Gly27 from both chains A and B) of the protease is covered by two flaps (residues 43-58) from each chain 9,10 . F11.2.32 is a monoclonal antibody (mAb) raised against the HIV-1 protease. The peptide P36-46 ( 36 MNLPGRWKPKM 46 ) corresponding to the epitope/elbow region of the protease binds to the complementarity determining regions (CDRs) of ...

No abstract

Recently we performed molecular dynamics (MD) simulations on the folding of the hairpin peptide DTVKLMYKGQPMTFR from staphylococcal nuclease in explicit water. We found that the peptide folds into a hairpin conformation with native and nonnative hydrogen-bonding patterns. In all the folding events observed in the folding of the hairpin peptide, loop formation involving the region YKGQP was an important event. In order to trace the origins of the loop propensity of the sequence YKGQP, we performed MD simulations on the sequence starting from extended, polyproline II and native type I' turn conformations for a total simulation length of 300 ns, using the GROMOS96 force field under constant volume and temperature (NVT) conditions. The free-energy landscape of the peptide YKGQP shows minima corresponding to loop conformation with Tyr and Pro side-chain association, turn and extended conformational forms, with modest free-energy barriers separating the minima. To elucidate the role of Gly in facilitating loop formation, we also performed MD simulations of the mutated peptide YKAQP (Gly --> Ala mutation) under similar conditions starting from polyproline II conformation for 100 ns. Two minima corresponding to bend/turn and extended conformations were observed in the free-energy landscape for the peptide YKAQP. The free-energy barrier between the minima in the free-energy landscape of the peptide YKAQP was also modest. Loop conformation is largely sampled by the YKGQP peptide, while extended conformation is largely sampled by the YKAQP peptide. We also explain why the YKGQP sequence samples type II turn conformation in these simulations, whereas the sequence as part of the hairpin peptide DTVKLMYKGQPMTFR samples type I' turn conformation both in the X-ray crystal structure and in our earlier simulations on the folding of the hairpin peptide. We discuss the implications of our results to the folding of the staphylococcal nuclease.

The Gal4p mediated transcriptional activation of GAL genes requires the interaction between Gal3p bound with ATP and galactose and Gal80p. Though numerous studies suggest that galactose and ATP activate Gal3p/Gal1p interaction with Gal80p, neither the mechanism of activation nor the interacting surface that binds to Gal80p is well understood. In this study we investigated the dynamics of Gal3p and Gal1p in the presence and absence of ligands ATP and galactose to understand the role played by dynamics in the function of these proteins through molecular dynamics simulation and protein-protein docking studies. We performed simulations totaling to 510 ns on both Gal1p and Gal3p proteins in the presence and absence of ligands ATP and galactose. We find that, while binding of ligands ATP and galactose to Gal3p/Gal1p do not affect the global conformation of proteins, some local conformational changes around upper-lip helix including insertion domain are observed. We observed that only in the presence of ATP and galactose, Gal3p displays opening and closing motion between the two domains. And because of this motion, a binding interface, which is largely hydrophobic, opens up on the surface of Gal3p and this surface can bind to Gal80p. From our simulation studies we infer probable docking sites for Gal80p on Gal3p/Gal1p, which were further ascertained by the docking of Gal80p on to ligand bound Gal1p and Gal3p proteins, and the residues at the interface between Gal3p and Gal80p are identified. Our results correlate quite well with the existing body of literature on functional and dynamical aspects of Gal1p and Gal3p proteins.

The folding of the sequence (21)DTVKLMYKGQPMTFR(35) from staphylococcal nuclease into a β-hairpin, nucleated by the turn region YKGQP, is known to be an early folding event. With YKGQ being the shortest sequence for a β-turn model and in view of its importance to the folding of staphylococcal nuclease, we investigated the thermodynamics of turn formation at a range of temperatures from 280 to 380 K, with a regular interval of 10 K. Eleven independent molecular dynamics simulations (under NPT conditions) were performed using the GROMACS package of programs and the OPLS-AA/L all-atom force field, each for a time period of 1 μs. Turn formation is supported by enthalpy at lower temperatures, while entropy supports it at higher temperatures. There are modest free energy barriers between turn and extended conformational ensembles. The turn propensity persists even at elevated temperatures. The role of proline in driving the turn formation has been re-examined, and it is inferred that the absence of proline does not affect turn propensity.