Equilibrium landscape of ingress/egress channels and gating residues of the Cytochrome P450 3A4.

The Cytochrome P450 (CYP) enzymes metabolize a variety of drugs, which may potentially lead to toxicity or reduced efficacy when drugs are co-administered. These drug-drug interactions are often manifested by CYP3A4, the most prevalent of all CYP isozymes. We carried out multiple MD simulations employing CAVER to quantify the channels, and Hidden Markov Models (HMM) to characterize the behavior of the gating residues. We discuss channel properties, bottleneck residues with respect to their likelihood to deem the respective channel ingress or egress, gating residues regarding their open or closed states, and channel location relative to the membrane. Channels do not display coordinated motion and randomly transition between different conformations. Gateway residues also behave in a random fashion. Our findings shed light on the equilibrium behavior of the gating residues and channels in the apo state.

Molecular dynamics simulations of an engineered T4 lysozyme exclude helix to sheet transition, and provide insights into long distance, intra-protein switchable motion.

An engineered variant of T4 lysozyme serves as a model for studying induced remote conformational changes in a full protein context. The design involves a duplicated surface helix, flanked by two loops, that switches between two different conformations spanning about 20 Å. Molecular dynamics simulations of the engineered protein, up to 1 µs, rule out α-helix to ß-sheet transitions within the duplicated helix as suggested by others. These simulations highlight how the use of different force fields can lead to radical differences in the structure of the protein. In addition, Markov state modeling and transition path theory were employed to map a 6.6 µs simulation for possible early intermediate states and to provide insights into the onset of the switching motion. The putative intermediates involve the folding of one helical turn in the C-terminal loop through energy driven, sequential rearrangement of nearby salt bridges around the key residue Arg63. These results provide a first step towards understanding the energetics and dynamics of a rather complicated intra-protein motion.

The NS4A Cofactor Dependent Enhancement of HCV NS3 Protease Activity Correlates with a 4D Geometrical Measure of the Catalytic Triad Region.

We are developing a 4D computational methodology, based on 3D structure modeling and molecular dynamics simulation, to analyze the active site of HCV NS3 proteases, in relation to their catalytic activity. In our previous work, the 4D analyses of the interactions between the catalytic triad residues (His57, Asp81, and Ser139) yielded divergent, gradual and genotype-dependent, 4D conformational instability measures, which strongly correlate with the known disparate catalytic activities among genotypes. Here, the correlation of our 4D geometrical measure is extended to intra-genotypic alterations in NS3 protease activity, due to sequence variations in the NS4A activating cofactor. The correlation between the 4D measure and the enzymatic activity is qualitatively evident, which further validates our methodology, leading to the development of an accurate quantitative metric to predict protease activity in silico. The results suggest plausible "communication" pathways for conformational propagation from the activation subunit (the NS4A cofactor binding site) to the catalytic subunit (the catalytic triad). The results also strongly suggest that the well-sampled (via convergence quantification) structural dynamics are more connected to the divergent catalytic activity observed in HCV NS3 proteases than to rigid structures. The method could also be applicable to predict patients' responses to interferon therapy and better understand the innate interferon activation pathway.

Comparative molecular dynamics simulation of Hepatitis C Virus NS3/4A protease (Genotypes 1b, 3a and 4b) predicts conformational instability of the catalytic triad in drug resistant strains.

The protease domain of the Hepatitis C Virus (HCV) nonstructural protein 3 (NS3) has been targeted for inhibition by several direct-acting antiviral drugs. This approach has had marked success to treat infections caused by HCV genotype 1 predominant in the USA, Europe, and Japan. However, genotypes 3 and 4, dominant in developing countries, are resistant to a number of these drugs and little progress has been made towards understanding the structural basis of their drug resistivity. We have previously developed a 4D computational methodology, based on 3D structure modeling and molecular dynamics simulation, to analyze the active sites of the NS3 proteases of HCV-1b and 4a in relation to their catalytic activity and drug susceptibility. Here, we improved the methodology, extended the analysis to include genotype 3a (predominant in South Asia including Pakistan), and compared the results of the three genotypes (1b, 3a and 4a). The 4D analyses of the interactions between the catalytic triad residues (His57, Asp81, and Ser139) indicate conformational instability of the catalytic site in HCV-3a and 4a compared to that of HCV-1b NS3 protease. The divergence is gradual and genotype-dependent, with HCV-1b being the most stable, HCV-4a being the most unstable and HCV-3a representing an intermediate state. These results suggest that the structural dynamics behavior, more than the rigid structure, could be related to the altered catalytic activity and drug susceptibility seen in NS3 proteases of HCV-3a and 4a.

A 3D structural model and dynamics of hepatitis C virus NS3/4A protease (genotype 4a, strain ED43) suggest conformational instability of the catalytic triad: implications in catalysis and drug resistivity.

Egypt has the highest prevalence of hepatitis C virus (HCV) infection worldwide with a frequency of 15%. More than 90% of these infections are due to genotype 4, and the subtype 4a (HCV-4a) predominates. Moreover, due to the increased mobility of people, HCV-4a has recently spread to several European countries. The protease domain of the HCV nonstructural protein 3 (NS3) has been targeted for inhibition by several drugs. This approach has had marked success in inhibiting genotype 1 (HCV-1), the predominant genotype in the USA, Europe, and Japan. However, HCV-4a was found to resist inhibition by a number of these drugs, and little progress has been made to understand the structural basis of its drug resistivity. As a step forward, we sequenced the NS3 HCV-4a protease gene (strain ED43) and subsequently built a 3D structural model threaded through a template crystal structure of HCV-1b NS3 protease. The model protease, HCV-4a, shares 83% sequence identity with the template protease, HCV-1b, and has nearly identical rigid structural features. Molecular dynamics simulations predict similar overall dynamics of the two proteases. However, local dynamics and 4D analysis of the interactions between the catalytic triad residues (His57, Asp81, and Ser139) indicate conformational instability of the catalytic site in HCV-4a NS3 protease. These results suggest that the divergent dynamics behavior, more than the rigid structure, could be related to the altered catalytic activity and drug resistivity seen in HCV-4a.