Molecular recognition between two or more molecular binding partners plays an essential role in many cellular processes and pathogen entry steps. It is widely accepted now that molecular recognition is mediated through non-covalent interactions (aka non-bonded interactions) like hydrogen bond (HB), salt bridge, π-system interactions like cation–π interaction, π–π stacking interaction, CH-π interaction, and XH-π interaction (X=N, S, O); and van der Waals interaction (VDW). In this dissertation, we have studied molecular recognition in three important biological systems, i.e., agonists/antagonists in G protein-coupled receptors, guanine in GTP-binding proteins, and spike proteins of novel coronavirus with the human ACE2 receptor. It is expected that the knowledge of molecular recognition in these selected proteins would help design of drugs targeting GPCRs and GTP-binding pockets, or design of neutralizing antibodies targeting the spike protein (RBD) of the novel coronavirus.In the first system, molecular recognition of ligands that binds to GPCRs, mainly A2A receptors, was studied. It has been reported that nearly 34% of FDA-approved drugs target GPCRs. Most drugs function as either an antagonist or an agonist of GPCRs. The molecular recognition of the ligands in A2A receptors represents a topic of great importance because of the significance of adenosine and cognate ligands in cellular physiology and the necessity to develop safe and effective medications for many pathophysiological issues. A total of 5 unique agonist-bound and 18 unique antagonist-bound complexes of adenosine A2A receptor were systematically studied. Interestingly, all agonists and antagonists feature a central heterocycle that is capable of hydrogen bonding and π-system interactions. For each complex, the interaction modes between an agonist/antagonist and its interacting residues were examined to identify all the non-covalent interactions necessary for molecular recognition. The interaction mode analysis revealed that three major interaction modes, i.e., hydrogen bond from N253, π-π stacking interactions from F168, and CH-π interaction from M177, L249 and I274, are responsible for the binding of the central heterocycle to the A2A receptor. Subsequently, strengths of intermolecular interactions are quantified by means of advanced quantum chemical calculations to determine the relative importance of different modes of non-covalent interactions in A2A-ligand binding. The interaction energies were calculated at the B2PLYP level of theory, in combination with Grimme’s D3BJ dispersion correction. It was found that these 5 residues (“key” residues) contribute significantly to ligand binding. The strength of π-π stacking interactions from F168 is the strongest. Residues M177, L249 and I274 offer some strong CH-π interactions. The strength of the hydrogen bond from N253 is comparable to that of CH-π interactions. These interactions from the “key” residues can be considered the major molecular determinant responsible for the molecular recognition of ligands in the A2A receptor. Our findings are in good agreement with site-directed mutagenesis data which showed that the mutation of these “key” residues to alanine leads to the total abolition of agonist or antagonist binding. In the second system, the molecular recognition in GTP-binding proteins, was studied. G proteins, including small monomeric GTPases (small G proteins) and heterotrimeric G proteins, are the family of proteins that act as molecular switches inside cells and play an essential role in cellular signal transduction. Thus, molecular recognition of GTP/GDP in G proteins is a subject of great importance for understanding the enzymatic mechanism and designing drugs targeting the GTP binding pocket.A large-scale data mining of Protein Data Bank was performed, which resulted in 302 GTP-binding protein complexes. To decipher specific interactions responsible for molecular recognition, the intermolecular interaction modes of hydrogen bond interaction, cation-π interaction, and π-π stacking interaction between guanine and its surrounding residues were systematically analyzed. A database of such interactions has been created.We discovered 6 different hydrogen bonding patterns in which the surrounding interacting residues prefer to interact with N2 and O6 of the guanine base via hydrogen bond interaction. This results in selectivity to the guanine base because these two atoms are absent in the adenine base. In 86.1 % (260 out of 302) of GTP-binding proteins, O6 is accepting an additional hydrogen bond from the main chain amino group of alanine and/or its succeeding residue from (T/G)(C/S)A sequence motif or any non-conserved residues. The N2, on the other hand, donates two hydrogens, N2H1, and N2H2, directly to the surrounding residue. In 90.1 % (272 out of 302) of complexes, N2H1 is being donated to surrounding residues, and in 21.9 % (66 out of 302) of the complexes, N2H2 is donated to the surrounding resides. High-level quantum chemical methods like B2PLYP and MP2 methods were applied to quantify the energetic contribution of each type of interaction. Overall, our analysis showed the significant contribution of cation-π interaction (between guanine and a positively charged residue) to the binding of guanine. In 86.4% of complexes (261 out of 302), at least one cation-π interaction does exist between the guanine base and positively charged side chains. The presence of π-π stacking interactions is less than cation-π interactions. In 54.3 % of complexes (164 out of 302), π-π stacking interactions exist between the guanine base and aromatic side chains. It was observed that the strength of the intermolecular cation-π interaction depends upon the combination of three factors, i.e., distance, the extent of overlap between the side chain of positively charged residue and the guanine ring, and multiple modes of interactions (cation-π and hydrogen bond). The strength of the intermolecular π-π stacking interaction is dependent upon a combination of three factors, i.e., distance, the angle between the aromatic rings, and the extent of π-π stacking.In order to get a detailed picture of molecular recognition in a whole GTP-binding protein, a complete analysis of the entire binding pocket was carried out on a representative GTP-binding protein, the p21-ras protein. The resulting interaction energies showed the most significant contributor to guanine binding energy is cation-π interaction which contributes almost half of the guanine binding energy, -18.2 kcal/mol (49.3 %). The remaining contribution is from hydrogen bond, -15.7 kcal/mol (42.5%), and π-π stacking interactions (8.2 %). The frequent presence of cation-π interaction (86.4% of complexes) and its significant energetic contribution to the binding of guanine suggest that cation-π interaction is one of the main molecular determinants for guanine binding in GTP-binding proteins besides hydrogen bonding.The final system studied was a protein-protein complex consisting of the spike protein of novel coronavirus (SARS-CoV-2) and the human ACE2 receptor. The development of therapies that target spike protein and the prevention of viral entry depend heavily on our understanding of molecular recognition between the spike protein and ACE2 or between spike and neutralizing antibodies. Binding of the spike protein with ACE2 initiates the viral entry, and the receptor binding domain (RBD) of the spike protein is directly involved. We identified all the residues at the binding interface of RBD, and their modes of interactions to unravel the molecular determinants responsible for recognition of spike protein. 13 hydrogen bonds, 2 salt bridges, 14 CH-π, 3 π-π, 1 cation-π, 1 NH-π and 1 S-π interaction are found at the binding interface. Our high-level quantum chemical analysis revealed three hot spot areas consisting of nine strong-binding residues, and one isolated strong-binding residue. Strong-binding residues K417, Y449, L455, F456, F486, N487, Y489, T500, N501, and Y505 possess interaction energies of greater than -4.0 kcal/mol. Remarkably, some of these strong-binding residues are among the epitope residues of several potent neutralizing antibodies, indicating that these residues also play a crucial role in strong antibody binding. The MM/PBSA method of binding free energy calculation is close to the experimental value, and it can be an efficient method to calculate the binding free energy of new variants or antibodies for which the complex structure is available from experiment or homology modeling.The impact of the N501Y mutation on the binding affinity between the spike protein and the ACE2 was investigated in another protein-protein system. The N501Y mutation occurs in several highly contagious alpha-, beta-, gamma-, and omicron-strains of SARS-CoV-2, and has been linked to enhanced transmission and immune evasion. There is an urgent need to understand this mutation at the molecular level so that effective vaccines and antibody therapies against quickly evolving variants can be designed. We performed a comparative examination of the molecular recognition of the wild-type and mutant spike proteins by the human cell surface receptor ACE2. It was found that the N501Y mutation results in a gain in new interaction modes and a gain in interaction energy by almost -2.0 kcal/mol. In addition, molecular dynamics simulation showed that Y501 is surrounded by a network of water molecules, facilitating the formation of a water mediated hydrogen bond with the hydroxyl group of Y501. It is concluded that the increased binding affinity in the N501Y mutants may be partially responsible for the greater transmissibility of the VOCs with the N501Y mutation.