Identification of potential anti-hyperglycemic compounds in Cordyceps militaris ethyl acetate extract: in vitro and in silico studies.

Cordyceps militaris has been long known for valuable health benefits by folk experience and was recently reported with diabetes-tackling evidences, thus deserving extending efforts on screening for component-activity relationship. In this study, experiments were carried out to find the evidence, justification, and input for computations on the potential against diabetes-related protein structures: PDB-4W93, PDB-3W37, and PDB-4A3A. Liquid chromatography identified 14 bioactive compounds in the ethyl acetate extract (1-14) and quantified the contents of cordycepin (0.11%) and adenosine (0.01%). Bioassays revealed the overall potential of the extract against α-amylase (IC50 = 6.443 ± 0.364 mg.mL-1) and α-glucosidase (IC50 = 2.580 ± 0.194 mg.mL-1). A combination of different computational platforms was used to select the most promising candidates for applications as anti-diabetic bio-inhibitors, i.e. 1 (ground state: -888.49715 a.u.; dipole moment 3.779 Debye; DS¯ -12.3 kcal.mol-1; polarizability 34.7 Å3; logP - 1.30), 10 (ground state: -688.52406 a.u.; dipole moment 5.487 Debye; DS¯ -12.6 kcal.mol-1; polarizability 24.9 Å3; logP - 3.39), and 12 (ground state: -1460.07276 a.u.; dipole moment 3.976 Debye; DS¯ -12.5 kcal.mol-1; polarizability 52.4 Å3; logP - 4.39). The results encourage further experimental tests on cordycepin (1), mannitol (10), and adenosylribose (12) to validate their in-practice diabetes-related activities, thus conducive to hypoglycemic applications.Communicated by Ramaswamy H. Sarma.

Molecular dynamics simulation of cancer cell membrane perforated by shockwave induced bubble collapse.

It has been widely accepted that cancer cells are softer than their normal counterparts. This motivates us to propose, as a proof-of-concept, a method for the efficient delivery of therapeutic agents into cancer cells, while normal cells are less affected. The basic idea of this method is to use a water jet generated by the collapse of the bubble under shockwaves to perforate pores in the cell membrane. Given a combination of shockwave and bubble parameters, the cancer membrane is more susceptible to bending, stretching, and perforating than the normal membrane because the bending modulus of the cancer cell membrane is smaller than that of the normal cell membrane. Therefore, the therapeutic agent delivery into cancer cells is easier than in normal cells. Adopting two well-studied models of the normal and cancer membranes, we perform shockwave induced bubble collapse molecular dynamics simulations to investigate the difference in the response of two membranes over a range of shockwave impulse 15-30 mPa s and bubble diameter 4-10 nm. The simulation shows that the presence of bubbles is essential for generating a water jet, which is required for perforation; otherwise, pores are not formed. Given a set of shockwave impulse and bubble parameters, the pore area in the cancer membrane is always larger than that in the normal membrane. However, a too strong shockwave and/or too large bubble results in too fast disruption of membranes, and pore areas are similar between two membrane types. The pore closure time in the cancer membrane is slower than that in the normal membrane. The implications of our results for applications in real cells are discussed in some details. Our simulation may be useful for encouraging future experimental work on novel approaches for cancer treatment.

Elastic moduli of normal and cancer cell membranes revealed by molecular dynamics simulations.

Recent studies indicate that there are mechanical differences between normal cells and cancer cells. Because the cell membrane takes part in a variety of vital processes, we test the hypothesis of whether or not two fundamental alterations in the cell membrane, i.e., the overexpression of phosphatidylserine lipids in the outer leaflet and a reduction in cholesterol concentration, could cause the softening in cancer cells. Adopting ten models of normal and cancer cell membranes, we carry out 1 µs all-atom molecular dynamics simulations to compare the structural properties and elasticity properties of two membrane types. We find that the overexpression of the phosphatidylserine lipids in the outer leaflet does not significantly alter the area per lipid, the membrane thickness, the lipid order parameters and the elasticity moduli of the cancer membranes. However, a reduction in the cholesterol concentration leads to clear changes in those quantities, especially decreases in the bending, tilt and twist moduli. This implies that the reduction of cholesterol concentration in the cancer membranes could contribute to the softening of cancer cells.

Effect of Surface Roughness on Aggregation of Polypeptide Chains: A Monte Carlo Study.

