Role of conformational fluctuations in the enzymatic reaction of HIV-1 protease

The folding free-energy surface of HIV-1 protease: insights into the thermodynamic basis for resistance to inhibitors

Spontaneous mutations at numerous sites distant from the active site of human immunodeficiency virus type 1 protease enable resistance to inhibitors while retaining enzymatic activity. As a benchmark for probing the effects of these mutations on the conformational adaptability of this dimeric beta-barrel protein, the folding free-energy surface of a pseudo-wild-type variant, HIV-PR(*), was determined by a combination of equilibrium and kinetic experiments on the urea-induced unfolding/refolding reactions. The equilibrium unfolding reaction was well described by a two-state model involving only the native dimeric form and the unfolded monomer. The global analysis of the kinetic folding mechanism reveals the presence of a fully folded monomeric intermediate that associates to form the native dimeric structure. Independent analysis of a stable monomeric version of the protease demonstrated that a small-amplitude fluorescence phase in refolding and unfolding, not included in the global analysis of the dimeric protein, reflects the presence of a transient intermediate in the monomer folding reaction. The partially folded and fully folded monomers are only marginally stable with respect to the unfolded state, and the dimerization reaction provides a modest driving force at micromolar concentrations of protein. The thermodynamic properties of this system are such that mutations can readily shift the equilibrium from the dimeric native state towards weakly folded states that have a lower affinity for inhibitors but that could be induced to bind to their target proteolytic sites. Presumably, subsequent secondary mutations increase the stability of the native dimeric state in these variants and, thereby, optimize the catalytic properties of the resistant human immunodeficiency virus type 1 protease.

Peptidomimetic therapeutic agents targeting the protease enzyme of the human immunodeficiency virus and hepatitis C virus

During the past two decades, great strides have been made in the design of peptidomimetic drugs for the treatment of viral infections, despite the stigma of poor drug-like properties, low oral absorption, and high clearance associated with such compounds. This Account summarizes the progress made toward overcoming such liabilities and highlights the drug discovery efforts that have focused specifically on human immunodeficiency virus (HIV) and hepatitis C virus (HCV) protease inhibitors. The arsenal against the incurable disease AIDS, which is caused by HIV infection, includes peptidomimetic compounds that target the virally encoded aspartic protease enzyme. This enzyme is essential to the production of mature HIV particles and plays a key role in maintaining infectivity. However, because of the rapid genomic evolution of viruses, an inevitable consequence in the treatment of all viral infections is the emergence of resistance to the drugs. Therefore, the incomplete suppression of HIV in treatment-experienced AIDS patients will continue to drive the search for more effective therapeutic agents that exhibit efficacy against the mutants raised by the earlier generation of protease inhibitors. Currently, a number of substrate-based peptidomimetic agents that target the virally encoded HCV NS3/4A protease are in clinical development. Mechanistically, these inhibitors can be generally divided into activated carbonyls that are transition-state mimics or compounds that tap into the feedback mode of enzyme-product inhibition. In the HCV field, there is justified optimism that a number of these compounds will soon reach commercialization as therapeutic agents for the treatment of HCV infections. Structural research has guided the successful design of both HIV and HCV protease inhibitors. X-ray crystallography, NMR, and computational studies have provided valuable insight in to the free-state preorganization of peptidomimetic ligands and their enzyme-bound conformation. Researchers have designed a variety of novel bioisosteric replacements of amino acids and short peptides that contain all of the required pharmacophore moieties and play a key role in inducing conformational changes to the overall molecule. The knowledge gained from these studies will undoubtedly guide the future design of therapeutic agents and further contribute to the success of this field.

A poke in the eye: inhibiting HIV-1 protease through its flap-recognition pocket

A novel mechanism of inhibiting HIV-1 protease (HIVp) is presented. Using computational solvent mapping to identify complementary interactions and the Multiple Protein Structure method to incorporate protein flexibility, we generated a receptor-based pharmacophore model of the flexible flap region of the semiopen, apo state of HIVp. Complementary interactions were consistently observed at the base of the flap, only within a cleft with a specific structural role. In the closed, bound state of HIVp, each flap tip docks against the opposite monomer, occupying this cleft. This flap-recognition site is filled by the protein and cannot be identified using traditional approaches based on bound, closed structures. Virtual screening and dynamics simulations show how small molecules can be identified to complement this cleft. Subsequent experimental testing confirms inhibitory activity of this new class of inhibitor. This may be the first new inhibitor class for HIVp since dimerization inhibitors were introduced 17 years ago.

Differential native state ruggedness of the two Ca2+-binding domains in a Ca2+ sensor protein

Characterization of near-native excited states of a protein provides insights into various biological functions such as co-operativity, protein-ligand, and protein-protein interactions. In the present study, we investigated the ruggedness of the native state of EhCaBP using nonlinear temperature dependence of backbone amide-proton chemical shifts. EhCaBP is a two-domain EF-hand calcium sensor protein consisting of two EF-hands in each domain and binds four Ca2+ ions. It has been observed that approximately 30% of the residues in the protein access alternative conformations. Theoretical modeling suggested that these low-energy excited states are within 2-3 kcal/mol from the native state. Further, it is interesting to note that the residues accessing alternative conformations are more dominated in the C-terminal domain compared with its N-terminal counterpart suggesting that the former is more rugged in its native state. These distinct characteristics of N- and C-terminal domains of a calcium sensor protein belonging to the super family of calmodulin would have implications for domain dependent Ca2+ signaling pathways.

