MDexciteR: Enhanced Sampling Molecular Dynamics by Excited Normal Modes or Principal Components Obtained from Experiments.

Molecular dynamics with excited normal modes (MDeNM) is an enhanced sampling method for exploring conformational changes in proteins with minimal biases. The excitation corresponds to injecting kinetic energy along normal modes describing intrinsic collective motions. Herein, we developed a new automated open-source implementation, MDexciteR (https://github.com/mcosta27/MDexciteR), enabling the integration of MDeNM with two commonly used simulation programs with GPU support. Second, we generalized the method to include the excitation of principal components calculated from experimental ensembles. Finally, we evaluated whether the use of coarse-grained normal modes calculated with elastic network representations preserved the performance and accuracy of the method. The advantages and limitations of these new approaches are discussed based on results obtained for three different protein test cases: two globular and a protein/membrane system.

Concerted conformational dynamics and water movements in the ghrelin G protein-coupled receptor.

There is increasing support for water molecules playing a role in signal propagation through G protein-coupled receptors (GPCRs). However, exploration of the hydration features of GPCRs is still in its infancy. Here, we combined site-specific labeling with unnatural amino acids to molecular dynamics to delineate how local hydration of the ghrelin receptor growth hormone secretagogue receptor (GHSR) is rearranged upon activation. We found that GHSR is characterized by a specific hydration pattern that is selectively remodeled by pharmacologically distinct ligands and by the lipid environment. This process is directly related to the concerted movements of the transmembrane domains of the receptor. These results demonstrate that the conformational dynamics of GHSR are tightly coupled to the movements of internal water molecules, further enhancing our understanding of the molecular bases of GPCR-mediated signaling.

Exploring free energy landscapes of large conformational changes: molecular dynamics with excited normal modes.

Proteins are found in solution as ensembles of conformations in dynamic equilibrium. Exploration of functional motions occurring on micro- to millisecond time scales by molecular dynamics (MD) simulations still remains computationally challenging. Alternatively, normal mode (NM) analysis is a well-suited method to characterize intrinsic slow collective motions, often associated with protein function, but the absence of anharmonic effects preclude a proper characterization of conformational distributions in a multidimensional NM space. Using both methods jointly appears to be an attractive approach that allows an extended sampling of the conformational space. In line with this view, the MDeNM (molecular dynamics with excited normal modes) method presented here consists of multiple-replica short MD simulations in which motions described by a given subset of low-frequency NMs are kinetically excited. This is achieved by adding additional atomic velocities along several randomly determined linear combinations of NM vectors, thus allowing an efficient coupling between slow and fast motions. The relatively high-energy conformations generated with MDeNM are further relaxed with standard MD simulations, enabling free energy landscapes to be determined. Two widely studied proteins were selected as examples: hen egg lysozyme and HIV-1 protease. In both cases, MDeNM provides a larger extent of sampling in a few nanoseconds, outperforming long standard MD simulations. A high degree of correlation with motions inferred from experimental sources (X-ray, EPR, and NMR) and with free energy estimations obtained by metadynamics was observed. Finally, the large sets of conformations obtained with MDeNM can be used to better characterize relevant dynamical populations, allowing for a better interpretation of experimental data such as SAXS curves and NMR spectra.

Conformational Equilibrium of CDK/Cyclin Complexes by Molecular Dynamics with Excited Normal Modes.

Cyclin-dependent kinases (CDKs) and their associated regulatory cyclins are central for timely regulation of cell-cycle progression. They constitute attractive pharmacological targets for development of anticancer therapeutics, since they are frequently deregulated in human cancers and contribute to sustained, uncontrolled tumor proliferation. Characterization of their structural/dynamic features is essential to gain in-depth insight into structure-activity relationships. In addition, the identification of druggable pockets or key intermediate conformations yields potential targets for the development of novel classes of inhibitors. Structural studies of CDK2/cyclin A have provided a wealth of information concerning monomeric/heterodimeric forms of this kinase. There is, however, much less structural information for other CDK/cyclin complexes, including CDK4/cyclin D1, which displays an alternative (open) position of the cyclin partner relative to CDK, contrasting with the closed CDK2/cyclin A conformation. In this study, we carried out normal-mode analysis and enhanced sampling simulations with our recently developed method, molecular dynamics with excited normal modes, to understand the conformational equilibrium on these complexes. Interestingly, the lowest-frequency normal mode computed for each complex described the transition between the open and closed conformations. Exploration of these motions with an explicit-solvent representation using molecular dynamics with excited normal modes confirmed that the closed conformation is the most stable for the CDK2/cyclin A complex, in agreement with their experimentally available structures. On the other hand, we clearly show that an open↔closed equilibrium may exist in CDK4/cyclin D1, with closed conformations resembling that captured for CDK2/cyclin A. Such conformational preferences may result from the distinct distributions of frustrated contacts in each complex. Using the same approach, the putative roles of the Thr(160) phosphoryl group and the T-loop conformation were investigated. These results provide a dynamic view of CDKs revealing intermediate conformations not yet characterized for CDK members other than CDK2, which will be useful for the design of inhibitors targeting critical conformational transitions.

Impact of M36I polymorphism on the interaction of HIV-1 protease with its substrates: insights from molecular dynamics.

BACKGROUND: Over the last decades, a vast structural knowledge has been gathered on the HIV-1 protease (PR). Noticeably, most of the studies focused the B-subtype, which has the highest prevalence in developed countries. Accordingly, currently available anti-HIV drugs target this subtype, with considerable benefits for the corresponding patients. RESULTS: Herein, we used molecular dynamics simulations to investigate the role of this polymorphism on the interaction of PR with six of its natural cleavage-sites substrates. CONCLUSIONS: With multiple approaches and analyses we identified structural and dynamical determinants associated with the changes found in the binding affinity of the M36I variant. This mutation influences the flexibility of both PR and its complexed substrate. The observed impact of M36I, suggest that combination with other non-B subtype polymorphisms, could lead to major effects on the interaction with the 12 known cleavage sites, which should impact the virion maturation.

Free Energy Profiles along Consensus Normal Modes Provide Insight into HIV-1 Protease Flap Opening.

Describing biological macromolecular energetics from computer simulations can pose major challenges, and often necessitates enhanced conformational sampling. We describe the calculation of conformational free-energy profiles along carefully chosen collective coordinates: "consensus" normal modes, developed recently as robust alternatives to conventional normal modes. In an application to the HIV-1 protease, we obtain efficient sampling of significant flap opening movements governing inhibitor binding from relatively short simulations, in close correspondence with experimental results.

New parameterization approaches of the LIE method to improve free energy calculations of PlmII-Inhibitors complexes.

The standard parameterization of the Linear Interaction Energy (LIE) method has been applied with quite good results to reproduce the experimental absolute binding free energies for several protein-ligand systems. However, we found that this parameterization failed to reproduce the experimental binding free energy of Plasmepsin II (PlmII) in complexes with inhibitors belonging to four dissimilar scaffolds. To overcome this fact, we developed three approaches of LIE, which combine systematic approaches to predict the inhibitor-specific values of α, ß, and γ parameters, to gauge their ability to calculate the absolute binding free energies for these PlmII-Inhibitor complexes. Specifically: (i) we modified the linear relationship between the weighted nonpolar desolvation ratio (WNDR) and the α parameter, by introducing two models of the ß parameter determined by the free energy perturbation (FEP) method in the absence of the constant term γ, and (ii) we developed a new parameterization model to investigate the linear correlation between WNDR and the correction term γ. Using these parameterizations, we were able to reproduce the experimental binding free energy from these systems with mean absolute errors lower than 1.5 kcal/mol.

Understanding the HIV-1 protease nelfinavir resistance mutation D30N in subtypes B and C through molecular dynamics simulations.

A major concern in the antiretroviral (ARV) treatment of HIV infections with protease inhibitors (PI) is the emergence of resistance, which results from the selection of distinct mutations within the viral protease (PR) gene. Among patients who do not respond to treatment with the PI nelfinavir (NFV), the D30N mutation is often observed. However, several reports have shown that D30N emerges with different frequencies in distinct HIV-1 genetic forms or subtypes. In the present work, we analyzed the binding of NFV and the Gag substrate CA/p2 to PR from HIV-1 subtypes B and C through molecular dynamics (MD) simulations. The wild-type and drug-resistant D30N mutants were investigated in both subtypes. The compensatory mutations N83T and N88D, observed in vitro and in vivo when subtype C acquires D30N, were also studied. D30N appears to facilitate conformational changes in subtype B PR, but not in that from subtype C, and this could be associated with disestablishment of an alpha-helical region of the PR. Furthermore, the total contact areas of NFV or the CA/p2 substrate with the mutant PR correlated with changes in the resistance patterns and replicative capacity. Finally, we observed in our MD simulations that mutant PR proteins show different patterns for hydrophobic/van der Waals contact. These findings suggest that different molecular mechanisms contribute to resistance, and we propose that a single mutation has distinct impacts on different HIV-1 subtypes.

How does heparin prevent the pH inactivation of cathepsin B? Allosteric mechanism elucidated by docking and molecular dynamics.

BACKGROUND: Cathepsin B (catB) is a promising target for anti-cancer drug design due to its implication in several steps of tumorigenesis. catB activity and inhibition are pH-dependent, making it difficult to identify efficient inhibitor candidates for clinical trials. In addition it is known that heparin binding stabilizes the enzyme in alkaline conditions. However, the molecular mechanism of stabilization is not well understood, indicating the need for more detailed structural and dynamic studies in order to clarify the influence of pH and heparin binding on catB stability. RESULTS: Our pKa calculations of catB titratable residues revealed distinct protonation states under different pH conditions for six key residues, of which four lie in the crucial interdomain interface. This implies changes in the overall charge distribution at the catB surface, as revealed by calculation of the electrostatic potential. We identified two basic surface regions as possible heparin binding sites, which were confirmed by docking calculations. Molecular dynamics (MD) of both apo catB and catB-heparin complexes were performed using protonation states for catB residues corresponding to the relevant acidic or alkaline conditions. The MD of apo catB at pH 5.5 was very stable, and presented the highest number and occupancy of hydrogen bonds within the inter-domain interface. In contrast, under alkaline conditions the enzyme's overall flexibility was increased: interactions between active site residues were lost, helical content decreased, and domain separation was observed as well as high-amplitude motions of the occluding loop - a main target of drug design studies. Essential dynamics analysis revealed that heparin binding modulates large amplitude motions promoting rearrangement of contacts between catB domains, thus favoring the maintenance of helical content as well as active site stability. CONCLUSIONS: The results of our study contribute to unraveling the molecular events involved in catB inactivation in alkaline pH, highlighting the fact that protonation changes of few residues can alter the overall dynamics of an enzyme. Moreover, we propose an allosteric role for heparin in the regulation of catB stability in such a manner that the restriction of enzyme flexibility would allow the establishment of stronger contacts and thus the maintenance of overall structure.

Predicting functional residues in Plasmodium falciparum plasmepsins by combining sequence and structural analysis with molecular dynamics simulations.

Plasmepsins are aspartic proteases involved in the initial steps of the hemoglobin degradation pathway, a critical stage in the Plasmodium falciparum life cycle during human infection. Thus, they are attractive targets for novel therapeutic compounds to treat malaria, which remains one of the world's biggest health problems. The three-dimensional structures available for P. falciparum plasmepsins II and IV make structure-based drug design of antimalarial compounds that focus on inhibiting plasmepsins possible. However, the structural flexibility of the plasmepsin active site cavity combined with insufficient knowledge of the functional residues and of those determining the specificity of parasitic enzymes is a drawback when designing specific inhibitors. In this study, we have combined a sequence and structural analysis with molecular dynamics simulations to predict the functional residues in P. falciparum plasmepsins. The careful analysis of X-ray structures and 3D models carried out here suggests that residues Y17, V105, T108, L191, L242, Q275, and T298 are important for plasmepsin function. These seven amino acids are conserved across the malarial strains but not in human aspartic proteases. Residues V105 and T108 are localized in a flap of an interior pocket and they only establish contacts with a specific non-peptide achiral inhibitor. We also observed a rapid conformational change in the L3 region of plasmepsins that closes the active site of the enzyme, which explains earlier experimental findings. These results shed light on the role of V105 and T108 residues in plasmepsin specificities, and they should be useful in structure-based design of novel, selective inhibitors that may serve as antimalarial drugs.

Molecular dynamics simulations applied to the study of subtypes of HIV-1 protease common to Brazil, Africa, and Asia.

Africa accounts for the majority of HIV-1 infections worldwide caused mainly by the A and C viral subtypes rather than B subtype, which prevails in the United States and Western Europe. In Brazil, B subtype is the major subtype, but F, C, and A also circulate. These non-B subtypes present polymorphisms, and some of them occur at sites that have been associated with drug resistance, including the HIV-1 protease (PR), one important drug target. Here, we report a Molecular Dynamics study of the B and non-B PR complexed with the inhibitor ritonavir to delineate the behavior of each subtype. We compare root mean squared deviation, binding free energy by linear interaction energy approach, hydrogen bonds, and intermolecular contact surface area between inhibitor and PR. From our results, we can provide a basis to understand the molecular mechanism of drug resistance in non-B subtypes. In this sense, we found a decrease of approx 4 kcal/mol in deltaG of binding between B and non-B subtypes. This corresponds to the loss of one hydrogen bond, which is in agreement with our H-bond analysis. Previous experimental affinity studies reported analogous results with inhibition constant values for non-B PR.

Interplay between parasite cysteine proteases and the host kinin system modulates microvascular leakage and macrophage infection by promastigotes of the Leishmania donovani complex.

Kinins, the vasoactive peptides proteolytically liberated from kininogens, were recently recognized as signals alerting the innate immune system. Here we demonstrate that Leishmania donovani and Leishmania chagasi, two etiological agents of visceral leishmaniasis (VL), activate the kinin system. Intravital microscopy in the hamster cheek pouch showed that topically applied promastigotes induced macromolecular leakage (FITC-dextran) through postcapillary venules. Peaking at 15 min, the parasite-induced leakage was drastically enhanced by captopril (Cap), an inhibitor of angiotensin-converting enzyme (ACE), a kinin-degrading metallopeptidase. The enhanced microvascular responses were cancelled by HOE-140, an antagonist of the B2 bradykinin receptor (B2R), or by pre-treatment of promastigotes with the irreversible cysteine proteinase inhibitor N-methylpiperazine-urea-Phe-homoPhe-vinylsulfone-benzene (N-Pip-hF-VSPh). In agreement with the above-mentioned data, the promastigotes vigorously induced edema in the paw of Cap-treated J129 mice, but not Cap-B2R-/- mice. Analysis of parasite-induced breakdown of high molecular weight kininogens (HK), combined with active site-affinity-labeling with biotin-N-Pip-hF-VSPh, identified 35-40 kDa proteins as kinin-releasing cysteine peptidases. We then checked if macrophage infectivity was influenced by interplay between these kinin-releasing parasite proteases, kininogens, and kinin-degrading peptidases (i.e. ACE). Our studies revealed that full-fledged B2R engagement resulted in vigorous increase of L. chagasi uptake by resident macrophages. Evidence that inflammatory macrophages treated with HOE-140 became highly susceptible to amastigote outgrowth, assessed 72 h after initial macrophage interaction, further suggests that the kinin/B2R activation pathway may critically modulate inflammation and innate immunity in visceral leishmaniasis.