The self-assembly of amyloidogenic peptides and proteins into fibrillar structures has been intensively studied for several decades, because it seems to be associated with a number of neurodegenerative diseases, such as Alzheimer's and Parkinson's disease. Therefore, understanding the molecular mechanisms of this phenomenon is important for identifying an effective therapy for the corresponding diseases. Protein aggregation in living organisms very often takes place on surfaces like membranes and the impact of a surface on this process depends not only on the surface chemistry but also on its topology. Our goal was to develop a simple lattice model for studying the role of surface roughness in the aggregation kinetics of polypeptide chains and the morphology of aggregates. We showed that, consistent with the experiment, an increase in roughness slows down the fibril formation, and this process becomes inhibited at a very highly level of roughness. We predicted a subtle catalytic effect that a slightly rough surface promotes the self-assembly of polypeptide chains but does not delay it. This effect occurs when the interaction between the surface and polypeptide chains is moderate and can be explained by taking into account the competition between energy and entropy factors.

Amyloid Oligomers: A Joint Experimental/Computational Perspective on Alzheimer's Disease, Parkinson's Disease, Type II Diabetes, and Amyotrophic Lateral Sclerosis.

Protein misfolding and aggregation is observed in many amyloidogenic diseases affecting either the central nervous system or a variety of peripheral tissues. Structural and dynamic characterization of all species along the pathways from monomers to fibrils is challenging by experimental and computational means because they involve intrinsically disordered proteins in most diseases. Yet understanding how amyloid species become toxic is the challenge in developing a treatment for these diseases. Here we review what computer, in vitro, in vivo, and pharmacological experiments tell us about the accumulation and deposition of the oligomers of the (Aß, tau), α-synuclein, IAPP, and superoxide dismutase 1 proteins, which have been the mainstream concept underlying Alzheimer's disease (AD), Parkinson's disease (PD), type II diabetes (T2D), and amyotrophic lateral sclerosis (ALS) research, respectively, for many years.

Remdesivir Strongly Binds to Both RNA-Dependent RNA Polymerase and Main Protease of SARS-CoV-2: Evidence from Molecular Simulations.

The outbreak of a new coronavirus SARS-CoV-2 (severe acute respiratory syndrome-coronavirus 2) has caused a global COVID-19 (coronavirus disease 2019) pandemic, resulting in millions of infections and thousands of deaths around the world. There is currently no drug or vaccine for COVID-19, but it has been revealed that some commercially available drugs are promising, at least for treating symptoms. Among them, remdesivir, which can block the activity of RNA-dependent RNA polymerase (RdRp) in old SARS-CoV and MERS-CoV viruses, has been prescribed to COVID-19 patients in many countries. A recent experiment showed that remdesivir binds to SARS-CoV-2 with an inhibition constant of µM, but the exact target has not been reported. In this work, combining molecular docking, steered molecular dynamics, and umbrella sampling, we examined its binding affinity to two targets including the main protease (Mpro), also known as 3C-like protease, and RdRp. We showed that remdesivir binds to Mpro slightly weaker than to RdRp, and the corresponding inhibition constants, consistent with the experiment, fall to the µM range. The binding mechanisms of remdesivir to two targets differ in that the electrostatic interaction is the main force in stabilizing the RdRp-remdesivir complex, while the van der Waals interaction dominates in the Mpro-remdesivir case. Our result indicates that remdesivir can target not only RdRp but also Mpro, which can be invoked to explain why this drug is effective in treating COVID-19. We have identified residues of the target protein that make the most important contribution to binding affinity, and this information is useful for drug development for this disease.

How Good is Jarzynski's Equality for Computer-Aided Drug Design?

Accurate determination of the binding affinity of the ligand to the receptor remains a difficult problem in computer-aided drug design. Here, we study and compare the efficiency of Jarzynski's equality (JE) combined with steered molecular dynamics and the linear interaction energy (LIE) method by assessing the binding affinity of 23 small compounds to six receptors, including ß-lactamase, thrombin, factor Xa, HIV-1 protease (HIV), myeloid cell leukemia-1, and cyclin-dependent kinase 2 proteins. It was shown that Jarzynski's nonequilibrium binding free energy ΔGneqJar correlates with the available experimental data with the correlation levels R = 0.89, 0.86, 0.83, 0.80, 0.83, and 0.81 for six data sets, while for the binding free energy ΔGLIE obtained by the LIE method, we have R = 0.73, 0.80, 0.42, 0.23, 0.85, and 0.01. Therefore, JE is recommended to be used for ranking binding affinities as it provides accurate and robust results. In contrast, LIE is not as reliable as JE, and it should be used with caution, especially when it comes to new systems.