Conformation of a coarse-grained protein chain (an aspartic acid protease) model in effective solvent by a bond-fluctuating Monte Carlo simulation

In a coarse-grained description of a protein chain, all of the 20 amino acid residues can be broadly divided into three groups: Hydrophobic (H) , polar (P) , and electrostatic (E) . A protein can be described by nodes tethered in a chain with a node representing an amino acid group. Aspartic acid protease consists of 99 residues in a well-defined sequence of H , P , and E nodes tethered together by fluctuating bonds. The protein chain is placed on a cubic lattice where empty lattice sites constitute an effective solvent medium. The amino groups (nodes) interact with the solvent (S) sites with appropriate attractive (PS) and repulsive (HS) interactions with the solvent and execute their stochastic movement with the Metropolis algorithm. Variations of the root mean square displacements of the center of mass and that of its center node of the protease chain and its gyration radius with the time steps are examined for different solvent strength. The structure of the protease swells on increasing the solvent interaction strength which tends to enhance the relaxation time to reach the diffusive behavior of the chain. Equilibrium radius of gyration increases linearly on increasing the solvent strength: A slow rate of increase in weak solvent regime is followed by a faster swelling in stronger solvent. Variation of the gyration radius with the time steps suggests that the protein chain moves via contraction and expansion in a somewhat quasiperiodic pattern particularly in strong solvent.

Second step of hydrolytic dehalogenation in haloalkane dehalogenase investigated by QM/MM methods

Mechanistic studies on the hydrolytic dehalogenation catalyzed by haloalkane dehalogenases are of importance for environmental and industrial applications. Here, Car-Parrinello (CP) and ONIOM hybrid quantum-mechanical/molecular mechanics (QM/MM) are used investigate the second reaction step of the catalytic cycle, which comprises a general base-catalyzed hydrolysis of an ester intermediate (EI) to alcohol and free enzyme. We focus on the enzyme LinB from Sphingomonas paucimobilis UT26, for which the X-ray structure at atomic resolution is available. In agreement with previous proposals, our calculations suggest that a histidine residue (His272), polarized by glutamate (Glu132), acts as a base, accepting a proton from the catalytic water molecule and transferring it to an alcoholate ion. The reaction proceeds through a metastable tetrahedral intermediate, which shows an easily reversed reaction to the EI. In the formation of the products, the protonated aspartic acid (Asp108) can easily adopt conformation of the relaxed state found in the free enzyme. The overall free energy barrier of the reaction calculated by potential of the mean force integration using CP-QM/MM calculations is equal to 19.5 +/- 2 kcal . mol(-1). The lowering of the energy barrier of catalyzed reaction with respect to the water reaction is caused by strong stabilization of the reaction intermediate and transition state and their preorganization by electrostatic field of the enzyme.

Microseconds dynamics simulations of the outer-membrane protease T

Conformational fluctuations of enzymes may play an important role for substrate recognition and/or catalysis, as it has been suggested in the case of the protease enzymatic superfamily. Unfortunately, theoretically addressing this issue is a problem of formidable complexity, as the number of the involved degrees of freedom is enormous: indeed, the biological function of a protein depends, in principle, on all its atoms and on the surrounding water molecules. Here we investigated a membrane protease enzyme, the OmpT from Escherichia coli, by a hybrid molecular mechanics/coarse-grained approach, in which the active site is treated with the GROMOS force field, whereas the protein scaffold is described with a Go-model. The method has been previously tested against results obtained with all-atom simulations. Our results show that the large-scale motions and fluctuations of the electric field in the microsecond timescale may impact on the biological function and suggest that OmpT employs the same catalytic strategy as aspartic proteases. Such a conclusion cannot be drawn within the 10- to 100-ns timescale typical of current molecular dynamics simulations. In addition, our studies provide a structural explanation for the drop in the catalytic activity of two known mutants (S99A and H212A), suggesting that the coarse-grained approach is a fast and reliable tool for providing structure/function relationships for both wild-type OmpT and mutants.

A molecular dynamics study comparing a wild-type with a multiple drug resistant HIV protease: Differences in flap and aspartate 25 cavity dimensions

HIV proteases can develop resistance to therapeutic drugs by mutating specific residues, but still maintain activity with their natural substrates. To gain insight into why mutations confer such resistance, long ( approximately 70 ns) Molecular Dynamics simulations in explicit solvent were performed on a multiple drug resistant (MDR) mutant (with Asn25 in the crystal structure mutated in silico back to the catalytically active Asp25) and a wild type (WT) protease. HIV proteases are homodimers, with characteristic flap tips whose conformations and dynamics are known to be important influences of ligand binding to the aspartates that form the catalytic center. The WT protease undergoes a transition between 25 and 35 ns that is absent in the MDR protease. The origin of this distinction is investigated using principal component analysis, and is related to differences in motion mainly in the flap region of each monomer. Trajectory analysis suggests that the WT transition arises from a concerted motion of the flap tip distances to their catalytic aspartate residues, and the distance between the two flap tips. These distances form a triangle that in the WT expands the active site from an initial (semi-open) form to an open form, in a correlated manner. In contrast, the MDR protease remains in a more closed configuration, with uncorrelated fluctuations in the distances defining the triangle. This contrasting behavior suggests that the MDR mutant achieves its resistance to drugs by making its active site less accessible to inhibitors. The migration of water to the active site aspartates is monitored. Water molecules move in and out of the active site and individual waters hydrogen bond to both aspartate carboxylate oxygens, with residence times in the ns time regime.