7-Methoxytacrine and 2-Aminobenzothiazole Heterodimers: Structure-Mechanism Relationship of Amyloid Inhibitors Based on Rational Design.

The formation and accumulation of amyloid aggregates are the phenomena that accompany amyloidoses, which are currently untreatable and include Alzheimer's and Parkinson's diseases, diabetes mellitus, non-neuropathic lysozyme systemic amyloidosis, and others. One of the very promising therapeutic approaches seems to be an inhibition of amyloid formation and/or clearance of amyloid aggregates. Small molecules have a great potential to interfere with amyloid fibrillation of peptides and polypeptides, which can be improved by connection of cyclic structures into single multicyclic molecules and their dimerization. In our study, we focused on heterodimers consisting of 7-methoxytacrine (7-MEOTA) and 2-aminobenzothiazole (BTZ) parent molecules connected by an aliphatic linker. Using in vitro and in silico methods, we investigated the ability of studied compounds to inhibit the amyloid aggregation of hen egg white lysozyme. Heterodimerization led to significant improvement of inhibitory activity compared to that of the parent molecules. The efficiency of the heterodimers varied; the most effective inhibitor contained the longest linker, eight carbons long. We suggest that binding of a heterodimer to a lysozyme blocks the interaction between the ß-domain and C-helix region essential for the formation of amyloid cross-ß structure. Elongation of the linker ultimately enhances the compound's ability to prevent this interaction by allowing the BTZ part of the heterodimer to bind more effectively, increasing the compound's binding affinity, and also by greater steric obstruction. This study represents an important contribution to the recent rational design of potential lead small molecules with anti-amyloid properties, and the heterodimers studied are prospective candidates for the treatment of systemic lysozyme amyloidosis and other amyloid-related diseases.

G37V mutation of Aß42 induces a nontoxic ellipse-like aggregate: An in vitro and in silico study.

The glycine zipper motif at the C-terminus of the ß-amyloid (Aß) peptide have been shown to strongly influence the formation of neurotoxic aggregates. A previous study showed that the G37L mutation dramatically reduces the Aß toxicity in vivo and in vitro. However, the primary cause and mechanism of the glycine zipper motif on Aß properties remain unknown. To gain molecular insights into the impact of glycine zipper on Aß properties, we substituted the residue 37 of Glycine by Valine and studied the structural and biochemical properties of G37V mutation, Aß42(37V), by using in vitro and in silico approaches. Unlike G37L mutation, the G37V mutation reduced toxicity substantially but did not significantly accelerate the aggregation rate or change the content of secondary structures. Further TEM analyses showed that the G37V mutation formed an ellipse-like aggregate rather than a network-like fibril as wild type or G37L mutation of Aß42 form. This different aggregation morphology may be highly linked with the reduction of toxicity. To gain the insight for the different properties of Aß42(37V), we studied the structure of Aß42 and G37V mutation using the replica exchange molecular dynamics simulation. Our results demonstrate that although the overall secondary structure population is similar with Aß42 and Aß42(G37V), Aß42(G37V) shows an increase in the ß-turn and ß-hairpin at residues 36-37 and the flexibility of the Asp23-Lys28 salt bridge. These unique structural features may be the possible reason to account for the ellipse-like morphology.

Bexarotene cannot reduce amyloid beta plaques through inhibition of production of amyloid beta peptides: in silico and in vitro study.

Recently, it has been reported that anti-cancer drug bexarotene can remarkably destroy amyloid beta (Aß) plaques in mouse models suggesting therapeutic potential for Alzheimer's disease. However, the effect of bexarotene on clearance of plaques has not been seen in some mouse models. One of the possible mechanisms explaining this phenomenon is that bexarotene levels up expression of apolipoprotein 4 (ApoE4) leading to intracellular clearance of Aß peptide. Therefore, an interesting question emerges of whether bexarotene can destroy Aß plaques by direct interaction with them or by preventing production of Aß peptides. In our previous work we have shown that bexarotene cannot clear amyloid aggregates due to their weak interaction using in silico and in vitro experiments. Here we explore the possibility of inhibiting Aß production through bexarotene binding to ß-secretase which can cleave Aß peptides from amyloid precursor protein. Using the molecular mechanics-Poisson-Boltzmann surface area method and all-atom simulations we have shown that bexarotene has a very low binding affinity to ß-secretase. This result has been also confirmed by our in vitro experiment implying that bexarotene cannot clear amyloid plaques through inhibition of Aß production. We have also shown that bexarotene tightly binds to both peroxisome proliferator-activated receptor γ (PPAR-γ) and retinoid X receptors (RXRs). Thus, our result does not contradict the hypothesis that the reduction of Aß plaques occurs due to bexarotene-induced overexpression of ApoE4.

Kinetics and mechanical stability of the fibril state control fibril formation time of polypeptide chains: A computational study.

Fibril formation resulting from protein misfolding and aggregation is a hallmark of several neurodegenerative diseases such as Alzheimer's and Parkinson's diseases. Despite much progress in the understanding of the protein aggregation process, the factors governing fibril formation rates and fibril stability have not been fully understood. Using lattice models, we have shown that the fibril formation time is controlled by the kinetic stability of the fibril state but not by its energy. Having performed all-atom explicit solvent molecular dynamics simulations with the GROMOS43a1 force field for full-length amyloid beta peptides Aß40 and Aß42 and truncated peptides, we demonstrated that kinetic stability can be accessed via mechanical stability in such a way that the higher the mechanical stability or the kinetic stability, the faster the fibril formation. This result opens up a new way for predicting fibril formation rates based on mechanical stability that may be easily estimated by steered molecular dynamics.

Generalization of the elastic network model for the study of large conformational changes in biomolecules.

The elastic network (EN) is a prime model that describes the long-time dynamics of biomolecules. However, the use of harmonic potentials renders this model insufficient for studying large conformational changes of proteins (e.g. stretching of proteins, folding and thermal unfolding). Here, we extend the capabilities of the EN model by using a harmonic approximation described by Lennard-Jones (LJ) interactions for far contacts and native contacts obtained from the standard overlap criterion as in the case of Go-like models. While our model is validated against the EN model by reproducing the equilibrium properties for a number of proteins, we also show that the model is suitable for the study of large conformation changes by providing various examples. In particular, this is illustrated on the basis of pulling simulations that predict with high accuracy the experimental data on the rupture force of the studied proteins. Furthermore, in the case of DDFLN4 protein, our pulling simulations highlight the advantages of our model with respect to Go-like approaches, where the latter fail to reproduce previous results obtained by all-atom simulations that predict an additional characteristic peak for this protein. In addition, folding simulations of small peptides yield different folding times for α-helix and ß-hairpin, in agreement with experiment, in this way providing further opportunities for the application of our model in studying large conformational changes of proteins. In contrast to the EN model, our model is suitable for both normal mode analysis and molecular dynamics simulation. We anticipate that the proposed model will find applications in a broad range of problems in biology, including, among others, protein folding and thermal unfolding.

Erythromycin leads to differential protein expression through differences in electrostatic and dispersion interactions with nascent proteins.

The antibiotic activity of erythromycin, which reversibly binds to a site within the bacterial ribosome exit tunnel, against many gram positive microorganisms indicates that it effectively inhibits the production of proteins. Similar to other macrolides, the activity of erythromycin is far from universal, as some peptides can bypass the macrolide-obstructed exit tunnel and become partially or fully synthesized. It is unclear why, at the molecular level, some proteins can be synthesized while others cannot. Here, we use steered molecular dynamics simulations to examine how erythromycin inhibits synthesis of the peptide ErmCL but not the peptide H-NS. By pulling these peptides through the exit tunnel of the E.coli ribosome with and without erythromycin present, we find that erythromycin directly interacts with both nascent peptides, but the force required for ErmCL to bypass erythromycin is greater than that of H-NS. The largest forces arise three to six residues from their N-terminus as they start to bypass Erythromycin. Decomposing the interaction energies between erythromycin and the peptides at this point, we find that there are stronger electrostatic and dispersion interactions with the more C-terminal residues of ErmCL than with H-NS. These results suggest that erythromycin slows or stalls synthesis of ErmCL compared to H-NS due to stronger interactions with particular residue positions along the nascent protein.

Bexarotene Does Not Clear Amyloid Beta Plaques but Delays Fibril Growth: Molecular Mechanisms.