Insight into catalytically relevant correlated motions in human purine nucleoside phosphorylase

The catalytic site of the homotrimeric enzyme human purine nucleoside phosphorylase enzyme (hPNP) features residue F200 and the 241-265 loop directly skirting the purine base and a residue belonging to the adjacent monomer, F159, immediately conterminous to the ribosyl moiety. Crystallographic B-factors of apo human purine nucleoside phosphorylase, and hPNP complexed with substrate/transition state (TS) analogues, show that residue E250 is the centroid of a highly mobile loop region. Furthermore, superimposition of apo hPNP and hPNP complexed with TS analogue Immucillin-H shows a tightening of the active site, caused by the ligand-dependent 241-265 loop rearrangement taking place upon substrate/inhibitor binding, suggesting a putative dynamic role of the loop in binding/catalysis. However, crystallographic structures reveal only average atomic positions, and more detailed information is needed to discern the dynamic behavior of hPNP. The Essential Dynamics (ED) method is used here to investigate the existence of correlated motions in hPNP and consequently proposes mutagenesis assays to estimate the relative importance of these motions in the phosphorolytic efficiency of the reaction catalyzed by hPNP. We compare the concerted motions obtained from multiple molecular dynamics simulations of apo and Michaelis complex of hPNP both in vacuo and in solution. The results of the principal component analysis for the apo hPNP indicate the existence of strong correlations predominantly in the vicinity of residue F159. However, for the Michaelis complex, concerted motions are seen mostly around both active site residue F200 and loop residue E250. Additionally, for a simulation depicting the relaxation of tight complexed hPNP with a TS analogue, toward its relaxed apo form (after removal of the TS analog), a combination of the apo hPNP and Michaelis complex motions is found, with prominent concerted modes centered around neighboring residues F159, F200, and E250. Finally, we probed the extent to which these concerted motions bear an intrinsic catalytic role by performing experimental site-directed mutagenesis on some residues, followed by kinetic analysis. The F159G and F200G mutants displayed a strong increase in K(M) and modest decrease in k(cat), suggesting that these concerted motions may provide dynamical roles in substrate binding and/or catalysis. However, further structural data for the hPNP mutants are needed to confirm our hypothesis.

Resistant mechanism against nelfinavir of human immunodeficiency virus type 1 proteases

Inhibitors against human immunodeficiency virus type-1 (HIV-1) proteases are finely effective for anti-HIV-1 treatments. However, the therapeutic efficacy is reduced by the rapid emergence of inhibitor-resistant variants of the protease. Among patients who failed in the inhibitor nelfinavir (NFV) treatment, D30N, N88D, and L90M mutations of HIV-1 protease are often observed. Despite the serious clinical problem, it is not clear how these mutations, especially nonactive site mutations N88D and L90M, affect the affinity of NFV or why they cause the resistance to NFV. In this study, we executed molecular dynamics simulations of the NFV-bound proteases in the wild-type and D30N, N88D, D30N/N88D, and L90M mutants. Our simulations clarified the conformational change at the active site of the protease and the change of the affinity with NFV for all of these mutations, even though the 88th and 90th residues are not located in the NFV-bound cavity and not able to directly interact with NFV. D30N mutation causes the disappearance of the hydrogen bond between the m-phenol group of NFV and the 30th residue. N88D mutation alters the active site conformation slightly and induces a favorable hydrophobic contact. L90M mutation dramatically changes the conformation at the flap region and leads to an unfavorable distortion of the binding pocket of the protease, although 90M is largely far apart from the flap region. Furthermore, the changes of binding energies of the mutants from the wild-type protease are shown to be correlated with the mutant resistivity previously reported by the phenotypic experiments.

Atomic resolution crystal structures of HIV-1 protease and mutants V82A and I84V with saquinavir

Saquinavir (SQV), the first antiviral HIV-1 protease (PR) inhibitor approved for AIDS therapy, has been studied in complexes with PR and the variants PR(I) (84V) and PR(V) (82A) containing the single mutations I84V and V82A that provide resistance to all the clinical inhibitors. Atomic resolution crystal structures (0.97-1.25 A) of the SQV complexes were analyzed in comparison to the protease complexes with darunavir, a new drug that targets resistant HIV, in order to understand the molecular basis of drug resistance. PR(I) (84V) and PR(V) (82A) complexes were obtained in both the space groups P2(1)2(1)2 and P2(1)2(1)2(1), which provided experimental limits for the conformational flexibility. The SQV interactions with PR were very similar in the mutant complexes, consistent with the similar inhibition constants. The mutation from bigger to smaller amino acids allows more space to accommodate the large group at P1' of SQV, unlike the reduced interactions observed in darunavir complexes. The residues 79-82 have adjusted to accommodate the large hydrophobic groups of SQV, suggesting that these residues are intrinsically flexible and their conformation depends more on the nature of the inhibitor than on the mutations in this region. This analysis will assist with development of more effective antiviral inhibitors.