In 2012, it was reported that anticancer drug bexarotene reduced amyloid plaque and improved mental functioning in a small sample of mice engineered to exhibit Alzheimer's like symptoms. It has been suggested that bexarotene stimulates expression of apolipoprotein E (ApoE) leading to intracellular clearance of amyloid beta (Aß). However, the effect of bexarotene on clearance of plaques has not been seen in some mouse models. Two interesting questions include whether bexarotene can destroy Aß fibrils via direct interaction with them and how this compound impacts the lag phase in the fibril growth process. By the Thioflavin T fluorescence assay and atomic force microscopy, we have shown that bexarotene prolongs the lag phase, but it does not degrade Aß fibrils. The impotence of bexarotene in destroying fibrils means that this compound is weakly bound to Aß. On the other hand, the weak binding would prevent bexarotene from prolonging the lag phase. Thus, our two main in vitro observations seem to contradict each other. In order to settle this problem at the atomic level, we have performed all-atom molecular dynamics simulations in explicit water. We have demonstrated that bexarotene is not capable to reduce amyloid deposits due to weak binding to Aß fibrils. However, it delays the self-assembly through reduction of the ß-content of Aß monomers at high enough ligand concentrations. Bexarotene is the first compound which displays such an unusual behavior. We have also shown that bexarotene has a low binding propensity to Aß monomer and dimer.

Fullerenol C60(OH)16 prevents amyloid fibrillization of Aß40-in vitro and in silico approach.

The generation of Aß amyloid aggregates in the form of senile plaques in the brain is one of the pathological hallmarks of Alzheimer's disease (AD). There is no cure for AD and one of the recent treatment strategies is focused on the inhibition of amyloid fibrillization of Aß peptide. Fullerene C60 has been proposed as a candidate for destroying Aß aggregates but it is not soluble in water and its toxicity to cells remains largely ambiguous. To overcome these drawbacks, we synthesized and studied the effect of water-soluble fullerenol C60(OH)16 (fullerene C60 carrying 16 hydroxyl groups) on the amyloid fibrillization of Aß40 peptide in vitro. Using a Thioflavin T fluorescent assay and atomic force microscopy it was found that C60(OH)16 effectively reduces the formation of amyloid fibrils. The IC50 value is in the low range (µg ml(-1)) suggesting that fullerenol interferes with Aß40 aggregation at stoichiometric concentrations. The in silico calculations supported the experimental data. It was revealed that fullerenol tightly binds to monomer Aß40 and polar, negatively charged amino acids play a key role. Electrostatic interactions dominantly contribute to the binding propensity via interaction of the oxygen atoms from the COO(-) groups of side chains of polar, negatively charged amino acids with the OH groups of fullerenol. This stabilizes contact with either the D23 or K28 of the salt bridge. Due to the lack of a well-defined binding pocket fullerenol is also inclined to locate near the central hydrophobic region of Aß40 and can bind to the hydrophobic C-terminal of the peptide. Upon fullerenol binding the salt bridge becomes flexible, inhibiting Aß aggregation. In order to assess the toxicity of fullerenol, we found that exposure of neuroblastoma SH-SY5Y cells to fullerenol caused no significant changes in viability after 24 h of treatment. These results suggest that fullerenol C60(OH)16 represents a promising candidate as a therapeutic for Alzheimer's disease.

Ligand binding to anti-cancer target CD44 investigated by molecular simulations.

CD44 is a cell-surface glycoprotein and receptor for hyaluronan, one of the major components of the tumor extracellular matrix. There is evidence that the interaction between CD44 and hyaluronan promotes breast cancer metastasis. Recently, the molecule F-19848A was shown to inhibit hyaluronan binding to receptor CD44 in a cell-based assay. In this study, we investigated the mechanism and energetics of F-19848A binding to CD44 using molecular simulation. Using the molecular mechanics/Poisson Boltzmann surface area (MM-PBSA) method, we obtained the binding free energy and inhibition constant of the complex. The van der Waals (vdW) interaction and the extended portion of F-19848A play key roles in the binding affinity. We screened natural products from a traditional Chinese medicine database to search for CD44 inhibitors. From combining pharmaceutical requirements with docking and molecular dynamics simulations, we found ten compounds that are potentially better or equal to the F-19848A ligand at binding to CD44 receptor. Therefore, we have identified new candidates of CD44 inhibitors, based on molecular simulation, which may be effective small molecules for the therapy of breast cancer.

Alzheimer's Amyloid-ß Sequesters Caspase-3 in Vitro via Its C-Terminal Tail.

Amyloid-ß (Aß), the main constituent in senile plaques found in the brain of patients with Alzheimer's disease (AD), is considered as a causative factor in AD pathogenesis. The clinical examination of the brains of patients with AD has demonstrated that caspase-3 colocalizes with senile plaques. Cellular studies have shown that Aß can induce neuronal apoptosis via caspase-3 activation. Here, we performed biochemical and in silico studies to investigate possible direct effect of Aß on caspase-3 to understand the molecular mechanism of the interaction between Aß and caspase-3. We found that Aß conformers can specifically and directly sequester caspase-3 activity in which freshly prepared Aß42 is the most potent. The inhibition is noncompetitive, and the C-terminal region of Aß plays an important role in sequestration. The binding of Aß to caspase-3 was examined by cross-linking and proteolysis and by docking and all-atom molecular dynamic simulations. Experimental and in silico results revealed that Aß42 exhibits a higher binding affinity than Aß40 and the hydrophobic C-terminal region plays a key role in the caspase-Aß interaction. Overall, our study describes a novel mechanism demonstrating that Aß sequesters caspase-3 activity via direct interaction and facilitates future therapeutic development in AD.

Picosecond infrared laser-induced all-atom nonequilibrium molecular dynamics simulation of dissociation of viruses.

Since the discovery of the plant pathogen tobacco mosaic virus as the first viral entity in the late 1800s, viruses traditionally have been mainly thought of as pathogens for disease-resistances. However, viruses have recently been exploited as nanoplatforms with applications in biomedicine and materials science. To this aim, a large majority of current methods and tools have been developed to improve the physical stability of viral particles, which may be critical to the extreme physical or chemical conditions that viruses may encounter during purification, fabrication processes, storage and use. However, considerably fewer studies are devoted to developing efficient methods to degrade or recycle such enhanced stability biomaterials. With this in mind, we carry out all-atom nonequilibrium molecular dynamics simulation, inspired by the recently developed mid-infrared free-electron laser pulse technology, to dissociate viruses. Adopting the poliovirus as a representative example, we find that the primary step in the dissociation process is due to the strong resonance between the amide I vibrational modes of the virus and the tuned laser frequencies. This process is determined by a balance between the formation and dissociation of the protein shell, reflecting the highly plasticity of the virus. Furthermore, our method should provide a feasible approach to simulate viruses, which is otherwise too expensive for conventional equilibrium all-atom simulations of such very large systems. Our work shows a proof of concept which may open a new, efficient way to cleave or to recycle virus-based materials, provide an extremely valuable tool for elucidating mechanical aspects of viruses, and may well play an important role in future fighting against virus-related diseases.

Picosecond melting of peptide nanotubes using an infrared laser: a nonequilibrium simulation study.

Self-assembled functional peptide biomaterials are emerging with a wide range of envisioned applications in the field of nanotechnology. Currently, methods and tools have been developed to control and manipulate as well as to explore new properties of self-assembled structures. However, considerably fewer studies are being devoted to developing efficient methods to degrade or recycle such extremely stable biomaterials. With this in mind, here we suggest a theoretical framework, inspired by the recent developed mid-infrared free-electron laser pulse technology, to dissociate peptide nanotubes. Adopting a diphenylalanine channel as a prototypical example, we find that the primary step in the dissociation process occurs due to the strong resonance between the carboxylate bond vibrations of the diphenylalanine peptides and the tuned laser frequencies. The effects of laser irradiation are determined by a balance between tube formation and dissociation. Our work shows a proof of concept and should provide a motivation for future experimental developments with the final aim to open a new and efficient way to cleave or to recycle bio-inspired materials.

In silico screening for potent inhibitors against the NS3/4A protease of hepatitis C virus.

Hepatitis C virus (HCV) infections are a serious viral health problem globally, causing liver cirrhosis and inflammation that can develop to hepatocellular carcinoma and death. Since the HCV NS3/4A protease complex cleaves the scissile peptide bond in the viral encoded polypeptide to release the non-structural proteins during the viral replication process, this protease is then an important target for drug design. The computer-aided drug design and screening targeted at NS3/4A protease of HCV were reviewed. In addition, using steered molecular dynamics simulations, potent inhibitors of the NS3/4A complex were searched for by screening the ZINC database based upon the hypothesis that a high rupture force indicates a high binding efficiency. Nine top-hit compounds (59500093, 59784724, 13527817, 26660256, 29482733, 25977181, 28005928, 13527826 and 13527826) were found that had the same or a greater maximum rupture force (and so assumed binding strength and inhibitory potency) than the four current drugs and so are potential candidates as anti- HCV chemotherapeutic agents. In addition, van der Waals interactions were found to be the main contribution in stabilizing the ligand- NS3/4A complex.