HIV-1 protease substrate binding and product release pathways explored with coarse-grained molecular dynamics

We analyze the encounter of a peptide substrate with the native HIV-1 protease, the mechanism of substrate incorporation in the binding cleft, and the dissociation of products after substrate hydrolysis. To account for the substrate, we extend a coarse-grained model force field, which we previously developed to study the flap opening dynamics of HIV-1 protease on a microsecond timescale. Molecular and Langevin dynamics simulations show that the flaps need to open for the peptide to bind and that the protease interaction with the substrate influences the flap opening frequency and interval. On the other hand, release of the products does not require flap opening because they can slide out from the binding cleft to the sides of the enzyme. Our data show that in the protease-substrate complex the highest fluctuations correspond to the 17- and 39-turns and the substrate motion is anticorrelated with the 39-turn. Moreover, the active site residues and the flap tips move in phase with the peptide. We suggest some mechanistic principles for how the flexibility of the protein may be involved in ligand binding and release.

Targeting structural flexibility in HIV-1 protease inhibitor binding

HIV-1 protease remains an important anti-AIDS drug target. Although it has been known that ligand binding induces large conformational changes in the protease, the dynamic aspects of binding have been largely ignored. Several computational models describing protease dynamics have been reported recently. These have reproduced experimental observations, and have also explained how ligands gain access to the binding site through dynamic behavior of the protease. Specifically, the transitions between three different conformations of the protein have been modeled in atomic detail. Two of these forms were determined by crystallography, and the third was implied by NMR experiments. Based on these computational models, it has been suggested that binding of inhibitors in allosteric sites might affect protease flexibility and disrupt its function.

In-situ Synthesis of a Tacrine-Triazole-Based Inhibitor of Acetylcholinesterase: Configurational Selection Imposed by Steric Interactions

Recently, researchers have used acetylcholinesterase (AChE) as a reaction vessel to synthesize its own inhibitors. Thus, 1 (syn-TZ2PA6), a femtomolar AChE inhibitor, which is formed in a 1:1 mixture with its anti-isomer by solution phase reaction from 3 (TZ2) and 4 (PA6), can be synthesized exclusively inside the AChE gorge. Our computational approach based on quantum mechanical/molecular mechanical (QM/MM) calculations, molecular dynamics (MD), and targeted molecular dynamics (TMD) studies answers why 1 is the sole product in the AChE environment. Ab initio QM/MM results show that the reaction in the AChE gorge occurs when 3/azide and 4/acetylene are extended in a parallel orientation. An MD simulation started from the final structure of QM/MM calculations keeps the azide's and acetylene's parallel orientations intact for 10 ns of simulation time. A TMD simulation applied on an antiparallel azide-acetylene conformation flips the acetylene easily to bring it to a position that is parallel to azide. A second set of QM/MM calculations performed on this flipped structure generates a similar minimum-energy path as obtained previously. Even a TMD simulation carried out on a parallel azide-acetylene conformation could not deform their parallel arrangement. All of these results, thus, imply that inside the AChE gorge, the azide group of 3 and the acetylene group of 4 always remain parallel, with the consequence that 1 is the only product. The architecture of the gorge plays an important role in this selective formation of 1.

Anisotropic network model: systematic evaluation and a new web interface

MOTIVATION

The Anisotropic Network Model (ANM) is a simple yet powerful model for normal mode analysis of proteins. Despite its broad use for exploring biomolecular collective motions, ANM has not been systematically evaluated to date. A lack of a convenient interface has been an additional obstacle for easy usage.

RESULTS

ANM has been evaluated on a large set of proteins to establish the optimal model parameters that achieve the highest correlation with experimental data and its limits of accuracy and applicability. Residue fluctuations in globular proteins are shown to be more accurately predicted than those in nonglobular proteins, and core residues are more accurately described than solvent-exposed ones. Significant improvement in agreement with experiments is observed with increase in the resolution of the examined structure. A new server for ANM calculations is presented, which offers flexible options for controlling model parameters and output formats, interactive animation of collective modes and advanced graphical features.

AVAILABILITY

ANM server (http://www.ccbb.pitt.edu/anm)

Restrained molecular dynamics simulations of HIV-1 protease: the first step in validating a new target for drug design

To test the anticorrelated relationship that was recently displayed in conventional molecular dynamics (MD) simulations, several different restrained MD simulations on a wild type and on the V82F/I84V drug-resistant mutant of HIV-1 protease were performed. This anticorrelated relationship refers to the observation that compression of the peripheral ear-to-cheek region of HIV protease (i.e., the elbow of the flap to the fulcrum and the cantilever) occurred as the active site flaps were opening, and, conversely, expansion of that ear-to-cheek region occurred as both flaps were closing. An additional examination of this anticorrelated relationship was necessary to determine whether it can be harnessed in a useful manner. Consequently, six different MD experiments were performed that incorporated pairwise distance restraints in that ear-to-cheek region (i.e., the distance between the alpha-carbons of Gly40 and Gln61 was restrained to either 7.7 or 10.5 A, in both monomers). Pushing the backbones of the ear and the cheek regions away from each other slightly did force the flaps that guard the active site to remain closed in both the wild type and the mutant systems-even though there were no ligands in the active sites. Thus, these restrained MD simulations provided evidence that the anticorrelated relationship can be exploited to affect the dynamic behavior of the flaps that guard the active site of HIV-1 protease. These simulations supported our hypothesis of the mechanism governing flap motion, and they are the first step towards validating that peripheral surface as a new target for drug design.

Convergent Dynamics in the Protease Enzymatic Superfamily

Proteases regulate various aspects of the life cycle in all organisms by cleaving specific peptide bonds. Their action is so central for biochemical processes that at least 2% of any known genome encodes for proteolytic enzymes. Here we show that selected proteases pairs, despite differences in oligomeric state, catalytic residues, and fold, share a common structural organization of functionally relevant regions which are further shown to undergo similar concerted movements. The structural and dynamical similarities found pervasively across evolutionarily distant clans point to common mechanisms for peptide hydrolysis.

An unusual functional group interaction and its potential to reproduce steric and electrostatic features of the transition states of peptidolysis

The donor-acceptor interaction between a tertiary amine and an aldehyde, first observed among a select class of alkaloids, was deliberately established in a peptidomimetic (1a-c) to mimic features of the two principal transition states of peptide hydrolysis. Compounds 1a-c show preferential adoption in methanol and water of a 'folded' conformation displaying the interaction. Proportions of the folded form in MeOH range from 45% to 70% and can reach 84% in buffer. Significantly, three tendencies for the folded/unfolded equilibrium are observed: increasing solubility and polarity of the medium and decreasing temperature results in a higher extent of folding. In the absence of any parameter set available for this weak bond, no modeling studies were conducted to aid in the design of 1a-c. The successful straightforward synthesis of 1 and its folding and inhibition results with HIV-1 peptidase using FRET technology encourage studies to further pre-organize candidate molecules and to screen the structure space by modeling and parallel combinatorial chemistry.

Drug Resistance of HIV-1 Protease Against JE-2147: I47V Mutation Investigated by Molecular Dynamics Simulation

Anti-retroviral therapies for acquired immunodeficiency syndrome (AIDS) patients are at risk because of drug resistance that has been identified with a number of currently marketed drugs. HIV-1 protease (HIV-pr), a well-validated AIDS therapeutic target, undergoes various mutations leading to resistance such existing drugs. However, the molecular basis of drug resistance of HIV-pr is still not fully understood. JE-2147, an experimental inhibitor of HIV-pr, shows a resistance profile different from that of known drugs. Noteworthy, it is less susceptible to several common mutations, but it is still susceptible to a few mutations, including I47V which appears to be specific for JE-2147. In this work, the molecular details of the effect of I47V mutation is investigated using molecular dynamics simulation. Four simulations of apo and complexed proteins in their wild type (WT) and mutant forms have been performed. It is found that the mobility of the side chain of mutant Val47 in chain B of HIV-pr about the inhibitor increases significantly relative to WT Ile47 in chain B. This is due to loss of optimized packing of the inhibitor to the residue 47 in chain B of the mutant when compared with WT enzyme. There also are subtle differences in motion involving residues in the flap region, which are more prominent in the apo form.

Molecular basis for substrate recognition and drug resistance from 1.1 to 1.6 angstroms resolution crystal structures of HIV-1 protease mutants with substrate analogs

HIV-1 protease (PR) and two drug-resistant variants--PR with the V82A mutation (PR(V82A)) and PR with the I84V mutation (PR(I84V))--were studied using reduced peptide analogs of five natural cleavage sites (CA-p2, p2-NC, p6pol-PR, p1-p6 and NC-p1) to understand the structural and kinetic changes. The common drug-resistant mutations V82A and I84V alter residues forming the substrate-binding site. Eight crystal structures were refined at resolutions of 1.10-1.60 A. Differences in the PR-analog interactions depended on the peptide sequence and were consistent with the relative inhibition. Analog p6(pol)-PR formed more hydrogen bonds of P2 Asn with PR and fewer van der Waals contacts at P1' Pro compared with those formed by CA-p2 or p2-NC in PR complexes. The P3 Gly in p1-p6 provided fewer van der Waals contacts and hydrogen bonds at P2-P3 and more water-mediated interactions. PR(I84V) showed reduced van der Waals interactions with inhibitor compared with PR, which was consistent with kinetic data. The structures suggest that the binding affinity for mutants is modulated by the conformational flexibility of the substrate analogs. The complexes of PR(V82A) showed smaller shifts of the main chain atoms of Ala82 relative to PR, but more movement of the peptide analog, compared to complexes with clinical inhibitors. PR(V82A) was able to compensate for the loss of interaction with inhibitor caused by mutation, in agreement with kinetic data, but substrate analogs have more flexibility than the drugs to accommodate the structural changes caused by mutation. Hence, these structures help to explain how HIV can develop drug resistance while retaining the ability of PR to hydrolyze natural substrates.

Coarse-grained model of proteins incorporating atomistic detail of the active site

We present a novel approach to explore the conformational space of globular proteins near their native state. It combines the advantages of coarse-grained models with those of all-atoms simulations, required to treat molecular recognition processes. The comparison between calculated structural properties with those obtained with all-atoms molecular dynamics simulations establishes the accuracy of the model. Our method has the potential to be extended to molecular recognition processes in systems whose characteristic size and time scale prevent an analysis based on all-atoms molecular dynamics.

Role of phosphorylated Thr160 for the activation of the CDK2/Cyclin A complex

The enzymatic activity of the CDK2/Cyclin A complex increases upon the specific phosphorylation of Thr160@CDK2. In the present study, we have performed a comparative molecular dynamics (MD) study of models of the complex CDK2/Cyclin A/Substrate, which differ for the presence or absence of the phosphate group bound to Thr160. The models are based on two X-ray structures available for CDK2/CyclinA and pCDK2/CyclinA/Substrate complexes. In this way, we analyze the influence of the phosphorylated Thr160 (pThr160) on both the flexibility of CDK2 activation loop (AL) and substrate binding in CDK2. Our calculations point to a decreased flexibility of the AL in the phosphorylated model, in fairly good agreement with experimental data, and to a key role of pThr160 for substrate recognition and stability. Multiple alignments of the CDKs sequences point to the very high conservation of the AL sequence among the CDKs, thus extending our results to all CDKs.

Is allostery an intrinsic property of all dynamic proteins?

Allostery involves coupling of conformational changes between two widely separated binding sites. The common view holds that allosteric proteins are symmetric oligomers, with each subunit existing in "at least" two conformational states with a different affinity for ligands. Recent observations such as the allosteric behavior of myoglobin, a classical example of a nonallosteric protein, call into question the existing allosteric dogma. Here we argue that all (nonfibrous) proteins are potentially allosteric. Allostery is a consequence of re-distributions of protein conformational ensembles. In a nonallosteric protein, the binding site shape may not show a concerted second-site change and enzyme kinetics may not reflect an allosteric transition. Nevertheless, appropriate ligands, point mutations, or external conditions may facilitate a population shift, leading a presumably nonallosteric protein to behave allosterically. In principle, practically any potential drug binding to the protein surface can alter the conformational redistribution. The question is its effectiveness in the redistribution of the ensemble, affecting the protein binding sites and its function. Here, we review experimental observations validating this view of protein allostery.

Disordered p27Kip1 exhibits intrinsic structure resembling the Cdk2/cyclin A-bound conformation

p27Kip1 (p27) influences cell division by regulating nuclear cyclin-dependent kinases. Before binding, p27 is at least partially disordered and folds upon binding its Cdk/cyclin targets. 30-40% of human proteins, including p27, are predicted to contain disordered segments, and have been termed intrinsically unstructured proteins (IUPs). Unfortunately, the inherent dynamics of IUPs hamper detailed analysis of their structure/function relationships. Here, we describe the use of molecular dynamics (MD) computations and solution NMR spectroscopy to reveal that several segments of the p27 kinase inhibitory domain (p27-KID), in addition to the previously characterized helical segment, exist as highly populated, intrinsically folded structural units (IFSUs). Several IFSUs resemble structural features of bound p27-KID, while another exhibits alternative conformations. Interestingly, the highly conserved, specificity determining segment of p27 is shown to be highly disordered. Elucidation of IFSUs within p27-KID allows consideration of their influences on the thermodynamics and kinetics of Cdk/cyclin binding. The degree to which IFSUs are populated within p27-KID is surprising and suggests that other putative IUPs contain IFSUs that may be studied using similar techniques.

Mimicking natural evolution in metallo-beta-lactamases through second-shell ligand mutations

Metallo-beta-lactamases (MBLs) represent the latest generation of beta-lactamases. The structural diversity and broad substrate profile of MBLs allow them to confer resistance to most beta-lactam antibiotics. To explore the evolutionary potential of these enzymes, we have subjected the Bacillus cereus MBL (BcII) to a directed evolution scheme, which resulted in an increased hydrolytic efficiency toward cephalexin. A systematic study of the hydrolytic profile, substrate binding, and active-site features of the evolved lactamase reveal that directed evolution has shaped the active site by means of remote mutations to better hydrolyze cephalosporins with small, uncharged C-3 substituents. One of these mutations is found in related enzymes from pathogenic bacteria and is responsible for the increase in that enzyme's hydrolytic profile. The mutations lowered the activation energy of the rate-limiting step rather than improved the affinity of the enzyme toward these substrates. The following conclusions can be made: (i) MBLs are able to expand their substrate spectrum without sacrificing their inherent hydrolytic capabilities; (ii) directed evolution is able to mimic mutations that occur in nature; (iii) the metal-ligand strength is tuned by second-shell mutations, thereby influencing the catalytic efficiency; and (iv) changes in the position of the second Zn(II) ion in MBLs affect the substrate positioning in the active site. Overall, these results show that the evolution of enzymatic catalysis can take place by remote mutations controlling reactivity.

Evolutionarily conserved functional mechanics across pepsin-like and retroviral aspartic proteases

The biological function of the aspartic protease from HIV-1 has recently been related to the conformational flexibility of its structural scaffold. Here, we use a multistep strategy to investigate whether the same mechanism affects the functionality in the pepsin-like fold. (i) We identify the set of conserved residues by using sequence-alignment techniques. These residues cluster in three distinct regions: near the cleavage-site cavity, in the four beta-sheets cross-linking the two lobes, and in a solvent-exposed region below the long beta-hairpin in the N-terminal lobe. (ii) We elucidate the role played by the conserved residues for the enzymatic functionality of one representative member of the fold family, the human beta-secretase, by means of classical molecular dynamics (MD). The conserved regions exhibit little overall mobility and yet are involved into the most important modes of structural fluctuations. These modes influence the substrate-catalytic aspartates distance through a relative rotation of the N- and C-terminal lobes. (iii) We investigate the effects of this modulation by estimating the reaction free energy at different representative substrate/enzyme conformations. The activation free energy is strongly affected by large-scale protein motions, similarly to what has been observed in the HIV-1 enzyme. (iv) We extend our findings to all other members of the two eukaryotic and retroviral fold families by recurring to a simple, topology-based, energy functional. This analysis reveals a sophisticated mechanism of enzymatic activity modulation common to all aspartic proteases. We suggest that aspartic proteases have been evolutionarily selected to possess similar functional motions despite the observed fold variations.

Quantum mechanical free energy barrier for an enzymatic reaction

We discuss problems related to in silico studies of enzymes and show that accurate and converged free energy changes for complex chemical reactions can be computed if a method based on a thermodynamic cycle is employed. The method combines the sampling speed of molecular mechanics with the accuracy of a high-level quantum mechanics method. We use the method to compute the free energy barrier for a methyl transfer reaction catalyzed by the enzyme catechol O-methyltransferase at the level of density functional theory. The surrounding protein and solvent are found to have a profound effect on the reaction, and we show that energies can be extrapolated easily from one basis set and exchange-correlation functional to another. Using this procedure we calculate a barrier of 69 kJ/mol, in excellent agreement with the experimental value of 75 kJ/mol.

Potential folding-function interrelationship in proteins

The possibility is addressed that protein folding and function may be related via regions that are critical for both folding and function. This approach is based on the building blocks folding model that describes protein folding as binding events of conformationally fluctuating building blocks. Within these, we identify building block fragments that are critical for achieving the native fold. A library of such critical building blocks (CBBs) is constructed. Then, it is asked whether the functionally important residues fall in these CBB fragments. We find that for over two-thirds of the proteins in our library with available functional information, the catalytic or binding site residues lie within the CBB regions. From the evolutionary standpoint, a folding-function relationship is advantageous, since the need to guard against mutations is limited to one region. Furthermore, conformationally similar CBBs are found in globally unrelated proteins with different functions. Hence, substituting CBBs may lead to designed proteins with altered functions. We further find that the CBBs in our library are conformationally unstable.

HIV-1 protease molecular dynamics of a wild-type and of the V82F/I84V mutant: possible contributions to drug resistance and a potential new target site for drugs

The protease from type 1 human immunodeficiency virus (HIV-1) is a critical drug target against which many therapeutically useful inhibitors have been developed; however, the set of viral strains in the population has been shifting to become more drug-resistant. Because indirect effects are contributing to drug resistance, an examination of the dynamic structures of a wild-type and a mutant could be insightful. Consequently, this study examined structural properties sampled during 22 nsec, all atom molecular dynamics (MD) simulations (in explicit water) of both a wild-type and the drug-resistant V82F/I84V mutant of HIV-1 protease. The V82F/I84V mutation significantly decreases the binding affinity of all HIV-1 protease inhibitors currently used clinically. Simulations have shown that the curling of the tips of the active site flaps immediately results in flap opening. In the 22-nsec MD simulations presented here, more frequent and more rapid curling of the mutant's active site flap tips was observed. The mutant protease's flaps also opened farther than the wild-type's flaps did and displayed more flexibility. This suggests that the effect of the mutations on the equilibrium between the semiopen and closed conformations could be one aspect of the mechanism of drug resistance for this mutant. In addition, correlated fluctuations in the active site and periphery were noted that point to a possible binding site for allosteric inhibitors.

The role of hydrogen bonding in the enzymatic reaction catalyzed by HIV-1 protease

The hydrogen-bond network in various stages of the enzymatic reaction catalyzed by HIV-1 protease was studied through quantum-classical molecular dynamics simulations. The approximate valence bond method was applied to the active site atoms participating directly in the rearrangement of chemical bonds. The rest of the protein with explicit solvent was treated with a classical molecular mechanics model. Two possible mechanisms were studied, general-acid/general-base (GA/GB) with Asp 25 protonated at the inner oxygen, and a direct nucleophilic attack by Asp 25. Strong hydrogen bonds leading to spontaneous proton transfers were observed in both reaction paths. A single-well hydrogen bond was formed between the peptide nitrogen and outer oxygen of Asp 125. The proton was diffusely distributed with an average central position and transferred back and forth on a picosecond scale. In both mechanisms, this interaction helped change the peptide-bond hybridization, increased the partial charge on peptidyl carbon, and in the GA/GB mechanism, helped deprotonate the water molecule. The inner oxygens of the aspartic dyad formed a low-barrier, but asymmetric hydrogen bond; the proton was not positioned midway and made a slightly elongated covalent bond, transferring from one to the other aspartate. In the GA/GB mechanism both aspartates may help deprotonate the water molecule. We observed the breakage of the peptide bond and found that the protonation of the peptidyl amine group was essential for the peptide-bond cleavage. In studies of the direct nucleophilic mechanism, the peptide carbon of the substrate and oxygen of Asp 25 approached as close as 2.3 A.

Uses of enthalpy-entropy compensation in protein research

Cooperative systems of proteins and small molecules form most of biology but are so weakly linked that conventional mass-law formalism requiring exact stoichiometry is inapplicable. The weaknesses cannot be eliminated but using selected families of reactions useful fragmentation of those quantities is often possible. Extra-thermodynamic treatments based on linear-free-energy relationships (LFE) are developed to utilize enthalpy, entropy and volume information not otherwise reliable Linkage systems build around mesophilic proteins are well suited to enforced marriage of linear equations and scaled molecule detail because the ratio of substructure sizes on which folded stability depends is independent of total number of amino-acid residues. Conformational changes in physiological function usually no greater than 0.5 A closely scale to linear thermodynamic changes. The formalisms for use of LFE and compensation relationships are modified to eliminate complications that have previously arisen from incorrect inclusion of the thermal parts of enthalpy and entropy changes in free energy changes. The results are used to remove current confusion about the basis of folded stability in proteins and to minimize the quantitative errors arising from classical treatments of denaturation data. The enthalpy to entropy ratio given by the slope of a compensation plot (its 'compensation temperature') is used to characterize protein construction and function so as to extract machine descriptions of protein linkage systems. In this way the 'fragile' nature of the free-energy surfaces of the myoglobin proteins and the 'strong' character of those surfaces of most other mesophiles can be deduced very simply from the Debye-Waller factors obtained in diffraction studies. The major evolutionary achievement in making proteins big is their crystallike phase behavior. That makes entropy exactly as important as enthalpy so the scalar quantities of small-molecule chemistry can be replaced by the vector quantities that appear necessary to make biology possible.

Accurate and efficient description of protein vibrational dynamics: Comparing molecular dynamics and Gaussian models

Current all-atom potential based molecular dynamics (MD) allows the identification of a protein's functional motions on a wide-range of timescales, up to few tens of nanoseconds. However, functional, large-scale motions of proteins may occur on a timescale currently not accessible by all-atom potential based MD. To avoid the massive computational effort required by this approach, several simplified schemes have been introduced. One of the most satisfactory is the Gaussian network approach based on the energy expansion in terms of the deviation of the protein backbone from its native configuration. Here, we consider an extension of this model that captures in a more realistic way the distribution of native interactions due to the introduction of effective side-chain centroids. Since their location is entirely determined by the protein backbone, the model is amenable to the same exact and computationally efficient treatment as previous simpler models. The ability of the model to describe the correlated motion of protein residues in thermodynamic equilibrium is established through a series of successful comparisons with an extensive (14 ns) MD simulation based on the AMBER potential of HIV-1 protease in complex with a peptide substrate. Thus, the model presented here emerges as a powerful tool to provide preliminary, fast yet accurate characterizations of protein near-native motion.

Conformational fluctuations and electronic properties in myoglobin

In this article we use the recently developed perturbed matrix method (PMM) to investigate the effect of conformational fluctuations on the electronic properties of heme in Myoglobin. This widely studied biomolecule has been chosen as a benchmark for evaluating the accuracy of PMM in a large and complex system. Using a long, 80-ns, molecular dynamics simulation and unperturbed Configuration Interaction (CISD) calculations in PMM, we reproduced the main spectroscopic features of deoxy-Myoglobin. Moreover, in line with our previous results on a photosensitive protein, this study reveals a clear dynamical coupling between electronic properties and conformational fluctuations, suggesting that this correlation could be a general feature of proteins.

Substrate binding to mononuclear metallo-β-lactamase from Bacillus cereus

Structure and dynamics of substrate binding (cefotaxime) to the catalytic pocket of the mononuclear zinc-beta-lactamase from Bacillus cereus are investigated by molecular dynamics simulations. The calculations, which are based on the hydrogen-bond pattern recently proposed by Dal Peraro et al. (J Biol Inorg Chem 2002; 7:704-712), are carried out for both the free and the complexed enzyme. In the resting state, active site pattern and temperature B-factors are in agreement with crystallographic data. In the complexed form, cefotaxime is accommodated into a stable orientation in the catalytic pocket within the nanosecond timescale, interacting with the enzyme zinc-bound hydroxide and the surrounding loops. The beta-lactam ring remains stable and very close to the hydroxide nucleophile agent, giving a stable representation of the productive enzyme-substrate complex.